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ABSTRACT 

Imaging  techniques  have  been  used  for  a  long  time  in  many  applications  such  as 
nondestructive  testing,  medical  diagnosis,  civil  infrastructure  assessment,  and  homeland  security. 
Advances  in  signal  and  image  processing  techniques  have  made  considerable  progress  in  a  way 
that  the  multi-dimensional  images  of  the  target  object  are  successfully  reconstructed  from 
electromagnetic  (EM)  scattering;  the  shape  of  the  object  is  clearly  identified  by  visual  inspection. 

The  dissertation  explores  the  diffraction  tomographic  imaging  of  a  target  object  using 
ultra-wideband  (UWB)  white  Gaussian  noise  (WGN)  waveforms  in  radar  system.  The  principles 
and  validation  of  UWB  noise  radar  system  for  tomography  application  are  presented  including 
the  theoretical  analysis,  numerical  simulations  and  experimental  results. 

Using  WGN  waveform  as  the  transmit  signal  ensures  the  low  probability  of  intercept 
(LPI)  since  the  transmitted  noise  waveform  is  constantly  varying  and  never  repeats  exactly. 
Ideally,  such  a  signal  is  a  stochastic  process  with  a  flat  spectral  density  over  a  wide  frequency 
range;  however,  the  spectral  density  of  WGN  cannot  be  uniform  and  consistent  for  all  frequencies 
in  reality.  In  order  to  obtain  a  flat  spectral  density  over  the  desired  frequency  range,  a  method  of 
averaging  multiple  independent  and  identically  distributed  (iid)  noise  waveforms  is  proposed. 
From  the  simulation  results,  tomographic  images  of  various  target  scenarios  are  successfully 
reconstructed  with  random  noise  waveforms. 

In  addition,  the  hardware  system  of  bistatic  UWB  noise  radar  for  tomographic  imaging 
application  is  designed  and  implemented.  The  tomographic  image  processing  technique  for  large 
sets  of  collected  data  samples  are  also  proposed.  Both  metallic  and  dielectric  target  objects  are 
correctly  imaged  based  on  measured  scattering  data,  and  the  experimental  results  compare  well 
with  the  numerical  simulation  results. 
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Chapter  1 
Introduction 

1.1  Motivation  and  Background 

Advances  in  signal  and  imaging  processing  techniques  of  radar  systems  have  progressed 
so  that  multi-dimensional  representations  of  the  target  object  can  be  obtained  from 
electromagnetic  (EM)  scattering.  One  such  imaging  technique  is  tomography  [1].  The 
fundamental  idea  of  tomography  is  to  solve  an  inverse  scattering  problem  governed  by  wave 
equations  [2].  While  prominently  found  in  the  medical  and  industrial  scanning  areas,  tomography 
has  more  recently  found  applications  in  radar  systems.  For  example,  tomography  has  proven  to 
be  a  vital  technique  for  the  imaging  of  target  objects  via  microwave  interaction  with  ground 
penetrating  radar  (GPR)  [3]. 

In  radar  imaging  systems,  images  of  a  target  are  reconstructed  based  on  range  profiles  of 
the  distributed  sensor  elements.  When  a  transmitter  radiates  a  waveform,  spatially  distributed 
receivers  collect  samples  of  the  scattered  field  due  to  the  material  properties  of  the  target.  In  the 
next  iteration,  a  different  transmitter  is  activated,  and  the  scattered  field  collecting  process  is 
repeated.  Finally  all  collected  scattered  field  data  is  relayed  for  signal  processing  and  subsequent 
imaging.  The  optimal  waveform  and  bandwidth  enhance  the  performance  of  the  radar  system  [4]. 

While  most  prior  research  has  been  accomplished  using  coherent  electromagnetic  waves 
in  the  microwave  and  millimeter-wave  regimes,  random  field  cross-correlation  is  a  new 
promising  technique  which  bypasses  shortcomings  of  usual  active  methods.  While  coherent 
source  environments  can  generate  anomalies  in  the  correlation  function,  random  uniformly 
distributed  sources  can  cause  such  ghosts  in  the  image  to  disappear.  The  idea  is  to  use  time 
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correlations  of  noise  for  recovering  the  Green’s  function  and  from  it,  the  impulse  response.  In  the 
case  of  a  completely  random  waveform,  the  cross-correlation  of  signals  recorded  between  two 
points  converges  to  the  complete  Green’s  function  of  the  medium,  including  all  reflection, 
scattering  and  propagation  modes.  Using  an  array  of  probes  located  around  the  region  to  be 
imaged,  the  detected  noise  correlation  can  be  shown  to  be  a  function  of  the  spatial  overlap  of  the 
electromagnetic  fields  at  the  probes  and  the  spatial  distribution  of  the  dielectric  constant  and/or 
conductivity  of  the  sample.  A  three-dimensional  conductivity  map  has  been  successfully 
generated  using  passive  thermal  noise  emissions  to  image  breast  tumors  and  buried  isolated 
scatterers  [5]. 

In  this  dissertation,  the  main  research  is  focused  on  tomographic  imaging  for  UWB  noise 
radar,  which  is  considered  as  a  promising  technique  for  high-resolution  imaging  without  alerting 
the  presence  of  the  radar  system.  Theoretical  analysis,  numerical  simulations  and  practical 
hardware  implementation  of  UWB  noise  radar  tomography  are  presented  in  detail. 


1.2  Contributions  of  the  Dissertation 

The  research  that  I  have  performed  during  my  Ph.D.  studies  at  the  Pennsylvania  State 
University  has  led  to  several  developments  in  the  field  of  UWB  noise  radar  tomography. 

•  White  Gaussian  noise  waveform  modeling  and  the  practical  implementation  of 
bandlimited  UWB  noise  waveforms  for  tomographic  imaging  application  are 
proposed.  More  specifically,  multiple  transmissions  of  independent  and  identically 
distributed  (iid)  WGN  noise  waveforms  are  suggested  to  bypass  the  shortcomings  of 
using  UWB  random  noise  waveform. 
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•  The  principle  of  tomography  based  noise  radar  imaging  is  proposed.  The 
tomographic  image  of  the  target  object  is  successfully  reconstructed  using  the  Fourier 
diffraction  theorem. 

•  Numerical  simulations  are  performed  to  validate  the  theory  of  Fourier  diffraction  for 
various  target  scenarios.  More  specifically,  multiple  cylindrical  targets  in  different 
sizes  are  correctly  imaged  with  a  linearly  positioned  receiver  array  under  various 
SNR  settings. 

•  The  image  quality  of  the  reconstructed  tomographic  images  is  quantitatively  analyzed. 
Furthermore,  the  image  quality  measure  for  various  SNR  settings  can  be  used  as  a 
guideline  for  the  practical  implementation  of  UWB  noise  radar  system  for 
tomographic  imaging  applications. 

•  Theoretical  analysis  of  noise  radar  tomography  is  proposed  for  sequentially  rotated 
target  objects.  The  numerical  simulations  of  a  bistatic  imaging  radar  system  are 
performed  to  obtain  the  complete  cross-sectional  images  of  various  shaped  target 
objects. 

•  Hardware  implementation  of  UWB  noise  radar  system  for  tomography  application  is 
proposed.  The  reconstructed  tomographic  images  based  on  experimental  results  can 
be  compared  with  those  obtained  based  on  numerical  simulation  results,  showing  a 
good  agreement  between  simulation  and  experiment  results  for  both  metallic  and 
dielectric  target  object  cases. 


1.3  Overview 

Chapter  2  presents  the  theoretical  foundations  of  diffraction  tomography  using  a  white 
Gaussian  noise  (WGN)  waveform  as  a  transmitted  signal  for  the  cross-section  images  of  multiple 
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metallic  target  objects  inside  a  rectangular  geometry  for  the  region  of  interest.  The  empirical 
solution  to  bypass  the  shortcomings  of  using  UWB  random  noise  waveform  and  the  image  quality 
measures  of  the  reconstructed  images  under  various  SNR  environments  are  discussed  as  well. 

Chapter  3  focuses  on  the  principle  and  numerical  simulation  results  of  tomographic 
imaging  of  various  target  scenarios  using  UWB  noise  radar.  Theoretical  analysis  of  the  image 
reconstruction  for  a  bistatic  imaging  radar  system  with  a  two-dimensional  rotating  cylinder  using 
the  Fourier  diffraction  theorem  under  the  assumption  of  plane  wave  illumination  is  presented. 
Also,  the  numerical  simulation  results  of  diffraction  tomography  for  UWB  noise  radar  are 
presented. 

Chapter  4  shows  the  hardware  implementation  and  experimental  investigation  of  UWB 
noise  radar.  Design  considerations  and  performance  for  UWB  noise  radar  using  WGN 
waveforms  are  presented.  The  effective  tomographic  image  processing  technique  for  large  sets  of 
collected  data  samples  is  proposed  as  well.  The  reconstructed  tomographic  images  based  on 
experiments  can  be  compared  with  those  obtained  based  on  numerical  simulation  results  for 
various  target  objects. 

Chapter  5  discusses  about  conclusions  and  suggestions  for  future  research  work  regarding 
UWB  noise  radar  tomography. 
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Chapter  2 

Diffraction  Tomography  using  UWB  Random  Noise  Waveforms 

2.1  Overview 

Research  on  the  use  of  random  or  pseudorandom  noise  transmit  signals  in  radar  has  been 
conducted  since  the  1950s  [1,2].  Noise  radar  has  been  considered  a  promising  technique  for  the 
covert  identification  of  target  objects  due  to  several  advantages,  such  as  excellent  electronic 
countermeasure  (ECM),  low  probability  of  detection  (LPD),  low  probability  of  interception  (LPI) 
features,  and  relatively  simple  hardware  architectures  [3-5].  Also,  advances  in  signal  and 
imaging  processing  techniques  in  radar  systems  have  progressed  so  that  multidimensional 
representations  of  the  target  object  can  be  obtained  [6]. 

In  general,  radar  imaging  tends  to  be  formulated  in  the  time  domain  to  exploit  efficient 
back-projection  algorithms,  generate  accurate  shape  features  of  the  target  object,  and  provide 
location  data  [7].  For  multistatic  radar  systems,  the  images  of  a  target  are  reconstructed  based  on 
range  profiles  obtained  from  the  distributed  sensor  elements.  When  a  transmitter  radiates  a 
waveform,  spatially  distributed  receivers  collect  samples  of  the  scattered  field  which  are  related 
to  the  electrical  parameters  of  the  target  object.  For  the  next  iteration,  a  different  transmitter  is 
activated,  and  the  scattered  field  collection  process  is  repeated.  Finally,  all  collected  scattered 
field  data  are  relayed  for  signal  processing  and  subsequent  image  formation  algorithms. 

Tomography -based  radar  imaging  algorithms  have  been  developed  based  on  the 
microwave  image  reconstruction  method  [9],  characterizing  the  material  property  profiles  of  the 
target  object  in  the  frequency  domain  and  reconstructing  specific  scattering  features  inside  the 
interrogation  medium  by  solving  the  inverse  scattering  problem.  The  capability  of  microwave 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


21 


7 


imaging  techniques  has  been  found  attractive  in  malignant  breast  cancer  detection  [10-13],  civil 
infrastructure  assessment  [14-16],  and  homeland  security  [17-19].  The  quality  of  the 
reconstructed  image  for  different  values  of  the  electrical  contrast  for  a  UWB  imaging  system  was 
investigated  and  published  for  both  low-contrast  and  high-contrast  object  cases.  For  low-contrast 
objects,  the  obtained  target  image  using  a  single  frequency  achieves  a  good  reconstruction  of  the 
electrical  contrast  that  is  almost  equivalent  to  the  one  obtained  with  the  entire  UWB  frequency 
range.  For  the  high-contrast  case,  while  the  formation  of  a  Moire  pattern  affects  the  single 
frequency  reconstruction,  this  artifact  does  not  appear  in  the  UWB  frequency  image  [20].  Thus, 
UWB  radar  tomography  is  expected  to  provide  advantages  over  the  single  or  narrow  band 
frequency  operation  in  terms  of  resolution  and  accuracy  for  any  target  object. 

This  chapter  is  organized  as  follows.  First,  the  characteristics  of  UWB  random  noise 
signal  and  the  shortcomings  of  using  such  noise  signal  as  a  radar  transmit  waveform  in 
tomographic  image  reconstruction  process  are  discussed  in  section  2.2.  The  empirical  solution  to 
bypass  the  shortcoming  of  using  UWB  random  noise  waveform  is  also  proposed.  In  section  2.3, 
the  formulations  of  the  image  reconstruction  of  two-dimensional  scattering  geometry  using  the 
Fourier  diffraction  theorem  under  the  assumption  of  plane  wave  illumination  is  presented.  The 
numerical  simulation  results  of  diffraction  tomography  imaging  for  various  arrangements  of 
circular  PEC  cylindrical  objects  are  presented  in  section  2.4.  The  image  quality  measures  of  the 
reconstructed  images  and  SNR  effects  for  multiple  transmissions  of  UWB  random  noise 
waveforms  are  discussed  in  section  2.5.  Conclusions  are  presented  in  section  2.6. 
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2.2  Analysis  of  White  Gaussian  Noise  Model 


The  main  advantage  of  transmitting  a  random  noise  waveform  is  to  covertly  detect  and 
image  a  target  without  alerting  others  about  the  presence  of  the  radar  system.  Such  LPI 
characteristics  of  the  noise  radar  are  guaranteed  because  the  transmitted  random  noise  waveform 
is  constantly  varying  and  never  repeats  itself  exactly  [20].  The  random  noise  waveform  can  be 
experimentally  generated  simply  by  amplifying  the  thermal  noise  generated  in  resistors  or  noise 
diodes  while  maintaining  relatively  flat  spectral  density  versus  frequency  [21].  Hence,  relatively 
simple  hardware  designs  can  be  achieved  for  noise  radars  compared  to  the  conventional  radar 
systems  using  complicated  signal  modulation  schemes. 

For  a  random  noise  waveform  model,  let  x\n\  be  a  discrete  time  WSS  and  ergodic  random 
process  and  a  sequence  of  independent  and  identically  distributed  (iid)  random  variable  drawn 
from  a  Gaussian  distribution,  A(0,  a^).  x\n\  defined  herein  as  white  Gaussian  noise,  that  is,  its 
probability  density  function  follows  a  Gaussian  distribution  and  its  power  spectral  density  is 
ideally  a  non-zero  constant  for  all  frequencies.  However,  the  finite  number  of  random  noise 
amplitude  samples  must  be  chosen  for  waveform  generation  for  any  numerical  simulations  and 
practical  experiments. 

Assume  that  a  sequence  of  only  /  samples  of  x\n\  is  selected  for  generating  a  white 
Gaussian  noise.  In  this  case,  the  estimate  for  the  power  spectral  density,  S^{co) ,  is  given  by  [22] 

I,{a))  =  S,{a))=  X  (2.1) 

where  r\m\  is  the  estimate  for  the  autocorrelation  sequence.  Ii{cd)  is  defined  as  the  periodogram 
estimate,  and  the  rigorous  analysis  of  the  expected  value  and  variance  of  the  periodogram 
estimate  for  any  arbitrary  co  is  described  in  [22-24].  The  expected  value  of  the  periodogram 
estimate  is  [22,  23] 
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=  z  1-7 


which  suggests  that  (of)  is  a  biased  estimator.  However,  it  is  considered  to  be  asymptotically 
unbiased  as  /  approaches  infinity.  In  this  case,  the  expected  valued  of  (of)  becomes  a  constant 


such  that 


E[/^((^))]n  S^(of)  =  c7^  as  /^oo.  (2.3; 

The  variance  of  the  periodogram  estimate  of  the  white  Gaussian  noise  waveform  formed 
by  a  sequence  of  /  samples  is  given  by  [22,  23] 


VAR[/,(^)]  =  5,(®r  1 


sin(^u/) 
/  sin(^u) 


which  is  proportional  to  the  square  of  the  power  spectrum  density  and  does  not  approach  zero  as  / 
increases.  In  order  to  decrease  the  variance  of  Ii(of) ,  the  periodogram  averaging  method  has 
been  proposed  by  Bartlett  [25].  The  average  of  K  independent  and  identically  distributed 
periodograms  on  samples  of  size  /  is  given  by 


^  /=o 


and  the  expected  value  of  the  average  with  K  iid  periodogram  estimate  is  written  as  [22] 


E[l,^ico)]  =  E[l,ico)] 


which  is  considered  to  be  asymptotically  unbiased  as  /  approaches  infinity.  Also  the  variance  of 


the  averaged  periodogram  estimate  is  given  by  [22] 
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VAR[/,J^y)]  =  lvAR[/,(®)] 
K 


(2.7) 


The  variance  of  the  averaged  periodogram  estimate  is  inversely  proportional  to  the 
number  of  iid  periodograms  K,  and  consequently  the  variance  approaches  zero  as  K  approaches 
infinity.  Equations  (2.6)  and  (2.7)  are  used  to  conclude  that  the  expected  value  remains 
unchanged,  but  only  the  variance  of  white  Gaussian  noise  decreases  for  averaging  K  iid 
periodogram  estimates.  Increasing  the  number  of  K  in  averaging  periodogram  estimate  truly 
flattens  the  spectral  density,  and  the  successful  tomographic  image  can  be  achieved  by 
transmitting  K  multiple  random  noise  waveforms  with  a  large  sequence  size  /.  For  the  numerical 
simulations  performed  in  this  chapter,  a  total  of  10  iid  random  noise  waveforms  are  transmitted, 
and  each  iid  noise  waveform  is  generated  with  500  random  amplitude  samples  drawn  from  A^(0, 
(7^).  The  tomographic  image  is  formed  based  on  the  dataset  from  41  discrete  frequencies  chosen 
uniformly  within  X-band  from  8  GHz  to  10  GHz  in  steps  of  50  MHz. 


2.3  Formulation  of  the  Scattering  Properties  for  Cylindrical  PEC  Object 

In  this  section,  the  scattering  properties  for  two-dimensional  cylindrical  impenetrable 
conducting  object  in  the  bistatic  scattering  arrangement  are  discussed,  and  the  Fourier  diffraction 
tomography  algorithm  is  applied  to  reconstruct  the  image  of  the  object  based  on  the  bistatic 
scattering  properties.  The  Fourier  diffraction  theorem  has  been  extensively  applied  in  the  area  of 
acoustical  imaging  [26-28].  The  goal  of  diffraction  tomography  is  to  reconstruct  the  properties  of 
a  slice  of  an  object  from  the  scattered  field.  For  planar  geometry,  an  object  is  illuminated  with  a 
plane  wave,  and  the  scattered  fields  are  calculated  or  measured  over  a  straight  line  parallel  to  the 
incident  plane  wave.  The  mathematical  formulation  and  proof  of  validity  of  the  Fourier 
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diffraction  theorem  shown  in  Figure  2-1  are  not  presented  in  this  section  since  they  were  already 
stated  and  published  [29]. 


space  D«>inaiii 


Frequency  Domnin 


Figure  2-1.  The  Fourier  diffraction  theorem  relates  the  Fourier  transform  of  a  diffracted 
projection  to  the  Fourier  transform  of  the  object  along  a  semicircular  arc.  An  arbitrary  object  is 
illuminated  by  a  plane  wave  propagating  along  the  unit  vector,  and  the  coordinate  system  is 
rotated. 


Two-dimensional  backward  scattering  geometry  for  a  cylindrical  conducting  object  is 
shown  in  Figure  2-2.  Prior  research  on  microwave  imaging  has  already  proven  that  the  image 
reconstructed  in  the  backward  scattering  case  is  better  than  that  obtained  in  the  forward  scattering 
case,  based  on  numerical  results  [30]. 

As  shown  in  Figure  2-2,  a  single  cylindrical  perfect  electric  conductor  (PEC)  object  is 
located  at  the  center  of  the  simulation  scene.  The  height  of  the  cylindrical  object  is  assumed  to  be 
infinitely  long  along  the  z-axis  and  the  incident  z-polarized  plane  wave  is  illuminated  from  -x 
direction.  The  receiver  spacing.  Ay,  is  calculated  to  avoid  aliasing  effect,  and  is  given  by 


Ay< 


2 


(2.8) 
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where  2min  is  the  minimum  wavelength.  A  linear  receiving  array  is  located  x  =  -d  from  the 
center  of  the  cylinder  object,  collecting  scattered  field  for  reconstruction  of  the  object  image.  The 
frequency  sampling  interval,  A/  is  also  considered  to  ensure  an  image  that  is  free  of  false  aliases, 
and  is  given  by 

A/<A  (2.9) 

2-r 

where  c  is  the  speed  of  wave  propagation  in  free  space,  and  r  is  the  radius  of  the  cylindrical  PEC 
object. 


Incident 
plane  wave 


E 

I 


y 


x  =  -d  y 


Figure  2-2.  Two-dimensional  backward  scattering  geometry  for  a  cylindrical  conducting  object. 
Red  dots  and  green  circle  represent  a  linear  receiving  array  and  the  PEC  cylinder,  respectively. 


As  shown  in  the  previous  section,  the  frequency  response  of  the  transmitted  white 
Gaussian  noise  waveform  shows  that  the  field  amplitude  is  a  constant  non-zero  amplitude  value 
for  all  possible  frequencies.  Thus,  the  z-polarized  incident  plane  wave  of  a  single  transmitted 
white  Gaussian  noise  waveform  for  N  discrete  frequencies  takes  the  form 


(F)  =  z  ■  + . . .  +  ) 


N 


n=l 


(2.10) 
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where  k^=co^ !  c  is  the  wavenumber,  and  is  the  field  amplitude  of  the  transmitted  white 
Gaussian  noise  waveform  at  each  discrete  frequency  of  interest.  The  field  amplitude  of  the  noise 
waveform  becomes  the  non-zero  constant  value  as  the  sequence  of  the  noise  waveform 
approaches  infinity. 

The  analysis  is  started  from  a  single  frequency,  and  the  process  for  N  multiple 
frequencies  are  developed  by  summing  up  the  analysis  results  for  N.  Also  the  entire  analysis 
must  be  repeated  for  K  times  when  K  multiple  iid  noise  waveforms  are  transmitted.  The  z- 
polarized  incident  waveform  is  defined  as 


71 


(2.11) 


where  rj  =  yMo/^o  is  the  intrinsic  impedance  in  free  space.  If  the  object  consists  of  a  material 
having  a  certain  dielectric  constant  value,  the  equivalent  electric  current  distribution,  Jeq  ,  is 
calculated  for  the  scattered  field.  The  object  is  defined  as  PEC  so  that  the  scattered  field  observed 
at  the  linear  receiving  array  in  the  y-direction  located  3Xx  =  -d  is  calculated  by  applying  the 
physical  optics  approximation,  and  given  by 


,{K,x  =  -d,y)  =  -j(oju^  jj^^-G{r-r')dr' 


-j\vj  ( 2«(r ')  X  (r'))  ■G{r-r')dr' 


(2.12) 


where  S  is  the  boundary  of  scatterer,  n(r)  is  the  outward  unit  normal  vector  to  S,  and  G(r  -  r') 
is  the  Green’s  function  for  two-dimensional  geometry  defined  as 

G{r-r')  =  -jH,^^^{k„(r-r'))  (2.13) 

where  is  the  zeroth-order  Hankel  function  of  the  second  kind. 
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Similarly,  for  the  case  of  discrete  scattering  objects,  the  scattered  field  of  p  objects  over  a 
linear  array  in  the  y-direction  located  ?iix  =  -d  is  expressed  as 

^sca,.n.uMple(K’X  =  -d ,  y)  =  J^E^ca,,single(P)-  (2-14) 

P 

On  assuming  the  polarization  in  the  z-direction,  the  general  form  of  scattered  field,  , 
obtained  by  the  receivers  ?iix  =  -d  becomes 


^scat  (k„,x  =  -d,  y)  =  z-  {k„,x  =  -d,  y) 

=  -JKE„  J| '■'G (r  -  r') d^r' 


(2.15) 


where 


Oscat, single  (f')  =  -2«(F')  •  X^(5(F'))  (2.16) 

is  defined  as  the  scattering  object  function  of  a  single  PEC  object  which  is  related  to  the  object 
shape,  and  <J(5’(F'))  is  a  Dirac  delta  function  defined  as 


^(VF')) 


=  0  as  F'  €  5' 
^  0  elsewhere 


(2.17) 


Similarly,  the  scattering  function  shown  in  (2.16)  for  discrete  scattering  objects  can  be  written  as 


o 


scat , multiple 


ir')  =  YjOsca,,singleiP)  ■ 


(2.18) 


By  using  the  plane  wave  expansion  of  the  Green’s  function  [31],  the  one-dimensional 
Fourier  transform  of  ,  defined  in  (2.15),  in  y-direction  can  be  written  as 


u^cai  {K,x=-d,k,)=^  {-r-K,ky) 


(2.19) 


where 


r  =  \ 


\p. 


y 

2  FT 


as 


k\^k 


as,k^>K 


(2.20) 


If  is  defined  as 
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K  =  -r-K^  (2.21) 

the  Fourier  transform  of  the  two-dimensional  scattering  object  function,  ,  defined  in 

(2.19)  is  given  by 

.  (2.22) 

In  this  case,  the  arguments  of  are  related  by 

iK  ^Kf  ^  ■  (2.23) 

Equations  (2.19)  and  (2.23)  show  that  as  a  two-dimensional  scattering  object  is 
illuminated  by  a  plane  wave  at  the  frequency  co^ ,  the  one-dimensional  Fourier  transform  of  the 
scattered  field  yields  a  semicircle  centered  at  (0,-^^)  with  radius  k^  in  the  two-dimensional 
Fourier  space,  which  relates  to  the  Fourier  diffraction  theorem  depicted  in  Figure  2-1.  By  using  N 
number  of  frequencies,  the  radius  in  the  two-dimensional  Fourier  space  changes,  which  enhances 
the  resolution  and  accuracy  of  the  image.  Also  the  variance  of  the  frequency  response  is  reduced 
by  taking  average  of  the  frequency  responses  of  K  multiple  transmissions  of  iid  UWB  noise 
waveforms  [32,  33]. 


2.4  Numerical  Simulation  Results 

2.4.1  Diffraction  tomography  with  a  single  transmitted  WGN  waveform  for  a  single 
scattering  object 

As  shown  in  Figure  2-3,  two-dimensional  backward  scattering  geometry  with  a  single 
cylindrical  conducting  object  is  simulated  with  two  band-limited  iid  UWB  WGN  waveforms. 
The  cylindrical  PEC  object  with  a  radius  of  15  cm  is  located  at  the  origin,  and  the  cylinder  is 
assumed  to  be  infinitely  long  along  the  z-axis.  The  scattered  field  is  uniformly  sampled  at 
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receiving  array  Rxl  through  RxlOl  with  frequency  swept  within  X-band  from  8  GHz  to  10  GHz 
in  41  steps  of  50  MHz.  The  locations  of  Rxl  and  RxlOl  are  at  (-90  cm,  -75  cm)  and  (-90  cm, 
75  cm),  respectively,  and  the  receiver  spacing.  Ay,  is  set  to  1.5  cm  based  on  (2.8)  when/nax  is  10 
GHz.  The  maximum  frequency  stepping  interval.  A/  is  also  calculated  using  (2.9),  which  yields 
a  value  of  1  GHz;  however,  50-MHz  frequency  stepping  interval  enhances  the  quality  of 
tomographic  image  compared  to  the  maximum  frequency  stepping  interval  of  1  GHz. 


z 


Figure  2-3.  Two-dimensional  backward  scattering  simulation  geometry  for  a  cylindrical 
conducting  object.  A  single  PEC  cylinder  with  a  radius  of  15  cm  is  located  at  (0  cm,  0  cm),  and  a 
linear  receiving  array  is  located  90  cm  away  from  the  origin  in  -x  direction. 


Figure  2-4(a)  and  Figure  2-5(a)  show  the  two  iid  UWB  WGN  waveforms.  Each 
waveform  is  generated  with  500  random  amplitude  samples  drawn  from  A^(0,  cr^),  and  they  are 
transmitted  after  8-10  GHz  bandpass  filtering  operation  to  reconstruct  tomographic  images  of  the 
cylindrical  PEC  object  based  on  the  scattered  field  observed  at  Rxl  through  RxlOl. 

A  block  diagram  shown  in  Figure  2-6  displays  the  tomographic  image  reconstruction 
method  using  diffraction  tomography.  The  one-dimensional  Fourier  transformed  scattered  field 
data  collected  at  the  receiving  array  Rxl  through  RxlOl  is  Fourier  transformed  into  two- 
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dimensional  object  Fourier  space  data,  0{k^,ky) ,  by  using  (2.19).  Such  Fourier  space  data  is 
two-dimensional  inverse  Fourier  transformed  to  obtain  the  scattering  object  function,  o(x,  y). 
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Input  Noise  Waveform  #1  for  Tomography  (500  amplitude  samples) 


Time  (ns) 


Figure  2-4.  (a)  The  first  UWB  WGN  waveform  generated  with  500  amplitude  samples  (/=500) 
drawn  fromA^(0,  a^).  Pulse  duration  is  2.4  ns.  (b)  The  frequency  spectrum  after  8-10  GHz 
bandpass  filtering  operation.  The  frequency  ranges  are  shown  from  DC  to  12  GHz  only. 


Input  Noise  Waveform  #2  for  Tomography  (500  amplitude  samples) 


Time  (ns) 


Figure  2-5.  (a)  The  second  UWB  WGN  waveform  generated  with  500  amplitude  samples  (/=500) 
drawn  fromA^(0,  a^).  Pulse  duration  is  2.4  ns.  (b)  The  frequency  spectrum  after  8-10  GHz 
bandpass  filtering  operation.  The  frequency  ranges  are  shown  from  DC  to  12  GHz  only. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


32 


18 


Figure  2-6.  The  image  reeonstmetion  method  using  diffraetion  tomography. 


Fourier  spaee  data  and  the  tomographie  images  of  the  single  eylindrieal  eondueting  objeet 
using  the  first  UWB  WGN  waveform  shown  in  Figure  2-4(a),  are  ealeulated  and  displayed  in 
Figure  2-7(a)  and  Figure  2-7,  respeetively.  The  tomographie  image  of  the  single  PEC  eylinder 
appears  to  be  sueeessfully  reeonstrueted  using  the  first  UWB  WGN  transmitted  waveform. 


-  -  -  -U.u  -U.H-  U 

[rad/m]  x[m] 


Figure  2-7.  (a)  The  normalized  magnitude  of  Fourier  spaee  data  of  the  single  eylindrieal 

eondueting  objeet,  and  (b)  the  normalized  tomographie  image  using  the  first  UWB  WGN 
waveform  shown  in  Figure  2-4(a). 


Similarly,  both  Fourier  spaee  data  and  the  tomographie  images  of  the  single  PEC  eylinder 
using  the  seeond  UWB  WGN  waveform  shown  in  Figure  2-5(a)  are  also  shown  in  Figure  2-8(a) 
and  Figure  2-8(b),  respeetively.  In  this  ease,  the  tomographie  image  is  signifieantly  affeeted  by 
the  unexpeeted  noteh  observed  in  the  frequeney  speetrum  at  9.1  GHz  in  Figure  2-5(b),  so  that  the 
tomographie  image  of  the  objeet  eannot  be  aehieved  eorreetly  eompared  to  the  previous  ease 
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where  the  first  WGN  waveform  is  transmitted.  Based  on  the  simulation  results  with  two  band- 
limited  iid  UWB  WGN  waveforms,  a  single  transmission  of  WGN  waveform  may  or  may  not  be 
suffieient  to  reeonstruet  a  sueeessful  tomographie  image  of  the  objeet  with  diffraetion 
tomography  algorithm  due  to  the  undesired  and  unexpeeted  notehes  in  the  speetral  density  of  the 
praetieal  UWB  WGN  transmitted  waveforms. 


Figure  2-8.  (a)  The  normalized  magnitude  of  Fourier  spaee  data  of  the  single  eylindrieal 

eondueting  objeet,  and  (b)  the  normalized  tomographie  image  using  the  first  UWB  WGN 
waveform  shown  in  Figure  2-5(a). 


2.4.2  Diffraction  tomography  with  multiple  transmitted  iid  WGN  waveforms  for  a  single 
scattering  object 

Figure  2-9  shows  K  iid  UWB  WGN  waveforms  being  transmitted  to  reeonstruet  a  final 
image  of  the  objeet  in  order  to  bypass  a  shorteoming  of  the  single  transmission  of  WGN 
waveform  displayed  in  Figure  2-8(b).  The  proposed  imaging  method  with  multiple  iid  WGN 
transmitted  waveforms  is  established  based  on  the  method  of  the  WGN  periodogram  averaging 
deseribed  in  seetion  2.2. 

As  shown  in  Figure  2-9,  the  final  tomographie  image  of  the  seattered  objeet  is 
reeonstrueted  via  sum  and  average  of  all  K  diserete  images  for  K-ih  band-limited  iid  UWB  WGN 
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transmitted  waveforms.  For  seattered  field  and  diffraetion  tomography  simulations  with 
multiple  iid  WGN  waveforms,  10  iid  UWB  WGN  waveforms  over  a  frequeney  range  from  8  GHz 
to  10  GHz  are  generated  with  500  random  amplitude  samples  drawn  from  A^(0,  cr^),  and 
transmitted  for  baekward  seattering  field  data. 


Figure  2-9.  The  image  reeonstruetion  method  with  K  multiple  iid  WGN  transmitted  waveforms 
using  diffraetion  tomography. 


Figure  2- 10(a),  Figure  2- 10(b),  Figure  2- 10(e),  and  Figure  2- 10(d)  display  the  four  final 
tomographie  images  when  one,  three,  seven  and  all  ten  diserete  images  are  summed  and  averaged, 
respeetively.  Sueeessful  tomographie  imaging  of  the  target  is  aehieved  after  averaging  all  ten 
images  by  visual  inspeetion  of  the  formed  images  as  shown  in  Figure  2- 10(d),  and  inereasing  the 
number  of  transmissions  of  the  iid  UWB  WGN  waveform  tends  to  enhanee  the  quality  of  final 
tomographie  image  of  the  objeet  by  redueing  the  varianee  of  the  speetral  response  of  WGN. 
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Figure  2-10.  The  final  tomographie  image  of  a  single  PEC  eylinder  loeated  at  (0,0)  after 
summing  and  averaging  proeess  with  the  (a)  one  transmitted  WGN  waveform  image,  (b)  three 
transmitted  WGN  waveform  images,  (e)  seven  transmitted  WGN  waveform  images,  and  (d)  all 
ten  transmitted  WGN  waveform  images.  The  eolorbar  indieates  the  normalized  magnitude  of 
seattering  objeet  funetion,  o(x,  y). 


2.4.3  Diffraction  tomography  with  multiple  transmitted  iid  WGN  waveforms  for  two 
symmetrically  distributed  scattering  objects 

Based  on  the  simulation  result  of  the  single  PEC  eylinder  ease  with  multiple  transmitted 
noise  waveforms,  transmitting  multiple  iid  UWB  WGN  waveforms  delivers  the  aeeeptable 
quality  of  image  of  the  objeet  for  diffraetion  tomography. 
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Incident  noise  waveform 


Figure  2-11.  Two-dimensional  backward  scattering  simulation  geometry  for  two  symmetrically 
distributed  cylindrical  conducting  objects.  Two  PEC  cylinders  with  radii  of  15  cm  are  located  at 
(0  cm,  -40  cm)  and  (0  cm,  40  cm),  and  a  linear  receiving  array  is  located  90  cm  away  from  the 
origin  in  -x  direction. 


Two-dimensional  backward  scattering  geometry  for  two  cylindrical  conducting  objects, 
which  is  shown  in  Figure  2-1 1,  is  also  simulated  with  10  iid  UWB  WGN  waveforms.  Two 
cylindrical  PEC  objects  with  radii  of  15  cm  are  located  at  (0  cm,  -40  cm)  and  (0  cm,  40  cm)  such 
that  they  are  symmetrically  positioned  with  respect  to  the  x-axis.  The  distance  from  the  center  of 
each  PEC  cylinder  to  the  x-axis,  a,  is  40  cm,  and  they  are  also  assumed  to  be  infinitely  long  along 
the  z-axis.  The  coordinates  of  all  101  receivers,  and  the  receiver  spacing.  Ay,  remained  the  same 
as  shown  in  Figure  2-3.  The  scattered  field  is  uniformly  sampled  at  Rxl  through  RxlOl  with 
same  frequency  swept  within  X-band  over  8  GHz  to  10  GHz  in  41  steps  as  well. 

For  scattered  field  and  diffraction  tomography  simulations  for  multiple  scattering  objects 
with  multiple  iid  WGN  waveforms,  the  same  10  iid  UWB  WGN  waveforms,  which  are  used  for 
the  simulation  with  a  single  PEC  object  in  the  previous  section,  are  transmitted  in  the  same 
sequence.  The  scattered  field  due  to  two  symmetrically  positioned  PEC  objects  is  calculated 
based  on  the  equations  derived  in  section  2.3,  and  the  final  tomographic  image  of  objects  is 
formed  via  the  proposed  image  reconstruction  method  shown  in  Figure  2-9. 
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Figure  2- 12(a),  Figure  2- 12(b),  Figure  2- 12(e),  and  Figure  2- 12(d)  display  four  final 
tomographie  images  for  two  symmetrieally  distributed  PEC  eylinders  ease  when  one,  three,  seven 
and  all  ten  diserete  images  are  summed  and  averaged,  respeetively.  As  shown  in  Figure  2- 12(d), 
tomographie  image  of  the  two  PEC  eylinders  is  sueeessfully  reeonstrueted  with  multiple 
transmitted  noise  waveforms  as  expeeted.  However,  the  mutual  eoupling  effeets  due  to  the 
multiple  PEC  objeets  are  also  imaged  in  all  tomographie  images. 


Figure  2-12.  The  final  tomographie  image  of  two  symmetrieally  distributed  PEC  eylinders 
loeated  at  (0  em,  -40  em)  and  (0  em,40  em)  after  summing  and  averaging  proeess  with  the  (a) 
one  transmitted  WGN  waveform  image,  (b)  three  transmitted  WGN  waveform  images,  (e)  seven 
transmitted  WGN  waveform  images,  and  (d)  all  ten  transmitted  WGN  waveform  images. 
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2.4.4  Diffraction  tomography  with  multiple  transmitted  iid  WGN  waveforms  for  two 
randomly  distributed  scattering  objects 

Two-dimensional  backward  scattering  geometry  for  two  randomly  distributed  cylindrical 
conducting  objects  is  shown  in  Figure  2-13.  The  radii  of  both  cylinders  are  15  cm,  and  the 
positions  of  two  cylinders  are  (-7.5  cm,  22.5  cm)  and  (35  cm,  -55  cm).  The  general 
configuration  of  simulation  geometry,  such  as  number  of  transmitted  UWB  WGN  waveforms, 
frequency  swept  ranges,  coordinates  of  the  linear  receiving  array,  and  the  receiver  spacing,  is 
identical  to  the  previous  cases. 


Figure  2-13.  Two-dimensional  backward  scattering  simulation  geometry  for  two  randomly 
distributed  cylindrical  conducting  objects.  PEC  cylinders  with  radii  of  15  cm  are  located  at 
(-7.5  cm,  22.5  cm)  and  (35  cm,  -55  cm). 


Four  final  tomographic  images  for  two  randomly  distributed  PEC  cylinders  are  displayed 
in  Figure  2-14(a),  Figure  2-14(b),  Figure  2-14(c),  and  Figure  2-14(d)  when  one,  three,  seven  and 
all  ten  discrete  images  are  summed  and  averaged,  respectively.  Again,  the  mutual  coupling  effects 
are  shown  in  the  tomographic  images  as  expected.  The  reconstructed  images  are  shown  to  be  in 
good  agreement  with  the  simulation  geometry  given  in  Figure  2-13. 
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Figure  2-14.  The  final  tomographie  image  of  two  randomly  distributed  PEC  eylinders  loeated  at 
(-7.5  em,  22.5  em)  and  (35  em,  -55  em)  after  summing  and  averaging  proeess  with  the  (a)  one 
transmitted  WGN  waveform  image,  (b)  three  transmitted  WGN  waveform  images,  (e)  seven 
transmitted  WGN  waveform  images,  and  (d)  all  ten  transmitted  WGN  waveform  images. 


2.4.5  Diffraction  tomography  with  multiple  transmitted  iid  WGN  waveforms  for  three 
randomly  distributed  scattering  objects  in  different  sizes 

As  shown  in  Figure  2-15,  three  PEC  eylinders  with  radii  of  7.5  em,  10  em,  and  15  em  are 
loeated  at  (-22.5  em,  -20  em),  (10  cm,  50  cm),  and  (25  cm,  -40  cm),  respectively.  Again,  the 
general  configurations  of  simulation  geometry  remain  unchanged  as  defined  in  previous  cases. 
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Figure  2-15.  Two-dimensional  backward  scattering  simulation  geometry  for  three  randomly 
distributed  cylindrical  conducting  objects  in  different  sizes.  Three  PEC  cylinders  with  radii  of 
7.5  cm,  10  cm,  and  15  cm  are  located  at  (-22.5  cm,  -20  cm),  (10  cm,  50  cm),  and  (25  cm,  -40  cm), 
respectively. 


Figure  2- 16(a),  Figure  2- 16(b),  Figure  2- 16(c),  and  Figure  2- 16(d)  show  four  final 
tomographic  images  for  the  scattering  geometry  given  in  Figure  2-15  when  one,  three,  seven  and 
all  ten  discrete  images  are  summed  and  averaged,  respectively.  As  shown  in  the  previous  cases, 
the  image  quality  of  reconstructed  tomographic  images  for  unevenly  distributed  multiple 
scattering  objects  is  affected  by  the  mutual  coupling  effects.  Successful  tomographic  image  of  the 
target  is  reconstructed  based  on  multiple  transmitted  noise  waveforms  as  shown  in  Figure  2- 16(d). 

Based  on  the  numerical  simulation  results  of  various  scattering  target  geometries  using 
multiple  band-limited  iid  UWB  WGN  waveforms,  increasing  the  number  of  transmissions  of  the 
iid  UWB  WGN  waveform  improves  the  quality  of  tomographic  image  by  reducing  the  variance 
of  the  spectral  response  of  WGN  as  stated  in  section  2.2. 
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Figure  2-16.  The  final  tomographie  image  of  three  randomly  distributed  PEC  eylinders  in 
different  sizes  after  summing  and  averaging  proeess  with  the  (a)  one  transmitted  WGN  waveform 
image,  (b)  three  transmitted  WGN  waveform  images,  (e)  seven  transmitted  WGN  waveform 
images,  and  (d)  all  ten  transmitted  WGN  waveform  images. 


2.5  Image  Quality  Measure 

As  eoneluded  in  seetion  2.4,  inereasing  the  number  of  transmissions  of  the  iid  WGN 
waveform  tends  to  enhanee  the  formed  tomographie  image.  The  image  quality  measures  (IQMs) 
are  diseussed  to  determine  the  quality  of  tomographie  images  based  on  quantitative  analysis.  A 
good  subjeetive  assessment  is  required  to  evaluate  the  image  quality  and  the  performanee  of 
imaging  systems.  Various  methods  to  measure  the  image  quality  and  investigate  their  statistieal 
performanee  have  been  studied  [34,  35].  The  most  frequently  used  measures  are  deviations 
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between  the  reference  and  reconstructed  image  with  varieties  of  the  mean  square  error  (MSE)  [36, 
37].  The  reasons  for  the  widespread  popularity  of  the  analysis  based  on  MSE  calculations  are 
their  mathematical  tractability  and  the  fact  that  it  is  often  straightforward  to  design  systems  that 
minimize  the  MSE. 


2.5.1  Reference  image  generation 

Prior  to  measuring  the  image  quality  of  reconstructed  tomographic  images  shown  in  the 
previous  section,  a  “reference”  image  is  obtained  via  diffraction  tomography  proposed  in  Figure 
2-6  using  the  first  derivative  Gaussian  waveform.  The  first  derivative  Gaussian  waveform 
possesses  the  desirable  property  in  that  it  has  zero  average  value.  Since  it  has  no  DC  component, 
the  first  derivative  Gaussian  waveform  is  the  most  popular  pulse  shape  for  various  UWB 
applications  [38,  39];  thus,  it  can  be  considered  to  be  an  appropriate  waveform  for  generating 
reference  images  to  compare  with  our  reconstructed  images.  Figure  2- 17(a)  displays  the  first 
derivative  Gaussian  input  waveform  with  a  pulse  width  of  0.1  ns,  and  the  frequency  spectrum  of 
the  pulse  is  relatively  flat  in  the  range  of  8  GHz  -10  GHz  as  shown  in  Figure  2- 17(b). 

The  Fourier  space  data  and  the  tomographic  image  of  four  scattering  target  geometries  in 
section  2.4  are  generated  with  the  first  derivative  Gaussian  input  waveform  shown  in  Figure 
2- 17(a),  and  they  are  displayed  in  Figure  2-18,  Figure  2-19,  Figure  2-20,  and  Figure  2-21.  These 
tomographic  images  are  considered  to  be  the  reference  images  for  IQM  analysis. 

The  Fourier  space  data  and  the  tomographic  image  of  four  scattering  target  geometries  in 
section  2.4  are  generated  with  the  first  derivative  Gaussian  input  waveform  shown  in  Figure 
2- 17(a),  and  they  are  displayed  in  Figure  2-18,  Figure  2-19,  Figure  2-20,  and  Figure  2-21.  These 
tomographic  images  are  considered  to  be  the  reference  images  for  IQM  analysis. 
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Figure  2-17.  (a)  The  first  derivative  Gaussian  input  waveform  with  a  pulse  width  of  0.1  ns.  (b) 
The  frequeney  speetrum  of  the  first  derivative  Gaussian  input  waveform  shown  in  Figure  2- 17(a). 
The  frequeney  speetrum  ranges  are  displayed  from  DC  to  12  GHz  only. 


[rad/m]  x  [m] 


Figure  2-18.  (a)  The  normalized  magnitude  of  Fourier  spaee  data  based  on  the  transmission  of 
the  first  derivative  Gaussian  waveform  shown  in  Figure  2- 17(a)  for  the  single  eylindrieal  PEC 
objeet  shown  in  Figure  2-3.  (b)  The  tomographie  image  of  the  single  eylindrieal  eondueting 
objeet  loeated  at  (0  em,  0  em).  The  eolorbar  indieates  the  normalized  magnitude  of  seattering 
objeet  funetion,  o(x,  y). 
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Figure  2-19.  (a)  The  magnitude  of  Fourier  spaee  data  based  on  the  transmission  of  the  first 
derivative  Gaussian  waveform  shown  in  Figure  2- 17(a)  for  two  symmetrieally  distributed 
eylindrieal  PEC  objeets  shown  in  Figure  2-11.  (b)  The  tomographie  image  of  two  symmetrieally 
distributed  eylindrieal  eondueting  objeets  loeated  at  (0  em,  -40  em)  and  (0  em,  40  em). 
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Figure  2-20.  (a)  The  magnitude  of  Fourier  spaee  data  based  on  the  transmission  of  the  first 
derivative  Gaussian  waveform  shown  in  Figure  2- 17(a)  for  two  randomly  distributed  eylindrieal 
PEC  objeets  shown  in  Figure  2-13.  (b)  The  tomographie  image  of  two  randomly  distributed 
eylindrieal  eondueting  objeets  loeated  at  (-7.5  em,  22.5  em)  and  (35  em,  -55  em). 
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Figure  2-21.  (a)  The  magnitude  of  Fourier  spaee  data  based  on  the  transmission  of  the  first 
derivative  Gaussian  waveform  shown  in  Figure  2- 17(a)  for  three  randomly  distributed  eylindrieal 
PEC  objeets  in  different  sizes  shown  in  Figure  2-15.  (b)  The  tomographie  image  of  three 
randomly  distributed  eylindrieal  eondueting  objeets  in  different  sizes  loeated  at  (-22.5  em, 
-20  em),  (10  em,  50  em),  and  (25  em,  -40  em). 


2.5.2  Pixel  difference-based  measure 

The  tomographie  images  shown  in  seetion  2.4  are  obtained  based  on  10  iid  UWB  WGN 
transmitted  waveforms,  and  the  image  enhaneement  teehnique  is  not  implemented  in  diffraetion 
tomography  algorithm.  Every  pixel  in  all  images  is  generated  based  on  the  seattering  objeet 
funetion  defined  in  (2.22).  MSE  is  the  eumulative  mean  squared  error  between  the  eorresponding 
pixels  of  the  referenee  images  and  the  reeonstrueted  images  based  on  multiple  iid  UWB  WGN 
transmitted  waveforms.  MSE  is  defined  as 

1  M-\N-\  2 

MSE  =  (2.24) 

MlM  y  =  Q  ^  =  0 

where  S(x^,y^)  represent  the  value  of  eaeh  pixel  of  the  referenee  image  and 

reeonstrueted  tomographie  images  with  multiple  iid  UWB  WGN  transmitted  waveforms, 
respeetively,  and  M  and  N  are  the  dimensions  of  the  images. 
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Each  pixel  of  tomographic  images  represents  the  normalized  magnitude  of  scattering 
object  function,  o(x, and  the  dimensions  of  images  are  equivalent  to  the  number  of  pixels 
across  the  x-axis  and  j;-axis.  By  using  (2.24),  MSE  values  are  calculated  to  evaluate  the 
deviation  of  the  pixel  values  from  those  of  the  corresponding  reference  image,  as  the  iid  noise 
images  are  summed,  averaged,  and  normalized. 


MSE  vs.  Number  of  Images 


Figure  2-22.  MSE  versus  number  of  reconstructed  images  after  normalization.  MSE  decreases  as 
the  number  of  reconstructed  images  is  overlapped. 


Figure  2-22  displays  the  MSE  versus  number  of  overlapping  reconstructed  images.  Such 
quantitative  analysis  shows  that  the  tomographic  images  of  the  objects  become  clear  and 
distinctive  as  the  number  of  reconstructed  images  based  on  iid  WGN  waveforms  increases.  To  be 
more  specific,  as  shown  in  Figure  2-22,  MSE  decreases  as  the  number  of  overlapped 
reconstructed  images  based  on  iid  UWB  WGN  transmitted  waveforms  increases. 
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2.5.3  Signal-to-noise  ratio  effects  on  image  quality 

All  tomographic  images  shown  above  are  formed  in  the  absence  of  noise  in  the  system; 
that  is,  signal-to-noise  ratio  (SNR)  is  equal  to  infinity.  In  practical  situations,  the  received  signal 
shows  unpredictable  perturbations  due  to  the  contributions  of  various  noise  sources  in 
communication  and  imaging  systems  [40,  41].  For  example,  thermal  noise  is  present  in  all 
electronic  devices  and  transmission  media  and  is  uniformly  distributed  across  the  frequency 
spectrum.  When  backscattering  data  acquisition  is  performed  with  a  signal  analyzer,  the  recorded 
data  consists  of  scattered  field  information  with  some  thermal  noise  from  the  equipment  and 
cables.  However,  the  collected  data  can  be  viewed  as  the  sum  of  the  scattered  field  and  the 
additive  white  Gaussian  noise  since  they  are  uncorrelated. 

In  this  section,  the  white  Gaussian  noise  is  added  accordingly  to  the  collected  scattered 
field  dataset,  establishing  simulation  environment  with  4  different  SNR  values:  -20  dB,  -10  dB, 

0  dB,  and  +10  dB.  Prior  to  measuring  the  quality  of  the  final  tomographic  images  with  various 
SNRs,  MSB  values  are  calculated  to  determine  the  deviation  of  the  pixel  values  from  that  of  the 
corresponding  reference  image,  as  the  iid  noise  images  generated  with  the  corresponding  SNR 
values  are  summed,  averaged,  and  normalized.  Based  on  the  previous  quantitative  analysis 
results  shown  in  Figure  2-22,  the  image  quality  of  the  tomographic  image  is  expected  to  be 
enhanced  as  the  number  of  reconstructed  images  increases  regardless  of  SNR  values.  Figure  2-23 
displays  the  MSB  versus  number  of  overlapping  reconstructed  images  with  4  different  SNRs  for 
four  scattering  target  geometries  shown  in  section  2.4. 
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MSE  vs.  Number  of  Images  (Single  Cylindrical  PEC  Object)  MSE  vs.  Number  of  Images  (2  Symmetrically  Distributed  PEC  Objects) 


MSE  vs.  Number  of  Images  (2  Randomly  Distributed  PEC  Objects)  MSE  vs.  Number  of  Images  (3  Cylindrical  PEC  Objects) 


Figure  2-23.  MSE  versus  number  of  reeonstrueted  images  after  normalization  with  4  different 
SNR  values  for  (a)  single  eylindrieal  PEC  objeet,  (b)  two  symmetrieally  distributed  eylindrieal 
PEC  objeets,  (e)  two  randomly  distributed  eylindrieal  PEC  objeets,  and  (d)  three  randomly 
distributed  eylindrieal  PEC  objeets  in  different  sizes.  For  all  eases,  MSE  deereases  as  the  number 
of  reeonstrueted  images  is  overlapped  regardless  of  SNRs. 


After  eompletion  of  total  10  transmissions  of  iid  noise  waveforms  for  4  different  SNR 
values  defined  above,  the  final  tomographie  images  of  all  four  seattering  geometries  are  aehieved 
and  displayed  in  Figure  2-24  through  Figure  2-27. 
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Figure  2-24.  The  final  tomographic  images  after  summing  and  averaging  process  with  the  10  iid 
noise  waveforms  with  various  SNRs  for  a  single  PEC  cylinder;  (a)  SNR  =  -20  dB,  (b)  SNR  = 
-10  dB,  (c)  SNR  =  0  dB,  and  (d)  SNR  =  +10  dB. 
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Figure  2-25.  The  final  tomographie  images  after  summing  and  averaging  proeess  with  the  10  iid 
noise  waveforms  with  various  SNRs  for  two  symmetrieally  distributed  PEC  eylinders:  (a)  SNR  = 
-20  dB,  (b)  SNR  =  -10  dB,  (e)  SNR  =  0  dB,  and  (d)  SNR  =  +10  dB. 
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Figure  2-26.  The  final  tomographie  images  after  summing  and  averaging  proeess  with  the  10  iid 
noise  waveforms  with  various  SNRs  for  two  randomly  distributed  PEC  eylinders:  (a)  SNR  = 
-20  dB,  (b)  SNR  =  -10  dB,  (e)  SNR  =  0  dB,  and  (d)  SNR  =  +10  dB. 
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Figure  2-27.  The  final  tomographie  images  after  summing  and  averaging  proeess  with  the  10  iid 
noise  waveforms  with  various  SNRs  for  three  randomly  distributed  PEC  eylinders  in  different 
sizes:  (a)  SNR  =  -20  dB,  (b)  SNR  =  -10  dB,  (e)  SNR  =  0  dB,  and  (d)  SNR  =  +10  dB. 


As  shown  in  Figure  2-24  through  Figure  2-27,  the  image  quality  of  the  final  tomographie 
image  after  summing,  averaging,  and  normalizing  proeess  with  transmitting  10  iid  noise 
waveforms  is  truly  affeeted  by  the  additive  white  Gaussian  noise.  By  visual  inspeetion  of  the 
formed  images,  the  image  degradation  due  to  the  additive  Gaussian  noise  is  elearly  displayed 
when  SNR  is  set  to  -20  dB.  In  eontrast,  no  sueh  image  degradation  is  observed  at  relatively 
higher  SNRs,  that  is,  SNR  of +10  dB. 

In  order  to  evaluate  the  overall  impaet  of  noise  effeet  on  the  formed  tomographie  images, 
MSE  is  ealeulated  at  the  SNR  range  from  -20  dB  to  +10  dB.  Figure  2-28  explains  the 
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relationship  between  the  image  quality  degradation  after  10  images  and  SNR.  The  rate  of  ehange 
in  MSB  is  the  maximum  at  the  SNR  range  between  -20  dB  and  -10  dB,  and  this  is  the  SNR  range 
where  the  final  tomographie  image  is  severely  degraded  by  unwanted  noise.  However,  the  rate  of 
ehange  in  MSB  starts  falling  as  SNR  inereases,  whieh  indieates  that  the  tomographie  image  is  less 
affeeted  by  the  noise  eontributions  at  high  SNR.  The  results  shown  in  Figure  2-28  are  the  rough 
measure  to  estimate  the  image  quality  on  the  final  tomographie  images  in  praetieal 
implementation  of  the  imaging  system  and  designs  for  aetual  experiment. 


SNR(dB) 


Figure  2-28.  MSB  versus  SNR  values  for  all  final  tomographie  images  shown  in  Figure  2-24(d) 
through  Figure  2-27(d).  MSB  deereases  as  the  level  of  SNR  inereases. 


2.6  Summary 

This  ehapter  shows  that  tomographie  images  of  seattering  objeets  are  sueeessfully 
aehieved  with  multiple  transmissions  of  random  noise  waveforms.  From  the  simulation  results,  a 
single  transmitted  UWB  WGN  may  not  be  suffieient  to  generate  the  eorreet  tomographie  image 
due  to  the  praetieal  implementation  of  a  transmitted  random  noise  waveform.  However,  multiple 
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iid  UWB  WGN  transmitted  waveforms  ean  bypass  the  shorteoming  of  a  single  transmission  of 
UWB  WGN  waveform,  forming  a  eorreet  image  of  target  objeets  by  summing  and  averaging 
diserete  objeet  images  based  on  eaeh  iid  WGN  waveform. 

Also  the  image  quality  of  the  tomographie  image  after  eompletion  of  multiple 
transmissions  of  iid  noise  waveforms  is  analyzed.  MSB  is  used  to  measure  the  image  quality  of 
all  tomographie  images  in  this  ehapter.  Image  quality  of  the  tomographie  images  based  on  the 
random  noise  waveform  is  enhaneed  as  the  number  of  iid  noise  waveform  transmissions  inereases. 
The  presenee  of  white  Gaussian  noise  degrades  the  image  quality;  however,  the  suppression  of 
unwanted  noise  eontributions  by  eontrolling  SNR  ean  help  to  aehieve  sueeessful  tomographie 
images  in  praetieal  radar  imaging  systems. 
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Chapter  3 

Principle  and  Simulations  of  UWB  Noise  Radar  Tomographic  System 

3.1  Overview 

Noise  radar  has  been  developed  and  refined  for  covert  detection  and  imaging  applications 
over  the  past  few  decades  due  to  several  advantages,  such  as  low  probability  of  interception  (LPI) 
[1],  electronic  countermeasure  (ECM)  [2],  and  relatively  simple  hardware  configurations  [3]. 
Further,  UWB  noise  radar  has  achieved  high-resolution  detection  and  imaging  in  foliage 
penetration  (FOPEN)  imaging  [4-5],  and  through-wall  imaging  (TWI)  [6-7]  by  employing  an 
ultra-wide  bandwidth. 

The  main  advantage  of  transmitting  a  random  noise  waveform  is  to  covertly  detect  and 
image  a  target  without  alerting  the  presence  of  radar  system.  Such  EPI  characteristic  of  the  noise 
radar  is  achieved  by  using  a  random  noise  waveform,  which  is  constantly  varying  and  never 
repeats  itself  exactly  [8].  The  random  noise  waveform  herein  is  mathematically  defined  as 
discrete  time  wide-sense  stationary  (WSS)  and  ergodic  random  process,  and  a  sequence  of 
independent  and  identically  distributed  random  variable  drawn  from  a  normal  distribution,  A^(0, 
cr^).  Also,  the  power  spectral  density  (PSD)  of  WGN  waveform  is  ideally  a  nonzero  constant  for 
all  frequencies.  In  practical  situations,  however,  the  PSD  of  the  random  noise  waveform  may  not 
be  uniform  across  the  wide  frequency  ranges  because  a  finite  number  of  random  noise  amplitude 
samples  must  be  chosen  for  WGN  waveform  generation.  In  order  to  bypass  such  a  shortcoming, 
multiple  sequences  of  iid  noise  waveforms  are  transmitted  to  flatten  the  spectral  density  without 
changing  the  mean  and  variance  [9]. 
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Radar  tomography  is  emerging  as  a  powerful  imaging  technique  for  obtaining  high- 
resolution  images  of  targets  [10].  The  principle  of  tomography  is  to  characterize  the  material 
property  profiles  of  the  target  objects  in  the  frequency  domain,  reconstructing  specific  scatters 
within  the  interrogation  medium  by  solving  the  inverse  scattering  problem  [11].  Since  the 
frequency  behavior  of  the  electrical  parameters  of  the  target  object  due  to  the  illuminating 
electromagnetic  fields  depends  on  the  frequency  used  in  the  system,  UWB  radar  imaging  clearly 
provides  advantages  over  single  or  narrow  band  frequency  operation  in  terms  of  resolution  and 
accuracy.  Furthermore,  transmitting  random  noise  waveform  inherits  the  aforementioned 
advantages  of  noise  radar. 

This  chapter  focuses  on  the  principle  and  numerical  simulation  results  of  tomographic 
imaging  of  various  target  scenarios  using  UWB  noise  radar.  Multiple  iid  UWB  noise  waveforms 
over  3GHz  to  5GHz  frequency  range  are  transmitted  to  obtain  the  tomographic  images  of  rotating 
cylindrical  objects.  In  section  3.2,  theoretical  analysis  of  the  image  reconstruction  for  a  bistatic 
imaging  radar  system  with  a  two-dimensional  rotating  cylinder  using  Fourier  diffraction  theorem 
under  the  assumption  of  plane  wave  illumination  is  presented.  The  numerical  simulation  results 
of  diffraction  tomography  using  UWB  random  noise  waveforms  are  presented  in  section  3.3. 
Conclusions  are  presented  in  section  3.4. 


3.2  Formulation  of  the  Scattering  Properties 

The  top  view  of  two-dimensional  backward  scattering  geometry  for  a  cylindrical  object  is 
shown  in  Figure  3-1.  The  incident  z-polarized  plane  wave  is  illuminated  from  -x  direction.  The 
receiver  spacing.  Ay,  and  the  frequency  sampling  intervals.  A/  are  calculated  accordingly  to  avoid 
aliasing  using  (2.8)  and  (2.9)  from  section  2.3.  A  linear  receiving  array  is  located  at  x  =  -<7  from 
the  center  of  the  object,  collecting  scattered  field  for  reconstruction  of  the  object  image. 
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Figure  3-1.  Two-dimensional  backward  scattering  geometry  for  a  cylindrical  object.  Red  dots 
and  green  circle  represent  a  linear  receiving  array  and  the  target  object,  respectively. 

The  incident  z-polarized  plane  wave  of  the  white  Gaussian  noise  waveform  for  N 
discrete  frequencies  takes  the  form 


E,  {k„,r)  =  z  + . . .  +  ) 


N 


(3.1) 


n=l 


where  k^^co^t  c  is  the  wavenumber,  rj  is  the  intrinsic  impedance  in  free  space,  and  is  the 

field  amplitude  of  the  transmitted  white  Gaussian  noise  waveform  at  each  discrete  frequency  of 
interest.  We  start  the  analysis  from  a  single  frequency  and  develop  the  process  for  N  multiple 
frequencies  by  summing  up  the  analysis  results  for  N.  The  entire  analysis  must  be  repeated  for  K 
times  when  K  multiple  sequences  of  iid  noise  waveforms  are  transmitted. 
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In  this  chapter,  the  target  object  is  assumed  to  be  a  dielectric  cylinder  of  the  complex 
permittivity  ^(r ) ,  and  the  object  is  surrounded  by  free  space.  The  incident  plane  wave  denoted 
in  (3.1)  induces  an  equivalent  electric  current  distribution,  and  is  defined  as 

dielectric  j6)[s{r)-s^]E(k^,r).  (3.2) 

E{k^,r)  in  (3.2)  is  the  total  field  due  to  the  presence  of  the  dielectric  object,  and  expressed  as  the 
sum  of  the  incident  field  and  the  scattered  field,  that  is, 

E{k„,r)  =  EXK,r)  +  EXK,r).  (3.3) 

The  general  form  of  the  scattered  field  {k^,r)  created  by  the  equivalent  induced 

current  is  expressed  as 

Es (K ’ =  -J(OMo I (k„,n-G{r-  r')dr’  (3 .4) 

where  S  is  the  boundary  of  scatterer,  and  G(r -r')  is  the  two-dimensional  Green’s  function 
defined  as 

G(r-r')  =  ^//o<^'(yt„(r-r'))  (3.5) 

where  is  the  zeroth-order  Hankel  function  of  the  second  kind. 

For  the  dielectric  cylinder,  the  first-order  Born  approximation  is  assumed  to  be  valid  [12]. 
In  this  case,  the  z-polarized  scattered  field,  dielectric  ’  receivers  at  x  =  -J,  is  given  by 

=  -d,y)  =  Z  •  ^.,dielectrio(^«»^  = 

where  is  the  complex  relative  dielectric  constant  =  s  I  s^) .  The  general  form  of  the  z- 
polarized  scattered  field,  ,  is  written  for  the  simplicity,  that  is, 

uXK,x^ -d,y)  =  (^11  G{r  -r')d^r'  (3.7) 
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where  y/  and  o{f)  represent  the  scaling  factor  outside  the  integral  and  the  scattering  object 
function  of  the  target  object,  respectively.  For  a  PEC  object,  and  o(r')  are  jk^E^  and 
[2/2(r')  ,  respectively,  as  previously  expressed  in  (2.15)  and  (2.16)  from  section  2.3. 

For  a  dielectric  object,  y/  and  o(r')  are  and  [^^(r')-l] ,  respectively,  as  given  in  (3.6).  The 

scattering  object  function  is  considered  to  be  the  electrical  contrast  between  the  target  object  and 
free  space. 

By  using  the  plane  wave  expansion  of  the  Green’s  function  [13],  the  one-dimensional 
Fourier  transform  of  the  scattered  field,  ,  is  given  by 


Us(K^^  =  -d.ky)  = 


-J¥ 
ly  ' 


-jyd 


o{-r-K,k) 


(3.8) 


where 


as  y\<k„ 
as  \k)^>k„ 


(3.9) 


If  we  define  k^=-y-k^,  the  scattering  object  function  in  the  Fourier-domain,  0{k^,ky),  is  a 
two-dimensional  Fourier  transform  of  the  scattering  function: 

d{k^,ky)-  jj dxdy .  (3.10) 

From  (3.9),  the  arguments  of  0{k^,ky)  are  related  by 

Equations  (3.8)  and  (3.11)  show  that  as  a  two-dimensional  scattering  object  is  illuminated  by  a 
plane  wave  at  the  frequency  cd^  ,  the  one-dimensional  Fourier  transform  of  the  scattered  field 
yields  a  semicircle  centered  at  {0,-kJ  with  radius  k^  in  the  two-dimensional  Fourier  space, 
which  relates  to  the  Fourier  diffraction  theorem  described  in  [14]. 
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The  scattering  object  function,  o(x,  y),  of  the  target  is  imaged  in  rectangular  coordinates 
after  performing  the  two-dimensional  inverse  Fourier  transform  of  (3.10).  With  the  multiple  N 

frequency  data  stepped  from  to  kj^ ,  the  radius  of  0{k^,ky)  in  the  Fourier-domain  changes 

accordingly,  and  N  sets  of  the  scattering  object  function  enhance  the  resolution  and  accuracy  of 
the  image  [15-16]. 


Figure  3-2.  The  image  reconstruction  method  with  K  multiple  iid  WGN  transmitted  waveforms 
for  a  sequentially  rotated  object  using  diffraction  tomography. 


As  the  object  is  sequentially  rotated  with  respect  to  the  rotational  axis  for  the  value  of 
rotational  angle  9,  the  corresponding  scattering  object  function  at  the  angle  9  is  obtained.  The 
complete  two-dimensional  tomographic  image  of  the  target  is  reconstructed  via  computing  the 
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sum  of  all  formed  images  at  each  rotational  angle.  A  block  diagram  shown  in  Figure  3-2  displays 
the  tomographic  image  reconstruction  method  for  multiple  transmitted  iid  noise  waveforms  for  a 
sequentially  rotated  object  using  diffraction  tomography. 


3.3  Numerical  Simulation  Results 


3.3.1  Image  quality  measure 

As  discussed  in  chapter  2,  the  image  quality  is  determined  by  various  factors  such  as  the 
number  of  transmissions  of  iid  noise  waveforms  and  the  number  of  random  amplitude  samples  of 
the  transmitted  noise  waveform.  In  this  chapter,  the  image  quality  of  the  complete  tomographic 
image  must  be  quantitatively  analyzed  for  different  values  of  the  rotational  angle  since  the 
complete  tomographic  image  is  achieved  based  on  the  rotation  of  the  target  object. 


First  Derivative  Gaussian  Waveform 


X  =  -68 


Figure  3-3.  Two-dimensional  backward  scattering  simulation  geometry  for  a  cylindrical 
conducting  object.  The  PEC  cylinder  with  a  radius  of  8  cm  is  located  at  (0  cm,  0  cm),  and  a 
linear  receiving  array  is  located  68  cm  away  from  the  center  of  the  object. 
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Figure  3-3  shows  the  two-dimensional  baekward  seattering  geometry  with  a  single 
eylindrieal  eondueting  objeet.  The  PEC  ey Under  with  a  radius  of  8  em  is  loeated  at  (0  em,  0  em), 
and  the  eylinder  is  assumed  to  be  infinitely  long  along  the  z-axis.  The  seattered  field  is  uniformly 
sampled  at  reeeiving  array  Rxl  through  Rx55  with  frequeney  swept  within  the  frequeney  ranges 
from  3-5  GHz.  The  loeations  of  Rxl  and  Rx55  are  at  (-68  em,  -40.5  em)  and  (-68  em,  40.5  em), 
respeetively,  and  the  reeeiver  spaeing  is  set  to  1.5  em. 

After  seattered  field  data  are  eolleeted  at  all  reeeiving  array  Rxl  through  Rx55,  the  objeet 
is  rotated  with  respeet  to  the  rotational  axis  (0  em,  0  em)  for  the  value  of  rotational  angle  9.  The 
simulation  repeats  until  the  objeet  is  fully  rotated  for  360°,  and  the  eomplete  seattering  data  for 
two-dimensional  tomographie  image  of  the  target  are  obtained.  The  minimum  rotational  angle  9 
is  1.8°,  whieh  is  the  eontrollable  rotational  angle  of  the  turntable  for  the  experiment.  The  value  of 
9  varies  from  3.6°  to  54°  to  determine  the  optimum  rotational  angle  by  ealeulating  mean  square 
error  (MSE). 

MSE  is  the  eumulative  mean  squared  error  between  the  eorresponding  pixels  of  the 
referenee  images  and  the  reeonstrueted  images.  MSE  is  defined  as 

1  M-\N-\ 

MSE  =  mXRix^^yn)-S{x„,y„)]^  (3.12) 

^=0  n=0 

where  R{x^,y^)  and  S(x^,y^)  represent  the  value  of  eaeh  pixel  of  the  referenee  image  and 

reeonstrueted  tomographie  images  based  on  the  transmitted  waveform,  respeetively,  and  M  and 
N  are  the  pixel  dimensions  of  images.  The  first  derivative  Gaussian  waveform  with  a  pulse  width 
of  0.225  ns  is  transmitted  for  all  MSE  ealeulations.  As  shown  in  Figure  3-4,  the  frequeney 
response  of  the  transmitted  Gaussian  waveform  in  the  ranges  of  3-5  GHz  is  relatively  fiat. 
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Figure  3-4.  (a)  The  first  derivative  Gaussian  input  waveform  with  a  pulse  width  of  0.225  ns.  (b) 
The  frequeney  speetrum  of  the  first  derivative  Gaussian  input  waveform  shown  in  Figure  3 -4(a). 
The  frequeney  speetrum  range  is  displayed  from  DC  to  8  GHz  only 


Figure  3-5  shows  the  tomographie  image  of  the  eylinder  with  a  rotational  angle  of  1.8°, 
whieh  ean  be  eonsidered  as  a  referenee  image  to  measure  the  deviations  of  image  quality  when 
the  value  of  the  rotational  angle  ehanges. 


X  (m) 

Figure  3-5.  The  reference  tomographic  image  of  the  circular  PEC  cylinder  shown  in  Figure  3-3 
iove=  1.8°. 
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X  (m)  X  (m) 


Figure  3-6.  Tomographic  images  of  the  circular  PEC  cylinder  with  6  different  rotational  angles: 
(a)  9  =  3.6°,  (b)  9  =  9°,  (c)  6/  =  18°,  (d)  9  =  36°,  (e)  9  =  45°,  and  {^)9  =  54°. 


Figure  3-6  shows  the  set  of  6  tomographic  images  when  the  rotational  angle  6  varies 
from  3.6°  to  54°.  As  shown  in  Figure  3-6,  the  successful  tomographic  images  of  the  circular 
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cylindrical  object  for  various  rotational  angles  are  aehieved  until  the  value  of  the  rotational  angle 
9  is  set  to  36°.  For  large  rotational  angles,  i.e.,  9  =  45°  and  9  =  54°,  the  tomographie  images  of 
the  objeet  beeome  blurry  and  slightly  inaeeurate  by  visual  inspeetion  beeause  the  are  of  the 
eylinder  eannot  be  preeisely  aehieved  due  to  the  limited  number  of  the  linear  reeeiving  array  and 
pixels  of  the  formed  images  [17]. 

MSB  vs.  the  rotational  angle  9  is  plotted  in  Figure  3-7.  Sueh  quantitative  analysis 
eoneludes  that  the  image  quality  deereases  as  the  value  of  rotational  angle  9  inereases  [18].  All 
images  shown  in  Figure  3-5  and  Figure  3-6  are  reeonstrueted  based  on  numerieal  simulations 
without  any  eonsideration  of  unwanted  additive  noise,  that  is,  the  signal-to-noise  ratio  (SNR)  is 
equal  to  infinity.  Thus,  MSB  purely  depends  on  the  value  of  rotational  angle  9. 


X  -10’^  \«.  Rotational  Angle 
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Rotational  Angle  (Degree) 

Figure  3-7.  MSB  vs.  rotational  angle  for  the  eireular  PBC  eylinder  simulation.  MSB  inereases  as 
the  value  of  the  rotational  angle  0  inereases. 

3.3.2  Multiple  transmitted  iid  WGN  waveforms  and  rotational  angle 

As  diseussed  in  the  previous  seetion,  the  tomographie  image  of  the  seattered  objeet  at 
eaeh  rotational  angle  9  is  reeonstrueted  via  sum  of  all  K  diserete  images  for  AT-th  iid  transmitted 
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UWB  noise  waveforms.  The  discrete  image  at  the  rotational  angle  6  is  obtained  based  on  the 
corresponding  scattering  object  function  at  the  angle  9,  and  the  complete  image  of  the  target 
cylinder  is  achieved  by  plotting  all  discrete  images  in  the  single  image  plane  as  shown  in  Figure 
3-2.  For  scattered  field  and  diffraction  tomography  simulations  with  multiple  iid  WGN 
waveforms,  10  iid  UWB  WGN  waveforms  over  a  frequency  range  from  3  GHz  to  5  GHz  are 
generated  with  500  random  amplitude  samples  drawn  from  A^(0,  cr^),  that  is,  10,  and  they  are 
sequentially  transmitted  for  backward  scattered  field  calculations.  The  rotational  angle  9  is  set  to 
9°  based  on  the  MSB  analysis  displayed  in  Figure  3-7  such  that  40  discrete  images  are  obtained 
after  full  360°  rotation.  The  complete  tomographic  image  of  the  target  is  achieved  by  combining 
40  discrete  images  into  one  single  image.  Table  3-1  shows  the  list  of  target  objects  for  the 
numerical  simulations. 


Table  3-1.  Target  objects  imaged  by  the  numerical  simulations. 


Object  Material 

Object  Description 

Dimension 

PEC 

Circular  Cylinder 

20  cm  diameter 
(7.9"  diameter) 

Height  of  the 
objects  are 
assumed  to  be 
infinitely  long 

Square  Box 

20  cm  X  20  cm 
(7.9"  X  7.9") 

Triangular  Prism 

20  cm  side 
(7.9"  side) 

Dielectric 

=5.8-70.75 

Circular  Cylinder 

1 6  cm  diameter 
(6.3"  diameter) 

Square  Box 

16  cm  X  16  cm 
(6.3"  X  6.3") 
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3.3.3  Tomographic  image  with  multiple  transmitted  iid  WGN  waveforms  for  a  circular 
PEC  cylinder 

As  shown  in  Figure  3-8,  two-dimensional  backward  scattering  geometry  with  a  circular 
PEC  cylindrical  object  is  simulated  with  10  iid  UWB  WGN  waveforms.  The  cylinder  with  a 
radius  of  10  cm  is  located  at  (0  cm,  0  cm),  and  the  height  of  the  cylinder  is  assumed  to  be 
infinitely  long  along  the  z-axis.  The  locations  of  Rxl  and  Rx55  are  at  (-68  cm,  -40.5  cm)  and 
(-68  cm,  40.5  cm),  respectively,  and  the  receiver  spacing  is  set  to  1.5  cm  in  order  to  avoid  any 
aliasing. 


Incident  noise  waveform 


Figure  3-8.  Two-dimensional  backward  scattering  simulation  geometry  for  a  circular  PEC 
cylinder.  The  cylinder  with  a  radius  of  10  cm  is  located  at  (0  cm,  0  cm),  and  a  linear  receiving 
array  is  located  68  cm  away  from  the  center  of  the  object. 


The  scattered  field  is  uniformly  sampled  at  receiving  array  Rxl  through  Rx55  with 
frequency  swept  within  the  frequency  ranges  from  3-5  GHz.  When  the  data  acquisition  is 
completed  for  all  receiving  elements  Rxl  through  Rx55,  the  object  is  rotated  at  9°  with  respect  to 
the  origin  (0  cm,  0  cm).  The  simulation  repeats  until  the  object  is  fully  rotated  for  360°. 
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Four  tomographic  images  of  the  circular  PEC  cylinder  shown  in  Figure  3-8  with  10 
transmitted  iid  WGN  waveform  are  shown  in  Figure  3-9(a),  Figure  3-9(b),  Figure  3-9(c)  and 
Figure  3 -9(d)  when  1,  3,  7  and  10  noise  waveforms  are  transmitted,  respectively. 


Figure  3-9.  Two-dimensional  tomographic  image  of  a  single  PEC  cylinder  with  a  radius  of  10  cm 
located  at  (0  cm,  0  cm)  after  summing  and  averaging  process  based  on  (a)  one  transmitted  WGN 
waveform,  (b)  three  transmitted  WGN  waveforms,  (c)  seven  transmitted  WGN  waveforms,  and 
(d)  all  ten  transmitted  WGN  waveforms.  The  rotational  angle  is  9°. 
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3.3.4  Tomographic  image  with  multiple  transmitted  iid  WGN  waveforms  for  a  square  PEC 
box 


Figure  3-10  shows  the  two-dimensional  backward  scattering  geometry  with  a  square  PEC 
box  located  at  (0  cm,  0  cm).  The  length  of  each  side  of  the  square  (a)  is  20  cm,  and  the  height  of 
the  box  is  assumed  to  be  infinitely  long  along  the  z-axis.  The  locations  of  all  55  receivers,  the 
parameters  for  receiver  spacing.  Ay  ,  and  the  rotational  angle,  9  ,  are  identical  to  the  previous  case. 

In  order  to  generate  the  complete  tomographic  image  of  the  square  PEC  box  shown  in 
Figure  3-10,  10  iid  random  noise  waveforms  are  sequentially  transmitted.  Figure  3-1 1(a),  Figure 
3-1 1(b),  Figure  3-1 1(c)  and  Figure  3-1 1(d)  display  four  tomographic  images  of  the  square  PEC 
box  when  1,  3,  7  and  10  noise  waveforms  are  transmitted,  respectively. 


Incident  noise  waveform 


x  =  -68 


Figure  3-10.  Two-dimensional  backward  scattering  simulation  geometry  for  a  square  PEC  box 
with  a  length  of  20  cm.  The  box  is  located  at  (0  cm,  0  cm),  and  a  linear  receiving  array  is  located 
68  cm  away  from  the  center  of  the  object. 
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Figure  3-11.  Two-dimensional  tomographie  image  of  a  square  PEC  box  with  a  length  of  20  em 
loeated  at  (0  em,  0  em)  after  summing  and  averaging  proeess  based  on  (a)  one  transmitted  WGN 
waveform,  (b)  three  transmitted  WGN  waveforms,  (e)  seven  transmitted  WGN  waveforms,  and 
(d)  all  ten  transmitted  WGN  waveforms.  The  rotational  angle  is  9°. 
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3.3.5  Tomographic  image  with  multiple  transmitted  iid  WGN  waveforms  for  an  equilateral 
triangular  PEC  prism 

Backward  scattering  simulation  geometry  with  an  equilateral  triangular  PEC  prism  is 
shown  in  Figure  3-12.  The  triangular  prism  is  located  at  (0  cm,  0  cm),  and  the  length  of  each  side 
of  the  triangle  {a)  is  20  cm.  The  target  object  is  also  considered  as  an  infinitely  long  prism  along 
the  z-axis.  Other  simulation  parameters,  such  as  the  position  of  receivers,  the  parameters  for 
receiver  spacing,  and  the  rotational  angle,  remain  unchanged.  Also,  four  tomographic  images  of 
the  equilateral  triangular  PEC  prism  with  1,3,7  and  10  transmitted  iid  noise  waveforms  are 
shown  in  Figure  3-13(a),  Figure  3-13(b),  Figure  3-13(c),  and  Figure  3-13(d),  respectively. 


Incident  noise  waveform 


x  =  -68 


Figure  3-12.  Two-dimensional  backward  scattering  simulation  geometry  for  an  equilateral 
triangular  PEC  prism  with  a  length  of  20  cm.  The  prism  is  located  at  (0  cm,  0  cm),  and  a  linear 
receiving  array  is  located  68  cm  away  from  the  center  of  the  object. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


77 


63 


Figure  3-13.  Two-dimensional  tomographie  image  of  an  equilateral  triangular  PEC  prism  with  a 
length  of  20  em  loeated  at  (0  em,  0  em)  after  summing  and  averaging  proeess  based  on  (a)  one 
transmitted  WGN  waveform,  (b)  three  transmitted  WGN  waveforms,  (e)  seven  transmitted  WGN 
waveforms,  and  (d)  all  ten  transmitted  WGN  waveforms.  The  rotational  angle  is  9°. 
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3.3.6  Tomographic  image  with  multiple  transmitted  iid  WGN  waveforms  for  a  circular 
dielectric  cylinder 

Figure  3-14  shows  two-dimensional  backward  scattering  geometry  with  a  circular 
dielectric  cylinder  is  simulated  with  10  iid  UWB  WGN  waveforms.  The  circular  cylinder  with  a 
radius  of  8  cm  is  located  at  (0  cm,  0  cm),  and  the  height  of  the  cylinder  is  assumed  to  be  infinitely 
long  along  the  z-axis.  The  locations  of  Rxl  and  Rx55  are  at  (-68  cm,  -40.5  cm)  and  (-68  cm, 
40.5  cm),  respectively,  and  the  receiver  spacing  is  set  to  1.5  cm  in  order  to  avoid  any  aliasing 
effect. 

The  dielectric  cylinder  is  assumed  to  be  a  homogeneous  concrete  block;  the  real  and 
imaginary  values  of  complex  relative  permittivity  are  5.8  and  0.75,  respectively,  at  the  evaluation 
frequency  of  3  GHz  [19].  The  moisture  contents  and  surface  roughness  of  the  dielectric  cylinder 
are  not  considered  for  the  simulations. 


z 


Incident  noise  waveform 
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Figure  3-14.  Two-dimensional  backward  scattering  simulation  geometry  for  a  circular  dielectric 
cylinder  with  a  radius  of  8  cm.  The  cylinder  is  located  at  (0  cm,  0  cm),  and  a  linear  receiving 
array  is  located  68  cm  away  from  the  center  of  the  object. 
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The  scattered  field  is  uniformly  sampled  at  receiving  elements  Rxl  through  Rx55  with 
frequency  swept  within  the  frequency  ranges  from  3-5  GHz.  When  the  data  acquisition  is 
completed  for  all  receiving  array  Rxl  through  Rx55,  the  object  is  rotated  at  9°  with  respect  to  the 
rotational  axis  (0  cm,  0  cm).  The  simulation  repeats  until  the  object  is  fully  rotated  for  360°. 

Four  tomographic  image  of  the  dielectric  circular  cylinder  based  on  10  transmitted  iid 
WGN  waveform  are  shown  in  Figure  3-15(a),  Figure  3-15(b),  Figure  3-15(c)  and  Figure  3-15(d) 
when  1,  3,  7  and  10  noise  waveforms  are  transmitted,  respectively. 


Figure  3-15.  Two-dimensional  tomographic  image  of  a  single  circular  dielectric  cylinder  with  a 
radius  of  8  cm  located  at  (0  cm,  0  cm)  after  summing  and  averaging  process  based  on  (a)  one 
transmitted  WGN  waveform,  (b)  three  transmitted  WGN  waveforms,  (c)  seven  transmitted  WGN 
waveforms,  and  (d)  all  ten  transmitted  WGN  waveforms.  The  rotational  angle  is  9°. 
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3.3.7  Tomographic  image  with  multiple  transmitted  iid  WGN  waveforms  for  a  square 
dielectric  box 

Figure  3-16  shows  the  two-dimensional  backward  scattering  geometry  with  a  square 
dielectric  box  located  at  (0  cm,  0  cm).  The  length  of  each  side  of  the  square  is  16  cm,  and  the 
height  of  the  box  is  assumed  to  be  infinitely  long  along  the  z-axis.  The  material  properties  of  the 
dielectric  square  box,  i.e.,  the  complex  dielectric  constant,  are  identical  to  the  previous  dielectric 
circular  cylinder  case. 

In  order  to  generate  the  complete  tomographic  image  of  the  dielectric  square  box,  10  iid 
random  noise  waveforms  are  transmitted.  Figure  3-17(a),  Figure  3-17(b),  Figure  3-17(c)  and 
Figure  3- 17(d)  display  four  tomographic  images  of  the  square  dielectric  box  when  1,3,7  and  10 
noise  waveforms  are  transmitted,  respectively. 


Incident  noise  waveform 


x  =  -68 


Figure  3-16.  Two-dimensional  backward  scattering  simulation  geometry  for  a  square  dielectric 
box  with  a  length  of  20  cm  is  located  at  (0  cm,  0  cm),  and  a  linear  receiving  array  is  located  68 
cm  away  from  the  center  of  the  object. 
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Figure  3-17.  Two-dimensional  tomographie  image  of  a  square  dieleetrie  box  with  a  length  of  20 
em  loeated  at  (0  em,  0  em)  after  summing  and  averaging  proeess  based  on  (a)  one  transmitted 
WGN  waveform,  (b)  three  transmitted  WGN  waveforms,  (e)  seven  transmitted  WGN  waveforms, 
and  (d)  all  ten  transmitted  WGN  waveforms.  The  rotational  angle  is  9°. 


3.4  Summary 

In  this  ehapter,  the  eomplete  tomographie  images  of  the  sequentially  rotated  objeets  are 
sueeessfiilly  aehieved  with  multiple  transmissions  of  iid  UWB  random  noise  waveforms.  Based 
on  the  simulation  results,  the  image  of  both  PEC  and  dieleetrie  objeets  are  eorreetly  reeonstrueted. 
The  image  quality  of  the  reeonstrueted  images  for  the  rotated  objeet  eases  depends  on  the 
rotational  angle.  MSB  analysis  quantitatively  explains  the  relationship  between  the  image  quality 
and  different  rotational  angles.  Image  quality  is  notieeably  degraded  as  the  value  of  the  rotational 
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angle  increases.  Also  the  multiple  transmissions  of  iid  UWB  WGN  waveforms  truly  enhance  the 
image  quality  for  the  rotated  object  cases  herein.  The  numerical  simulation  results  shown  in  this 
chapter  will  be  compared  to  the  experimental  results,  which  will  be  presented  in  the  following 
chapter. 
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Chapter  4 

Design  and  Implementation  of  UWB  Noise  Radar  Tomographic  System 

4.1  Overview 

Microwave  tomography  has  been  developed  and  refined  over  the  past  few  decades  to 
obtain  high-resolution  images  of  metallic  targets  and  dielectric  contrasts  embedded  within  a 
dielectric  medium.  The  objective  of  active  microwave  tomography  is  to  reconstruct  the  dielectric 
properties  of  a  body  illuminated  with  microwaves  from  a  measurement  of  the  scattered  fields. 
Conventional  microwave  tomography  systems  are  based  on  illuminating  the  body  with  a  plane 
wave  and  measuring  the  scattering  with  a  linear  array  of  probes.  A  practical  requirement  when 
using  this  method  is  the  need  for  mechanical  rotation  of  either  the  body  or  the  antenna  in  order  to 
obtain  measurements  in  different  views  [1].  Every  object,  when  inserted  into  an  electromagnetic 
field,  causes  a  well-defined  field  change,  and  diffraction  occurs  when  the  wavelength  of  the 
microwave  radiation  is  of  the  order  of  the  dimension  of  the  object.  A  one-to-one  relationship 
relating  the  scattered  field  to  the  object  complex  permittivity  can  be  obtained  within  the  Born 
approximation,  via  a  Fourier  transform  using  the  so-called  Fourier  Diffraction  Theorem  [2]. 

Several  hardware  systems  have  been  developed  and  implemented  over  the  years.  A 
microwave  tomographic  system  consisting  of  64  circularly  arranged  electronically  scanned 
antennas,  divided  into  32  transmitters  and  32  receivers,  operating  at  a  frequency  of  2.45  GHz  was 
developed  for  biological  tissue  imaging  [3].  A  C-system  operating  at  7.5  GHz  was  developed  for 
imaging  weakly  scattering  objects,  such  as  paper  cylinders  [4].  More  recently,  a  two-dimensional 
UWB  microwave  imaging  system  for  3-6  GHz  frequency  range  using  24  antenna  elements 
connected  to  a  vector  network  analyzer  via  a  2x24  port  matrix  switch  was  able  to  quantitatively 
reconstruct  dielectric  objects  [5].  Also  a  time-domain  UWB  tomographic  imaging  system 
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transmitting  a  differentiated  cosine-modulated  Gaussian  pulse  with  a  carrier  frequency  of  1  GHz 
was  used  to  reconstruct  plastic  phantoms  of  different  sizes  and  shapes  immersed  in  a  tank  filled 
with  pipe  water  [6]. 

In  this  chapter,  the  hardware  design  and  experimental  investigation  of  UWB  noise  radar 
tomographic  system  are  discussed.  First,  implementation  of  UWB  noise  radar  tomographic 
system  is  described  in  section  4.2.  The  data  acquisition  process  in  UWB  noise  radar  for 
diffraction  tomography  and  the  effective  tomographic  image  processing  technique  for  large  sets 
of  collected  data  samples  are  proposed  in  section  4.3.  In  section  4.4,  successful  experimental 
results  of  various  cross-sectional  shapes  of  metallic  and  dielectric  cylindrical  targets  are  presented. 
Conclusions  of  the  implementation  and  performance  of  UWB  noise  radar  tomographic  system  are 
discussed  in  section  4.5. 


4.2  Implementation  of  the  Hardware  System 

4.2.1  Experimental  configuration  of  UWB  noise  radar 

The  entire  experiment  is  sequentially  controlled  by  a  desktop  computer.  The  computer  is 
directly  connected  to  the  electrical  components:  an  arbitrary  waveform  generator  (AWG)  and  a 
digital  oscilloscope.  Before  the  data  acquisition  process  begins,  the  AWG  at  the  transmitter  side 
starts  generating  the  random  waveforms  based  on  a  given  specification  such  as  the  number  of 
amplitude  samples,  frequency  bandwidth,  and  sampling  rate.  The  AWG  is  designed  to  transmit 
the  random  noise  waveform  continuously  for  the  duration  of  the  measurement.  The  digital 
oscilloscope  is  configured  to  collect  and  save  the  scattering  data  at  the  receiver  side. 

Also  two  mechanical  components,  i.e.,  a  turntable  and  a  linear  scanner,  are  controlled  by 
the  desktop  computer.  For  data  acquisition,  the  target  object  is  placed  on  the  turntable,  and  a 
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stepper  motor  of  the  turntable  rotates  the  object  a  certain  number  of  degrees.  The  receiving 
antenna  is  mounted  on  the  linear  scanner;  the  antenna  moves  along  the  linear  scanner  and  stops  at 
the  designated  coordinate  to  collect  the  data  at  that  point.  The  receiving  antenna  is  connected 
directly  to  the  oscilloscope;  thus,  the  oscilloscope  acquires  time-domain  data  and  then  passes  the 
data  back  to  the  desktop  computer.  Figure  4-1  shows  the  proposed  hardware  implementation  of 
UWB  noise  radar  to  measure  the  scattering  data  of  the  target  object  in  a  backward  scattering 
arrangement. 


Figure  4-1.  Automated  UWB  noise  radar  tomographic  system. 
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4.2.2  Transmitter  and  receiver 

The  block  diagram  for  the  transmitter  and  the  receiver  portions  of  the  system  are  shown 
in  Figure  4-2.  The  AWG  is  Agilent  M8190A  system  mounted  on  a  M9602A  chassis.  For 
transmit  noise  waveforms,  the  AWG  generates  a  bandlimited  random  noise  waveform  from  of  3- 
5  GHz  at  an  output  power  of  0  dBm.  The  gain  of  ZVE-8G+  with  12  V  DC  input  power  across  the 
frequency  range  of  2-8GHz  is  plotted  in  Figure  4-3.  The  output  power  ranges  of  the  transmitted 
bandlimited  noise  waveform  after  the  amplifier  is  approximately  35±0.5  dBm  because  of  the 
tolerance  of  the  gain  of  ZVE-8G+. 


Rx  Antenna 


I 


Figure  4-2.  Block  diagram  of  transmitter  and  receiver. 


The  amplified  noise  waveform  is  transmitted  through  an  A-Info  dual-polarization  horn 
antenna.  The  dual-polarized  horn  antenna  operates  over  the  frequency  range  of  2-18  GHz,  and 
the  antenna  has  a  gain  of  approximately  10  dB  over  3-5  GHz  operating  range.  The  position  of 
the  transmitting  antenna  is  fixed  at  the  center  of  the  backward  scattering  geometry,  but  the 
receiving  antenna  is  designed  to  move  along  the  linear  scanner  to  collect  the  scattering  data.  The 
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receiving  side  of  the  system  uses  an  identical  antenna  for  collecting  scattering  data.  Both 
transmitting  and  receiving  antenna  are  installed  to  measure  only  vertically  polarized  scattering 
data.  The  receiving  antenna  is  connected  to  an  Agilent  Infiniium  DS090804A  oscilloscope 
which  has  a  maximum  sampling  rate  of  40  GSa/s.  The  oscilloscope  then  collects  the  scattering 
data  and  then  sends  it  to  the  computer  for  tomographic  image  processing. 
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Figure  4-3.  The  gain  vs.  frequency  plot  of  Mini-circuit  ZVE-8G+. 

4.2.3  Turntable  and  linear  scanner 

The  receiving  antenna  is  placed  on  a  linear  scanner  in  place  of  an  antenna  array,  as  shown 
in  Figure  4-2.  The  scanner  stops  at  55  locations  to  measure  the  vertically  polarized  scattering 
data.  Since  the  length  of  the  scanner  is  81  cm,  the  receiver  spacing.  Ay,  is  calculated  to  be  1.5  cm, 
satisfying  the  equation  (2.8).  Also  microwave  absorber  material  is  placed  behind  the  target  object 
in  order  to  attenuate  signal  reflected  off  of  the  back  wall,  as  well  as  to  minimize  the  multipath 
reflections. 


Gain  vs.  Frequency 
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Both  metallic  and  dielectric  cylinders  with  different  shapes  are  measured  for  tomographic 
image  reconstructions.  The  metallic  objects  are  constructed  from  sheets  of  aluminum.  For  the 
dielectric  target  objects,  they  are  assumed  to  be  homogeneous  concrete  blocks.  The  moisture 
contents  and  surface  roughness  of  the  concrete  blocks  are  not  considered  for  the  experiments. 

The  test  objects  are  placed  on  the  center  of  the  turntable  one  at  a  time.  Table  4-1  shows  the  target 
objects  for  which  the  tomographic  images  are  presented  herein. 


Table  4-1.  Target  objects  imaged  by  the  tomographic  noise  radar  system. 


Object  Material 

Object  Description 

Dimension 

Metal 

Circular  Cylinder 

20.3  cm  diameter  x  40.6  cm  height 
(8"  diameter  x  16"  height) 

Square  Box 

20.3  cm  X  20.3  cm  x  61  cm  height 

(8"  X  8"  X  24") 

Triangular  Prism 

20.3  cm  side  x  30.5  cm  height 
(8"  side  X  12"  height) 

Concrete 

Circular  Cylinder 

15.2  cm  diameter  x  30.5  cm  height 
(6"  diameter  x  12"  height) 

Square  Block 

15.2  cm  X  15.2  cm  x  28  cm  height 

(6"  X  6"  X  11") 

For  the  entire  measurement  process,  the  receiving  antenna  stops  in  order  to  ensure  that 
the  oscilloscope  has  sufficient  amount  of  time  to  collect  and  save  scattering  data  at  that  position. 
After  data  acquisition  for  all  55  points  of  data  are  complete,  the  linear  scanner  moves  back  to  its 
initial  position.  As  the  linear  scanner  is  moving  back  to  the  first  point,  the  turntable  rotates  by  9 
degrees.  The  linear  scanner  repeats  the  data  acquisition  for  same  55  scanning  points  as  before, 
and  such  measurement  process  are  continued  until  the  object  has  rotated  by  a  full  360  degrees. 
The  turntable  is  controlled  by  a  stepper  motor  and  is  able  to  rotate  in  any  multiple  of  1.8°.  Since 
the  rotational  angle  is  defined  as  9°  based  on  MSB  analysis  in  Section  3.3,  a  total  of  40  angular 
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measurement  dataset  are  obtained.  The  number  of  scanning  points  and  the  angle  of  rotation  were 
chosen  as  a  tradeoff  between  data  collection  time  as  well  as  image  quality. 


4.2.4  Data  acquisition  and  image  processing 

Figure  4-4  shows  the  actual  metallic  cylinder  and  the  complete  experimental  hardware 
configuration  with  the  cylinder.  Both  the  transmitting  and  the  receiving  antennas  point  towards 
the  foam  wall.  The  target  object  is  positioned  in  the  main  lobe  for  the  transmitting  and  the 
receiving  antenna.  The  data  collection  process  starts  with  the  Gaussian  noise  waveforms  being 
transmitted  to  the  target  object.  The  transmitted  noise  waveforms  are  not  identical  between  each 
measurement;  however,  the  iid  white  Gaussian  noise  waveforms  are  transmitted  for  the 
experiments.  Each  time  the  linear  scanner  stops,  the  oscilloscope  collects  a  total  of  363,000 
amplitude  samples.  As  previously  described,  the  turntable  rotates  9°  after  a  single  line  scan  has 
been  completed,  and  the  process  repeats  until  the  object  has  fully  rotated  for  360°.  The  total  data 
stored  of  a  single  target  object  with  complete  360°  rotation  are  in  40x363,000  matrix  array  for 
tomographic  image  processing  at  a  later  point  in  time. 

As  discussed  in  the  previous  chapter,  the  spectral  response  of  the  transmitted  random 
noise  waveform  is  directly  related  to  the  image  quality.  If  the  oscilloscope  collects  the 
insufficient  amount  of  scattering  data,  an  unrecognizable  image  of  the  target  may  be  formed  due 
to  the  fluctuation  of  the  spectrum  response.  Although  a  long  sample  length  or  the  averaging  of 
multiple  samples  is  necessary  to  achieve  a  flat  frequency  spectrum,  the  data  size  of  40x363,000 
matrix  array  may  not  be  processed  at  the  same  time  due  to  the  memory  limitations  of  the  desktop 
computer.  For  the  computational  efficiency,  each  363,000  point  sample  is  split  into  N  segments 
for  processing  then  summed  and  averaged  to  create  an  image.  A  total  of  N  images  are 
reconstructed  for  each  rotational  angle  where  the  turntable  stops,  and  the  final  tomographic  image 
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Figure  4-4.  A  metal  cylinder  placed  on  the  turntable.  Also  shown  are  the  transmitting  and 
receiving  antenna  and  microwave  absorber  behind  the  target  to  minimize  reflections  from  the 
back. 
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is  generated  via  sum  and  average  of  N  images.  After  measurement  with  full  360  °  rotations,  a 
total  40 xA/^  final  tomographie  images  are  ereated,  and  they  are  eombined  into  one  eomplete  image. 
The  entire  experiment  proeess  is  identieal  to  the  Figure  3-2. 

Figure  4-5  shows  an  example  of  363,000  measured  amplitude  samples  in  time  domain 
and  the  frequeney  response  when  both  transmitting  and  reeeiving  antennas  are  aligned  at  the 
eenter  of  the  target  objeet.  The  measured  seattering  data  sample  in  Figure  4-5  (a)  are  split  into 
200  segments  so  that  eaeh  data  segment  eontains  1815  amplitude  samples.  For  example,  the  first 
segment  in  the  red  box  #1  shown  in  Figure  4-5 (a)  and  the  frequeney  response  are  plotted  in 
Figure  4-6.  For  eaeh  rotational  angle  0,  a  total  of  200  images  are  generated  from  a  large  sequenee 
of  a  single  noise  waveform,  and  the  final  image  of  the  objeet  at  eaeh  rotational  angular  step  is 
reeonstrueted  via  sum  and  average  of  all  200  images. 


1  2  3  4  .  197  198  199  200 


Time  (microsecond) 


Figure  4-5.  (a)  The  measured  seattering  data  when  both  Tx  and  Rx  are  loeated  at  the  eenter  of  the 
eylinder,  and  (b)  the  frequeney  speetrum  of  the  measured  data.  The  frequeney  speetrum  ranges 
are  shown  from  DC  to  8  GHz  only. 
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Figure  4-6.  (a)  The  first  segment  of  the  measured  seattering  data  shown  in  Figure  4-5(a),  and  (b) 
the  frequeney  speetrum  from  DC  to  8  GHz  only. 

4.3  Experimental  Results 

4.2.1  Tomographic  image  for  a  circular  metallic  cylinder 

A  top  view  and  side  view  of  experiment  eonfiguration  for  a  eireular  metallie  eylinder  are 
shown  in  Figure  4-7.  The  measured  radius  (r)  and  height  (/z)  of  the  eylinder  are  approximately 
10.2  em  and  40.6  em,  respeetively,  and  the  measured  distanee  {d)  from  the  eenter  of  the  eylinder 
to  the  edge  of  a  transmitting  antenna  is  approximately  66  em.  The  measured  dimensions  of  the 
eireular  metallie  eylinder  are  shown  in  Figure  4-8.  The  transmitting  antenna  (Tx)  is  fixed  at  the 
eenter  of  the  objeet,  and  the  reeeiving  antenna  (Rx)  is  moving  along  they-axis  to  eolleet  the 
seattering  data.  Moving  interval  (Ay)  for  Rx  is  1.5  em,  i.e.,  seattering  data  are  measured  and 
reeorded  for  55  different  positions  along  y-axis  from  (-66  em,  -40.5  em)  and  (-66  em,  40.5  em). 
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Figure  4-7.  Top  and  side  view  of  experiment  configuration  for  a  circular  metallic  cylinder.  The 
cylinder  is  placed  on  the  turntable  for  rotation.  Also  the  microwave  absorber  behind  the  target  is 
shown  to  minimize  reflections  from  the  back. 


Figure  4-8.  The  metallic  cylinder  used  for  the  experiment.  The  dimensions  of  the  cylinder  are 
measured  in  inches. 


As  shown  in  Figure  4-7,  the  target  object  is  placed  on  the  center  of  the  turntable  in  front 
of  a  wall  made  of  microwave  absorbent  foam.  The  turntable  is  configured  to  rotate  9°  after  a 
single  line  scan  of  55  points  along  the  linear  scanner  so  that  the  scattering  data  is  measured  for  40 
different  angles.  Figure  4-9  shows  the  complete  tomographic  image  of  the  circular  metallic 
cylinder  based  on  the  measured  scattering  data.  As  shown  in  Figure  4-9,  the  tomographic  image 
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of  the  cylinder  is  correctly  reconstructed;  however,  the  image  quality  is  not  comparable  to  the 
simulated  image  displayed  in  Figure  3-9(d).  Each  pixel  of  the  tomographic  image  displayed 
herein  represents  the  normalized  magnitude  of  scattering  object  function,  which  is  proportional  to 
the  difference  between  the  scattering  due  to  the  target  object  and  the  surrounding  medium.  Since 
the  experiment  is  performed  in  the  room  surrounded  by  concrete  walls,  the  image  is  severely 
affected  by  the  additional  scattering  data  due  to  concrete  walls  where  no  microwave  absorbers  are 
present.  Also,  unwanted  noise  from  active  RF  components  and  equipment  degrades  the  image 
quality.  However,  despite  these  unwanted  reflections,  the  cylindrical  shape  is  well  recognized. 
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Figure  4-9.  The  final  tomographic  image  of  the  metallic  cylinder  based  on  the  measured  data. 


4.2.2  Tomographic  image  for  a  square  metallic  box 

Figure  4-10  displays  a  top  view  and  side  view  of  experiment  configuration  for  a  square 
metallic  box.  The  dimension  of  the  box  is  20.3x20.3x61  cm.  The  measured  distance  (d)  from  the 
center  of  the  cylinder  to  the  edge  of  a  transmitting  antenna  is  approximately  66  cm.  The 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


97 


83 


measured  dimensions  of  the  square  metallic  box  are  shown  in  Figure  4-11.  For  the  square 
metallic  box  measurement,  moving  interval  (A;;)  for  Rx  is  1.5  cm  so  that  scattering  data  are 
measured  and  recorded  for  55  different  positions  along  j^-axis  from  (-66  cm,  -40.5  cm)  and  (-66 
cm,  40.5  cm).  The  turntable  is  configured  to  rotate  9°  after  a  single  line  scan,  and  the  scattering 
data  is  measured  for  40  different  angles. 


Figure  4-10.  Top  and  side  view  of  experiment  configuration  for  a  square  metallic  box.  The  box 
is  placed  on  the  turntable  for  rotation.  Also  the  microwave  absorber  behind  the  target  is  shown  to 
minimize  reflections  from  the  back. 


Figure  4-11.  The  metallic  box  used  for  the  experiment.  The  dimensions  of  the  box  are  measured 
in  inches. 
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As  shown  in  Figure  4-12,  the  tomographie  image  of  the  box  is  eorreetly  reeonstmeted  as 
expeeted;  however,  the  image  quality  is  not  eomparable  to  the  simulation  results  shown  in  Figure 
3-1 1(d).  Similarly,  the  image  is  severely  affeeted  by  the  additional  seattering  data  due  to  eonerete 
walls  and  unwanted  noise  from  aetive  RF  eomponents. 
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Figure  4-12.  The  final  tomographie  image  of  the  metallie  box  based  on  the  measured  data. 


4.2.3  Tomographic  image  for  an  equilateral  triangular  metallic  prism 

Figure  4-13  shows  a  top  view  and  side  view  of  experiment  eonfiguration  for  an 
equilateral  triangular  metallie  prism.  The  length  of  all  sides  of  the  triangle  is  20.3  em,  and  the 
height  Qi)  of  the  prism  is  30.5  em.  The  measured  distanee  {d)  from  the  eenter  of  the  prism  to  the 
edge  of  a  transmitting  antenna  is  approximately  66  em.  All  other  experimental  eonfigurations, 
sueh  as  the  positions  of  55  seanning  points  and  the  rotational  angle,  remain  unehanged. 

Figure  4-14  and  Figure  4-15  show  the  aetual  triangular  prism  and  the  eomplete 
tomographie  image  with  a  full  360°  rotation  of  the  objeet,  respeetively.  The  image  quality  of 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


99 


85 

tomographic  image  shown  in  Figure  4-15  is  degraded  by  the  additional  scattering  data  due  to 
concrete  walls  and  unwanted  noise  from  active  RF  components.  However,  the  cross-section 
image  of  the  prism  can  be  clearly  inferred. 


Figure  4-13.  Top  and  side  view  of  experiment  configuration  for  an  equilateral  triangular  metallic 
prism.  The  prism  is  placed  on  the  turntable  for  rotation.  Also  the  microwave  absorber  behind  the 
target  is  shown  to  minimize  reflections  from  the  back. 


Figure  4-14.  The  metallic  prism  used  for  the  experiment.  The  dimensions  of  the  prism  are 
measured  in  inches. 
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Figure  4-15.  The  final  tomographie  image  of  the  metallie  prism  based  on  the  measured  data. 

4.2.4  Tomographic  image  for  a  circular  concrete  cylinder 

A  top  view  and  side  view  of  experiment  eonfiguration  for  a  eireular  eonerete  eylinder  are 
shown  in  Figure  4-16.  The  measured  radius  (r)  and  height  {h)  of  the  eylinder  are  approximately 
7.6  em  and  30.5  em,  respeetively,  and  the  measured  distanee  {d)  from  the  eenter  of  the  eylinder  to 
the  edge  of  a  transmitting  antenna  is  approximately  66  em.  The  measured  dimensions  of  the 
eireular  eonerete  eylinder  are  shown  in  Figure  4-17.  The  transmitting  antenna  (Tx)  is  fixed  at  the 
eenter  of  the  objeet,  and  the  reeeiving  antenna  (Rx)  is  moving  along  they-axis  to  eolleet  the 
seattering  data.  Moving  interval  (Ay)  for  Rx  is  3  em,  i.e.,  seattering  data  are  measured  and 
reeorded  for  29  different  positions  along  y-axis  from  (-66  em,  -40.5  em)  and  (-66  em,  40.5  em). 
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Figure  4-16.  Top  and  side  view  of  experiment  configuration  for  a  circular  concrete  cylinder.  The 
cylinder  is  placed  on  the  turntable  for  rotation.  Also  the  microwave  absorber  behind  the  target  is 
shown  to  minimize  reflections  from  the  back.. 


Figure  4-17.  The  concrete  cylinder  used  for  the  experiment.  The  dimensions  of  the  cylinder  are 
measured  in  inches. 


The  turntable  is  configured  to  rotate  9°  after  a  single  line  scan  of  29  points  along  the 
linear  scanner  so  that  the  scattering  data  is  measured  for  40  different  angles.  Figure  4-18  shows 
the  complete  tomographic  image  of  the  circular  concrete  cylinder  based  on  the  measured 
scattering  data.  Although  the  tomographic  image  of  the  dielectric  cylinder  is  successfully 
obtained,  the  image  is  not  truly  vivid  compared  to  the  metallic  circular  cylinder,  as  expected.  The 
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measured  seattering  field  intensity  due  of  eonerete  objeet  is  not  as  strong  as  the  metallie  objeet 
ease.  Also  the  surfaee  roughness  and  the  eomposition  of  eonerete  eylinder  affeet  the  seattering 
behavior,  whieh  degrades  the  image  quality  eompared  to  the  metallie  eylinder  ease  shown  in 
Figure  4-9. 


X  [m] 

Figure  4-18.  The  final  tomographie  image  of  the  eonerete  eylinder  based  on  the  measured  data. 


4.2.5  Tomographic  Image  for  a  square  concrete  block 

Figure  4-19  displays  a  top  view  and  side  view  of  experiment  eonfiguration  for  a  square 
eonerete  bloek.  The  dimension  of  the  bloek  is  15.2x15.2x28  em.  The  measured  distanee  {d) 
from  the  eenter  of  the  eylinder  to  the  edge  of  a  transmitting  antenna  is  approximately  66  em.  The 
measured  dimensions  of  the  square  metallie  box  are  shown  in  Figure  4-20.  For  the  eonerete 
bloek  measurement,  the  transmitting  antenna  (Tx)  is  fixed  at  the  eenter  of  the  objeet,  and  only 
reeeiving  antenna  (Rx)  is  moving  along  the  y-axis  to  eolleet  the  seattering  data  for  29  different 
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scanning  points,  i.e,  Isy  =3  cm,  from  (-66  cm,  -40.5  cm)  and  (-66  cm,  40.5  cm).  The  turntable  is 
configured  to  rotate  9°  after  a  single  line  scan,  and  the  scattering  data  is  measured  for  40  different 
angles. 
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Figure  4-19.  Top  and  side  view  of  experiment  configuration  for  a  square  concrete  block.  The 
box  is  placed  on  the  turntable  for  rotation.  Also  the  microwave  absorber  behind  the  target  is 
shown  to  minimize  reflections  from  the  back. 


Figure  4-20.  The  concrete  block  used  for  the  experiment.  The  dimensions  of  the  concrete  block 
are  measured  in  inches. 
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Figure  4-21  shows  the  eomplete  tomographie  image  with  a  full  360°  rotation  of  the  objeet, 
respeetively.  As  shown  in  Figure  4-21,  the  tomographie  image  of  the  box  is  eorreetly 
reeonstrueted.  Similarly,  the  image  of  the  eonerete  box  is  affeeted  by  the  measured  field  intensity 
from  the  experiment.  The  surfaee  of  the  eonerete  bloek  is  not  truly  smooth  eompared  to  the 
eireular  eonerete  eylinder  as  shown  in  Figure  4-20. 
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Figure  4-21.  The  final  tomographie  image  of  the  eonerete  bloek  based  on  the  measured  data. 


4.4  Summary 

A  bistatie  UWB  noise  tomography  radar  system  is  implemented,  demonstrating  various 
tomographie  imaging  seenarios  using  iid  WGN  waveforms.  The  eomplete  tomographie  images 
of  the  sequentially  rotated  objeets  are  sueeessfully  reeonstrueted  with  relatively  simple  hardware 
setup.  Based  on  experimental  results  presented  in  this  ehapter,  radar  tomography  using  WGN 
waveforms  is  aeeomplished  while  maintaining  the  advantages  of  UWB  noise  radar  sueh  as  LPI 
and  high-resolution  imaging  eapability.  The  eross-seetion  of  the  metallie  objeet  is  sueeessfully 
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imaged  compared  to  the  dielectric  object  cases;  however,  more  diffraction  occurs  from  the  sharp 
edges  of  the  metallic  objects  compared  to  the  dielectric  objects.  The  reconstructed  tomographic 
images  based  on  measured  scattering  data  can  be  compared  with  those  obtained  based  on 
numerical  simulation  results,  showing  a  good  agreement  between  simulation  and  experiment 
results  for  both  metallic  and  dielectric  target  object  cases. 


4.5  References 

[1]  A.  Broquetas,  R.  Jordi,  J.  M.  Rius,  A.  R.  Elias-Fuste,  A.  Cardama  and  L.  Jofre, 
“Cylindrical  geometry:  a  further  step  in  active  microwave  tomography,”  IEEE 
Transactions  on  Microwave  Theory  and  Techniques,  vol.  39,  pp.  836-844,  1991. 

[2]  J.  C.  Bolomey  and  C.  Pichot,  "Microwave  tomography:  from  theory  to  practical  imaging 
systems,"  International  Journal  of  Imaging  Systems  and  Technology,  vol.  2,  pp.  144-156, 
1990. 

[3]  S.  Y.  Semenov,  R.  H.  Svenson,  A.  E.  Boulyshev,  A.  E.  Souvorov,  V.  Y.  Borisov,  Y. 
Sizov,  A.  N.  Starostin,  K.  R.  Dezem,  G.  P.  Tatsis  and  V.  Y.  Baranov,  “Microwave 
tomography:  two-dimensional  system  for  biological  imaging,”  IEEE  Transactions  on 
Biomedical  Engineering,  vol.  43,  pp.  869-877,  1996. 

[4]  A.  Verity,  S.  Gavirilov,  T.  Adigiizel,  I.  Voynovskyy,  G.  Yiiceer  and  A.  O.  Salman,  “C- 
band  tomography  system  for  imaging  of  cylindrical  objects,”  in  Proc.  SPIE  Conference 
on  Subsurface  Sensors  and  Applications,  Denver,  CO,  1999,  pp.  224-230. 

[5]  C.  Gilmore,  P.  Mojabi,  A.  Zakaria,  M.  Ostadrahimi,  C.  Kaye,  S.  Noghanian,  L.  Shafai,  S. 
Pistorius  and  J.  LoVetri,  “An  ultra- wideband  microwave  tomography  system:  preliminary 
results,”  in  Proc.  Annual  International  Conference  of  the  IEEE  Engineering  in  Medicine 
and  Biology  Society  (EMBC  2009),  Minneapolis,  MN,  2009,  pp.  2288-2291. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


106 


92 


[7]  M.  Z.  Abdullah,  S.  A.  Binajjaj,  T.  F.  Zanoon  and  A.  J.  Peyton,  “High-resolution  imaging 
of  dielectric  profiles  by  using  a  time-domain  ultra  wideband  radar  sensor,”  Measurement, 
vol.  44,  pp.  859-870,  2011. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


107 


93 


Chapter  5 

Conclusions  and  Future  Work 

5.1  Conclusions 

The  dissertation  has  been  concerned  with  tomographic  imaging  for  UWB  noise  radar 
using  Fourier  diffraction  theorem.  This  concluding  chapter  gives  a  summary  of  the  results  from 
previous  chapters  of  the  dissertation. 

In  Chapter  2,  the  theoretical  foundations  of  diffraction  tomography  using  a  white 
Gaussian  noise  waveform  as  a  transmitted  signal  for  the  cross-section  images  of  multiple  metallic 
target  objects  inside  a  rectangular  geometry  are  presented.  Further,  the  empirical  solution  to 
bypass  the  shortcoming  of  using  UWB  random  noise  waveform  and  the  image  quality  measures 
of  the  reconstructed  images  under  various  SNR  environments  are  discussed. 

In  Chapter  3,  the  principle  and  numerical  simulation  results  of  tomographic  imaging  of 
various  target  scenarios  using  UWB  noise  radar  are  presented.  Theoretical  analysis  of  the  image 
reconstruction  for  a  bistatic  imaging  radar  system  with  a  two-dimensional  rotating  cylinder  using 
Fourier  diffraction  theorem  under  the  assumption  of  plane  wave  illumination  is  discussed.  Also, 
the  numerical  simulation  results  of  diffraction  tomography  for  UWB  noise  radar  are  shown  in  the 
chapter. 

In  Chapter  4,  the  hardware  implementation  and  experimental  investigation  of  UWB  noise 
radar  are  proposed.  Design  considerations  and  performance  for  UWB  noise  radar  using  WGN 
waveforms  are  presented  in  the  chapter.  Further,  tomographic  image  processing  technique  for 
large  sets  of  collected  data  samples  is  proposed.  The  reconstructed  tomographic  images  based  on 
experiments  for  various  target  objects  are  presented  as  well. 
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5.2  Suggestions  for  Future  Work 

5.2.1  UWB  noise  radar  tomography  for  X-band 

The  block  diagram  for  the  transmitter  and  the  receiver  portions  of  the  system  for  3-5 
GHz  frequency  ranges  are  shown  in  Figure  4-2.  Since  the  maximum  sampling  rate  of  AWG 
(Agilent  M8190A)  is  12  GSa/s,  the  maximum  frequency  of  the  output  waveform  is  up  to  6  GHz. 
For  the  operating  frequency  beyond  6  GHz,  additional  RF  components  are  required  for  frequency 
up-conversion  of  the  bandlimited  noise  waveforms  by  heterodyning  on  the  transmitter  block. 
Also,  frequency  down-conversion  should  be  considered  for  the  receiver  design  unless  the 
sampling  rate  of  the  oscilloscope  is  twice  of  the  frequency  of  the  received  signals. 

X-band  is  a  segment  of  frequency  spectrum  ranges  between  8-12  GHz,  which  has  been 
commonly  used  for  high-resolution  imaging  and  remote  sensing  applications  in  radars  [1,  2,  3]. 
Figure  5-1  shows  a  proposed  UWB  noise  radar  configuration  for  X-band  frequency  operation, 
which  is  the  modified  hardware  design  based  on  the  hardware  configuration  shown  in  Figure  4-2. 
As  shown  in  Figure  5-1,  the  mixer  and  local  oscillator  are  added  for  a  frequency  up-conversion 
purpose.  Also,  several  power  amplifiers  may  be  cascaded  in  order  to  provide  sufficient  gain  for 
higher  frequency  ranges.  The  beam  pattern  of  both  transmitting  and  receiving  antennas  should  be 
investigated  so  that  the  target  object  is  always  in  the  main  lobe  for  all  measurements.  In  general, 
the  image  quality  of  the  tomographic  image  can  be  degraded  due  to  the  effects  of  unwanted 
additive  noise  from  additional  active  RF  components. 
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Figure  5-1.  The  proposed  design  of  transmitter  and  receiver  for  X-band  operation. 


5.2.2  Moving  average  filter  for  noise  suppression 

The  moving  average  filter  (MAF)  operates  by  averaging  a  number  of  points  from  the 
input  signal  to  produce  each  point  in  the  output  signal.  The  general  equation  of  the  moving 
average  filter  is  written  as 

j  (M-l)/2 

=  Z  (5-1) 

where  x[-  ]  is  the  input  signal,  y[-  ]  is  the  output  signal,  and  M  is  the  number  of  points  in  the 
average.  The  measured  scattered  field  data  can  be  considered  as  the  input  signal,  x[*],  and  the 
obtained  signals  after  applying  the  moving  average  filter  is  considered  as  the  output  signal,  y[-]. 
The  value  of  M,  i.e.,  the  number  of  points  in  the  average,  is  changed  accordingly  to  check  how 
the  moving  average  filter  effectively  suppresses  the  additive  noise. 
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Figure  5-2  shows  a  top  view  and  side  view  of  experiment  configuration  for  a  circular 
metallic  cylinder.  The  frequency  ranges  of  the  transmitted  noise  waveforms  are  2-4  GHz,  and 
100  iid  noise  waveforms  with  amplitude  samples  of  1001  are  transmitted  and  received.  The 
measured  radius  (r)  and  height  {h)  of  the  cylinder  are  approximately  8  cm  and  30.5  cm, 
respectively.  The  measured  distance  {d)  from  the  center  of  the  cylinder  to  the  edge  of  a 
transmitting  antenna  is  approximately  71  cm.  The  transmitting  antenna  (Tx)  is  fixed  at  the  center 
of  the  object,  and  the  receiving  antenna  (Rx)  is  moving  along  they-axis  to  collect  the  scattering 
data.  Moving  interval  (Ay)  for  Rx  is  1.5  cm,  i.e.,  scattering  data  are  measured  and  recorded  for  55 
different  positions  along  y-axis  from  (-71  cm,  -40.5  cm)  and  (-71  cm,  40.5  cm).  The  target 
object  is  not  rotated  for  the  experiment  of  moving  average  filter. 
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Figure  5-2.  Top  and  side  view  of  experiment  configuration  for  a  circular  metallic  cylinder.  The 
cylinder  is  not  rotated.  Also  the  microwave  absorber  behind  the  target  is  shown  to  minimize 
reflections  from  the  back. 


Figure  5-3  shows  the  example  of  frequency  domain  scattering  data  plots  among  1001 
collected  scattering  data  when  both  Tx  and  Rx  antennas  are  aligned  at  the  center  of  the  linear 
scanner  with  different  numbers  of  points  in  the  average:  1-pt  (no  MAF),  3-pt,  5-pt  and  7-pt  MAF. 
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Figure  5-3.  Frequency  responses  of  the  collected  scattering  data  with  (a)  1-pt  (no  MAF),  (b)  3-pt, 
(c)  5-pt,  and  (d)  7-pt  MAF  applied. 

As  shown  in  Figure  5-3,  the  smoothing  action  of  the  moving  average  filter  decreases  the 
amplitude  of  the  random  noise  in  collected  waveforms,  but  also  reduces  the  sharpness  of  the 
original  waveform  as  the  number  of  average  increase.  In  general,  the  amount  of  noise  reduction 
is  equal  to  the  square  root  of  the  number  of  points  in  the  average. 

Figure  5-4  shows  the  reconstructed  tomographic  images  based  on  4  different  points  of 
moving  average  filter  shown  in  Figure  5-3.  As  shown  in  Figure  5-4,  tomographic  images  clearly 
show  how  the  moving  average  filters  effectively  work  by  suppressing  the  random  noise  in 
collected  waveforms. 
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Figure  5-4.  Tomographic  images  based  on  measured  scattering  data  when  (a)  no  MAF,  (b)  3-pt 
MAF,  (c)  5-pt  MAF,  and  (d)  7-pt  MAF  applied. 


5.2.3  Background  subtraction  method  for  image  enhancement 

Background  subtraction  is  a  widely  used  approach  in  the  image  processing,  identifying 
target  objects  from  the  clutter  or  surrounding  environment.  Additional  data  acquisition  and  image 
processing  are  required  in  order  to  perform  the  background  subtraction.  When  the  scattering  data 
acquisition  process  with  a  target  object  is  complete,  additional  scattering  data  must  be  collected 
without  the  presence  of  the  target  object  for  obtaining  a  background  image.  The  final 
tomographic  image  is  reconstructed  based  on  the  comparison  of  two  images,  showing  the 
intensity  for  the  pixel  which  has  changed  in  two  images. 
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Two  tomographic  images  shown  in  Figure  5-5  display  the  effectiveness  of  background 
subtraction.  All  experimental  configurations  for  background  subtraction  technique  are  identical 
to  the  previous  setup  shown  in  Figure  5-2.  As  shown  in  Figure  5-5,  the  shape  of  the  target  object 
is  clearly  extracted  and  imaged,  providing  a  noticeable  difference  in  image  after  background 
subtraction  is  correctly  implemented. 


X  [m] 


X  [m] 


Figure  5-5.  Tomographic  images  based  on  measured  scattering  data  when  (a)  no  background 
subtraction,  and  (b)  background  subtraction  is  implemented. 
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ABSTRACT 


A  hardware  system  has  been  developed  to  perform  ultrawideband  (UWB)  noise  radar 
tomographic  imaging  over  the  3-5  GHz  frequency  range.  The  system  has  been  tested  on 
a  variety  of  target  objects  that  have  been  concealed  in  both  cardboard  and  wood.  The 
system  utilizes  RF  hardware  to  transmit  multiple  independent  and  identically  distributed 
(iid)  UWB  random  noise  waveforms.  A  3-5  GHz  band-limited  signal  is  generated  using 
an  arbitrary  waveform  generator  and  the  waveform  is  then  amplified  and  transmitted 
through  a  horn  antenna.  A  linear  scanner  with  a  single  antenna  is  used  in  place  of  an 
antenna  array  to  collect  backscatter.  The  backscattered  data  are  collected  from  the 
transmission  of  each  waveform  and  reconstructed  to  form  an  image.  The  images  that  result 
from  each  scan  are  averaged  to  produce  a  single  tomographic  image  of  the  target.  After 
background  subtraction,  the  scans  are  averaged  to  improve  the  image  quality.  The 
experimental  results  are  compared  to  the  theoretical  predictions.  The  system  is  able  to 
successfully  image  metallic  and  dielectric  cylinders  of  different  cross  sections. 
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1.  INTRODUCTION  AND  BACKGROUND 


Microwave  tomography  has  been  developed  and  refined  over  the  past  few  decades  to  obtain 
high-resolution  images  of  metallie  targets  and  dieleetrie  contrasts  embedded  within  a 
dielectric  medium.  The  objective  of  active  microwave  tomography  is  to  reeonstruet  the 
dieleetrie  properties  of  a  body  illuminated  with  mierowaves  from  a  measurement  of  the 
scattered  fields.  Conventional  mierowave  tomography  systems  are  based  on  illuminating 
the  body  with  a  plane  wave  and  measuring  the  scattered  fields  with  a  linear  array  of  probes. 
A  practical  requirement  when  using  this  method  is  the  need  for  mechanieal  rotation  of 
either  the  body  or  the  antenna  in  order  to  obtain  measurements  in  different  views'.  Every 
object,  when  inserted  into  an  eleetromagnetic  field,  eauses  a  well-defined  field  change,  and 
diffraction  occurs  when  the  wavelength  of  the  microwave  radiation  is  of  the  order  of  the 
dimension  of  the  objeet.  A  one-to-one  relationship  relating  the  seattered  field  to  the  objeet 
complex  permittivity  ean  be  obtained  within  the  Bom  approximation,  via  a  Fourier 
transform  using  the  so-ealled  Fourier  Diffraction  Theorem^.  Microwave  imaging  has  been 
studied  for  malignant  breast  eaneer  deteetion^'^.  Several  hardware  systems  have  been 
developed  and  implemented  over  the  years.  A  microwave  tomographic  system  consisting 
of  64  circularly  arranged  electronically  scanned  antennas,  divided  into  32  transmitters  and 
32  reeeivers,  operating  at  a  frequency  of  2.45  GHz  was  developed  for  biological  tissue 
imaging®.  A  C-band  system  operating  at  7.5  GHz  was  developed  for  imaging  weakly 
scattering  objeets,  such  as  paper  cylinders^.  More  reeently,  a  four-port  imaging  system 
using  a  veetor  network  analyzer  (VNA)  operating  over  the  7.5-12.5  GHz  frequeney  range 
was  used  to  reconstract  images  of  eontinuous  and  discrete  conducting  objects  and  a  B-52 
model  aircraft*.  A  2D  UWB  microwave  imaging  system  3-6  GHz  frequency  range  using 
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24  antenna  elements  connected  to  a  VNA  via  a  2x24  port  matrix  switch  was  able  to 
quantitatively  reconstruct  dielectric  objects^.  A  time-domain  UWB  tomographic  imaging 
system  transmitting  a  differentiated  cosine-modulated  Gaussian  pulse  with  a  carrier 
frequency  of  1  GHz  was  used  to  reconstruct  plastic  phantoms  of  different  sizes  and  shapes 
immersed  in  a  tank  filled  with  pipe  water' Tomographic  reconstruction  has  been  used  in 
many  settings,  such  as  medical  and  infrastructure""'^. 

Noise  waveforms  have  been  used  in  radar  for  many  years*"'  Noise  waveforms  provide 
several  advantages  for  covert  applications  due  to  their  low  probability  of  intercept  (LPI) 
and  low  probability  of  detection  (LPD)  characteristics*^*^.  These  attributes  come  from  the 
fact  the  random  noise  signal  changes  constantly  and  does  not  repeat^**.  Radar  tomography 
using  an  UWB  Gaussian  noise  waveform  has  been  proposed  as  a  covert  tomographic 
imaging  technique.  Using  either  resistors  or  noise  diodes  and  then  amplifying  the  thermal 
noise  it  is  possible  to  cheaply  generate  the  noise^*.  This  process  relies  on  the  properties  of 
the  noise  used  in  the  waveform.  The  method  uses  band-limited  independent  and  identically 
distributed  (iid)  white  Gaussian  Noise  (WGN)  for  the  waveform.  The  approach  utilizes 
the  fact  that  WGN  has  a  flat  frequency  spectrum.  To  actually  achieve  the  flat  frequency 
spectrum  in  hardware,  multiple  WGN  waveforms  must  be  averaged  together.  Using 
Fourier  diffraction  theory,  it  has  been  shown  that  is  indeed  possible  to  construct  a 
tomographic  image  of  the  objecfr^. 

This  thesis  discusses  a  UWB  noise  tomographic  system  specially  developed  by  us  to 
generate  tomographic  images  of  various  objects  and  presents  the  experimental  results 
from  this  system.  Chapter  2  describes  the  hardware  implementation.  Chapter  3  presents 
the  data  collection  and  processing  approach.  Chapter  4  shows  experimental  tomographic 

2 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


127 


images  on  a  variety  of  metallic  and  dielectric  targets,  both  under  open  and  concealed 
conditions.  Chapter  5  presents  conclusions  and  future  work. 
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2.  SYSTEM  DESIGN  AND  IMPLEMENTATION 


The  data  collection  system  is  controlled  using  a  single  computer.  The  computer  is  directly 
connected  to  four  devices.  The  computer  directly  controls  two  Arduinos,  an  Arbitrary 
Waveform  Generator  (AWG),  and  an  oscilloscope.  Before  the  data  collection  process 
begins,  the  computer  uploads  a  noise  waveform  to  the  AWG  using  a  MATLAB  code.  The 
AWG  continuously  transmits  that  noise  waveform  for  the  duration  of  the  test.  The  two 
Arduinos  are  used  to  control  the  position  of  the  turntable  that  the  test  object  is  placed  on 
and  the  receiving  antenna.  The  Arduino  in  control  of  the  turntable  rotates  the  object  a 
certain  number  of  degrees  that  is  selected  using  MATLAB  using  a  stepper  motor.  The 
second  Arduino  controls  a  linear  scanner  that  has  the  receiving  antenna  mounted  on  it.  The 
scanner  stops  at  each  point  that  is  designated  in  MATLAB  so  data  can  be  collected  at  that 
point.  The  receiving  antenna  is  connected  directly  to  the  oscilloscope  for  sampling.  The 
oscilloscope  collects  time-domain  data  and  then  passes  that  data  back  to  the  computer. 


2.1  Transmission  and  Reception 

The  block  diagram  for  the  transmit  and  the  receive  portions  of  the  system  are  shown  in 
Figure  1 .  The  computer  uploads  an  iid  Gaussian  Noise  Waveform  to  the  AWG.  The  AWG 
is  an  Agilent  M8190A  system  mounted  on  a  M9602A  chassis.  The  AWG  continuously 
transmits  the  noise  waveform  at  an  output  power  of  0  dBm.  The  iid  Gaussian  noise 
waveform  is  then  amplified  through  a  Mini-Circuits  ZVE-8G+  amplifier  having  a  gain  of 
30  dB  across  the  operating  band  of  3-5  GHz.  The  output  of  the  amplifier  is  at  +30  dBm. 
The  cable  between  the  AWG  and  the  amplifier  is  1.52  m  (5  ft)  long  and  has  1.6  dB  of 
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insertion  loss  over  the  3-5  GHz  frequeney  range.  The  amplified  waveform  is  transmitted 
through  an  A-Info  dual-polarization  horn  antenna.  The  dual-polarized  horn  anterma 
operates  over  the  frequency  range  of  2-18  GHz.  The  antenna  has  a  gain  of  approximately 
10  dB  over  our  3-5  GHz  operating  range.  There  is  another  1.6  dB  of  loss  between  the 
amplifier  and  the  antenna  due  to  cable  loss.  The  receiving  side  of  the  system  uses  an 
identical  antenna  to  collect  the  scattered  power.  The  receive  antenna  is  connected  to  an 
Agilent  Infmium  DS090804A  oscilloscope  which  has  a  maximum  sampling  rate  of  40 
GSa/s.  For  this  case  the  oscilloscope  sampled  at  20  GSa/s.  The  cable  between  the  receiving 
antenna  and  the  oscilloscope  is  1.83  m  (6  ft)  long  and  adds  another  2  dB  of  insertion  loss. 
The  oscilloscope  then  collects  and  directly  samples  the  reflected  waveform  and  then  sends 
it  to  the  computer  for  processing.  A  comparison  of  the  output  of  the  AWG  and  the 
reflection  from  an  aluminum  cylinder  are  presented  in  Figure  2. 

TABLE  1 :  HARDWARE  Characteristics 


3  dB  Beamwidth  (deg) 

44.75-40.07 

Antenna  Gain 

10  dB 

Amplifier  Gain 

30  dB 

Cable  Insertion  Loss 

5.2  dB 

AWG  Output  Power 

0  dBm 
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Armpliliucle  (clB5 


(a) 


(b) 


Figure  1:  Block  diagram  of  (a)  transmitter  and  (b)  receiver. 


Frequency(Ghz)  x10^ 

(a) 
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Frequency  (Ghz)  x  1 0^ 


(b) 


Frequency  (Ghz)  x  1 0^ 

(c) 


Figure  2:  (a)  Noise  waveform  taken  directly  out  of  the  AWG,  (b)  averaged  noise  waveform  taken 
directly  from  the  AWG,  and  (c)  reflected  and  averaged  waveform  off  of  a  metal  cylinder. 

2.2  Turntable  and  Linear  Scanner 

The  receiving  antenna  is  placed  on  a  linear  scanner  in  place  of  an  antenna  array,  as  shown 
in  Figure  3.  The  scanner  stops  at  55  points  (locations)  for  the  majority  of  the  collected 
data.  We  also  experimented  with  29  points  for  the  antenna  but  the  image  was  noticeably 
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worse  than  the  55-point  image,  as  shown  in  Figure  4  for  a  metallic  rectangular  box. 
Therefore,  we  used  55  points  to  collect  all  subsequent  data. 


Figure  3 :  A  metal  cylinder  placed  on  the  turntable.  Also  shown  are  the  transmitting  and  receiving 
antenna  and  microwave  absorber  behind  the  target  to  minimize  reflections  from  the  back. 
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(a) 


(b) 


Figure  4:  Tomographic  image  of  a  metallic  rectangular  box  imaged  with  (a)  55  lateral  points,  and 
(b)  29  lateral  points.  Note  that  the  imaging  with  55  points  is  much  better  and  cleaner  than  the  one 
with  29  points.  All  of  the  data  presented  in  this  document  used  55  lateral  points  for  imaging. 


The  test  objects  consisted  of  simple  geometric  shapes  of  different  sizes.  The  test  objects 
were  made  out  of  two  different  materials.  The  metal  test  objects  were  constructed  from 
sheets  of  aluminum.  The  dielectric  objects  were  made  out  of  solid  concrete.  The  test 
objects  were  placed  on  the  center  of  the  turntable  one  at  a  time.  Table  I  shows  the  test 
objects  for  which  the  images  are  presented  herein.  The  photographic  images  of  the  test 
object  are  included  in  the  appendix. 
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TABLE  2:  Test  Objects  Imaged  by  the  Tomographic  Noise  System 


DESCRIPTION 

DIMENSIONS 

Metal  Rectangular  Box  1 

20.3  cm  X  20.3  cm  x  60.96  cm  (8”  x  8”  x  24”) 

Metal  Rectangular  Box  2 

30.48  cm  X  30.48  cm  x  30.48  cm  (12”  x  12”  x 

12”) 

Metal  Cylinder  1 

15.24  cm  diameter  x  60.96  cm  height  (6” 

diameter  x  24”  height) 

Metal  Cylinder  2 

20.32  cm  diameter  x  40.64  cm  height  (8” 

diameter  x  16”  height) 

Dielectric  Rectangular  Box 

15.24  cm  X  15.24  cm  X  30.48  cm  (6”  x  6”  x  12”) 

Dielectric  Cylinder 

15.24  cm  diameter  x  30.48  cm  height  (6” 

diameter  x  12”  height) 

Dielectric  Semicircle 

7.62  cm  radius  x  30.48  cm  height  (3”  radius  x 

12”  height) 

Metal  Cylinder  of  Triangular 

Cross  Section 

15.24"  cm  side  x  30.48  cm  height  (6”  side 

X  12”  height) 
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The  tomographic  images  for  a  metal  rectangular  box  with  dimensions  20.3  cm  x  20.3  cm 
X  60.96  cm  are  shown  in  Figure  4(a).  At  each  point,  the  antenna  stops  in  order  to  ensure 
that  the  oscilloscope  has  time  to  collect  a  sufficient  amount  of  data.  After  all  55  point  of 
data  are  collected,  the  linear  scanner  moves  back  to  its  initial  position.  As  the  linear  scanner 
is  moving  back  to  the  first  point,  the  turntable  rotates  by  9  degrees.  The  linear  scanner  then 
starts  and  stops  at  the  same  55  points  as  before.  This  rotation  and  then  linear  scan  process 
are  continued  until  the  object  has  rotated  by  a  full  360  degrees.  In  the  data  presented  in 
this  document,  a  total  of  40  angular  measurements  are  made  (360^9  =  40).  The  number  of 
scanning  points  and  the  angle  of  rotation  were  chosen  as  a  tradeoff  between  data  collection 
time  and  image  resolution.  The  turntable  is  controlled  by  a  stepper  motor  and  is  able  to 
rotate  in  any  multiple  of  1 .8  degrees. 


11 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


136 


3.  DATA  COLLECTION  AND  PROCESSING 


The  target  objeet  is  placed  on  the  center  of  a  turntable  in  front  of  a  wall  made  of  microwave 
absorbent  foam.  Both  the  transmitting  and  the  receiving  antennas  point  towards  the  foam 
wall.  The  turntable  is  pictured  with  a  metal  cylinder  on  it  in  Figure  3. 

The  data  collection  process  starts  with  the  Gaussian  noise  waveform  being  uploaded  to  the 
AWG.  The  AWG  automatically  transmits  the  most  recently  uploaded  waveform.  The  noise 
waveform  is  not  identical  between  tests.  The  waveform  is  also  not  saved.  The  processing 
only  requires  that  the  transmitted  waveform  be  iid  and  white  Gaussian.  The  waveform  is 
then  amplified  and  then  transmitted  from  a  stationary  antenna  located  at  the  center  of  the 
window.  The  linear  scanner  then  moves  to  each  one  of  the  55  designated  points  stopping 
at  each  one.  Each  time  the  linear  scanner  stops,  the  oscilloscope  collects  a  total  of  363,000 
amplitude  samples.  After  a  single  line  scan  has  been  taken,  the  turntable  rotates  9  degrees 
and  the  process  repeats.  The  data  are  stored  as  they  comes  in  and  they  are  then  saved  in  a 
40x363,000  array  for  processing  at  a  later  point  in  time. 

Each  363,000  point  sample  is  split  into  200  segments  for  processing.  The  collected 
waveform  has  to  be  split  up  due  to  the  memory  limitations  of  the  computer.  Each  of  the 
200  segments  is  then  averaged  and  processed  to  create  an  image.  An  image  is  created  for 
each  angle  at  which  the  turntable  stops.  For  the  data  presented  in  this  document,  there  were 
40  images  created  for  each  object.  Using  image  processing,  these  40  images  were 
combined  into  one  final  image  such  as  those  presented  in  this  document. 

A  long  sample  length  or  the  averaging  of  multiple  samples  is  necessary  to  achieve  a  flat 
frequency  spectrum.  Fourier  diffraction  theory  is  used  to  construct  a  tomography  image  of 
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the  object.  If  too  little  data  are  used,  then  the  spectrum  is  no  longer  flat  and  it  is  not  possible 
to  create  a  recognizable  image.  In  Figure  2(a),  the  Fourier  transform  of  a  single  sample  of 
length  1815  shows  that  the  spectrum  is  not  really  flat,  as  desired.  In  Figure  2(b),  200 
segments  are  averaged  together  to  create  a  much  flatter  frequency  spectrum,  which  has 
been  shown  to  be  absolutely  necessary  to  create  a  good  image'^. 
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4.  EXPERIMENTAL  RESULTS 


The  noise  tomography  system  was  used  to  aequire  data  from  targets  with  rectangular, 
cylindrical,  and  triangular  cross  sections.  The  images  of  the  objects  were  successfully 
reconstructed.  As  expected,  there  is  less  scattered  power  from  the  dielectric  test  objects 
compared  to  the  metallic  objects. 

Figure  5  shows  the  reconstructed  tomographic  image  of  a  metal  cylinder  and  a  dielectric 
(i.e.  concrete)  cylinder.  We  note  that  the  circular  cross-section  is  clearly  seen  in  both  cases. 
The  metal  cylinder  shows  a  much  stronger  and  clearer  image  with  very  low  noise  outside 
its  physical  extent.  In  the  case  of  the  dielectric  cylinder,  the  shape  can  be  clearly  inferred; 
however,  there  is  observed  noise  outside  its  physical  extent,  owing  to  its  lower  scattered 
power. 


(a)  (b) 

Figure  5:  Tomographic  image  of  (a)  aluminum  cylinder,  and  (b)  dielectric  cylinder. 
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The  system  also  successfully  imaged  rectangular  boxes  made  of  both  types  of  materials, 
as  seen  in  Figure  6.  The  rectangular  objects  have  a  large  amount  of  diffraction  from  their 
edges,  which  can  be  seen  as  plumes.  Again,  the  cross-section  of  the  metal  box  is  well 
imaged  compared  to  the  dielectric  box;  however,  there  is  more  diffraction  from  the  edges 
of  the  metal  box  compared  to  the  dielectric  box. 

Figure  7  shows  the  image  of  a  metal  cylinder  having  a  triangular  cross-section.  Its  shape 
can  be  seen  as  also  the  wave  diffraction  from  its  edges. 


(a) 


(b) 


Figure  6:  Tomographic  image  of  (a)  aluminum  rectangular  box,  and  (b)  dielectric  rectangular 
box. 
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x(m) 

Figure  7:  Tomographic  image  of  a  triangular  metal  cylinder. 


Some  objects  were  also  imaged  when  concealed  within  a  cardboard  box  of  size  50.8  cm  x 
40.64  cm  X  54.07  cm.  Figures  8(a)  and  8(b)  show  the  images  of  metal  triangular  and 
circular  cylinders  inside  the  cardboard  box,  respectively.  By  comparing  Figures  7  and  8(a) 
and  also  Figures  5(a)  and  8(b),  it  can  be  observed  that  that  while  the  cardboard  adds  a  small 
amount  of  noise,  the  target  object  is  still  clearly  visible  and  discernible. 
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1 


Figure  8:  Tomographic  image  of  (a)  triangular  metal  cylinder  in  a  cardboard  box,  and  (b)  circular 
metal  cylinder  in  a  cardboard  box. 


Objects  were  concealed  in  a  wooden  box  with  sides  that  are  1.27  cm  thick.  The 
dimensions  of  the  box  are  38. 1  cm  x  38. 1  cm  x  38. 1  cm.  Figure  9  a-e  shows  metal 
objects  in  the  wooden  box.  The  wooden  box  adds  more  noise  but  the  target  objects  are 
still  visible  and  recognizable. 
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1 


Figure  9:  Tomographic  image  of  (a)  metal  cylinder  1,  (b)  metal  cylinder  2,  (c)  metal  rectangular 
box  1,  (d)  metal  rectangular  box  2,  (e)  metal  cylinder  with  triangular  cross  section,  (f)  dielectric 
cylinder,  and  (g)  dielectric  rectangular  box.  These  objects  are  concealed  in  a  wooden  box. 
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5.  FUTURE  WORK  AND  CONCLUSIONS 


5.1  Future  Work 

Some  work  has  been  done  towards  developing  a  higher  frequency  version  of  the  noise 
radar  tomography  system.  The  system  operates  from  8-10  GHz.  The  block  diagram  for 
the  transmit  and  the  receive  portions  of  the  system  are  shown  in  Figures  10  and  1 1 .  The 
computer  uploads  an  iid  Gaussian  Noise  Waveform  to  the  AWG.  The  AWG  is  an 
Agilent  M8190A  system  mounted  on  a  M9602A  chassis.  The  AWG  continuously 
transmits  the  noise  waveform  at  an  output  power  of  0  dBm.  The  Gaussian  waveform  is 
transmitted  at  2-4  GHz  through  two  Minicircuit  filters.  The  waveform  in  then  mixed 
with  a  6  GHz  tone  from  a  signal  generator.  3  dB  attenuators  are  located  on  each  terminal 
of  the  mixer.  The  resulting  waveform  is  8-10  GHz.  The  waveform  is  then  filtered  again 
before  being  transmitted  through  two  cascaded  amplifiers.  The  amplifiers  have  a  7  dB 
attenuator  in  between  them  to  prevent  saturation  in  the  second  amplifier.  The  first 
amplifier  is  a  Mini-Circuits  ZX60-24-S+  with  a  gain  of  24  dB.  The  second  amplifier  is  a 
Miteq  AMF-5B-8012-29P  with  a  gain  of  28  dB.  The  waveform  is  finally  transmitted 
through  an  A-info  dual-polarized  horn  antenna. 

An  identical  antenna  is  used  to  collect  the  backscatter  off  of  the  target  object.  The 
waveform  is  amplified  using  a  Minicircuits  ZX60-24-S+  with  gain  24  dB.  The  waveform 
is  then  filtered  and  passed  to  the  Agilent  Infmium  DS090804A  oscilloscope. 
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Figure  10  :  High  frequency  transmitting  system. 


Figure  11:  High  frequency  receiving  system. 
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5.2  Conclusions 


The  UWB  noise  tomography  hardware  system  designed  by  us  was  able  to  successfully 
image  a  variety  of  shapes  and  materials  using  a  relatively  simple  setup.  It  was 
demonstrated  that  an  iid  WON  waveform  can  be  used  to  correctly  image  target  objects. 
This  shows  that  radar  tomography  can  be  accomplished  while  maintaining  LPI  and  LPD 
properties.  This  document  described  a  successful  data  collection  system  for  noise  radar 
tomography.  Even  with  obstructions  such  as  a  wooden  box,  the  system  was  able  to  image 
objects  quite  clearly. 
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Appendix  A 
Images  of  Test  Objects 


Figure  12:  Metal  Rectangular  Box  1. 
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Figure  13:  Metal  Rectangular  Box  2. 
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Figure  14:  Metal  Cylinder  1. 
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Figure  15:  Metal  Cylinder  2. 
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Figure  16:  Dielectric  Rectangular  Box. 
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Figure  1 7:  Dielectric  Circular  Cylinder. 
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Figure  18:  Dielectric  Semicircular  Cylinder. 
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Figure  20:  Metal  Cylinder  of  Triangular  Cross  Section. 
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Appendix  B 


Computer  Programs 


%  Original  code  written  by  Scott  Wilson  modified  for  radar 
tomography  by  Mark  Asmuth 

function  varargout  =  ScannerGUI (varargin) 

%  SCANNERGUI  M-file  for  ScannerGUI . fig 

%  SCANNERGUI,  by  itself,  creates  a  new  SCANNERGUI  or 

raises  the  existing 
%  singleton"^ . 

o, 

'o 

%  H  =  SCANNERGUI  returns  the  handle  to  a  new  SCANNERGUI 

or  the  handle  to 

%  the  existing  singleton"^. 

o, 

'o 

%  SCANNERGUI ( 'CALLBACK' ,hObject,eventData, handles, . . .) 

calls  the  local 

%  function  named  CALLBACK  in  SCANNERGUI. M  with  the 

given  input  arguments. 

o, 

'o 

%  SCANNERGUI ( ’ Property ’ ,  *  Value  *,...)  creates  a  new 

SCANNERGUI  or  raises  the 

%  existing  singleton"^.  Starting  from  the  left, 

property  value  pairs  are 

%  applied  to  the  GUI  before  ScannerGUI_OpeningFcn  gets 

called.  An 

%  unrecognized  property  name  or  invalid  value  makes 

property  application 

%  stop.  All  inputs  are  passed  to  ScannerGUI_OpeningFcn 

via  varargin. 

o, 

'o 

%  ^See  GUI  Options  on  GUIDE’S  Tools  menu.  Choose  "GUI 

allows  only  one 

%  instance  to  run  (singleton)". 

o, 

'o 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDLES 

%  Edit  the  above  text  to  modify  the  response  to  help 
ScannerGUI 


%  Last  Modified  by  GUIDE  v2 . 5  29-Oct-2013  14:16:00 


%  Begin  initialization  code  -  DO  NOT 

gui_Singleton  =  1; 

gui_State  =  struct (’ gui_Name ’ , 

’ gui_Singleton ’ , 

’ gui_OpeningFcn ’ , 


EDIT 

mfilename,  . . . 
gui_Singleton,  . . . 
@ScannerGUI_OpeningFcn, 
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' gui_OutputFcn ' ,  @ScannerGUI_OutputFcn, 

' gui_LayoutFcn ' ,  []  ,  ... 

' gui_Callback ' ,  []); 

if  nargin  &&  ischar ( varargin { 1 } ) 

gui_State . gui_Callback  =  str2func (varargin{ 1 } ) ; 

end 

if  nargout 

[ varargout { 1 : nargout } ]  =  gui_mainf cn (gui_State, 
varargin { : } ) ; 
else 

gui_mainfcn (gui_State,  varargin { : } ) ; 

end 

%  End  initialization  code  -  DO  NOT  EDIT 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 

%  -  Executes  just  before  ScannerGUI  is  made  visible. 

function  ScannerGUI_OpeningFcn (hOb ject^  eventdata^  handles^ 
varargin) 

%  This  function  has  no  output  args^  see  OutputFcn. 

%  hObject  handle  to  figure 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 


MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  varargin  command  line  arguments  to  ScannerGUI  (see 
VARARGIN) 


%  Choose  default  command  line  output  for  ScannerGUI 

handles . output  =  hObject; 

set (gcf^  * CIoseRequestFcn * , @CIoseGuiFcn) 

handles . root_dir  =  *  C : \HoIograms2 \ * ; 

axes (handles . ho I Axe s ) ; 
axis  image; 

%  colormap (bone) ; 
axes (handles . imAxes ) ; 
axis  image; 

%  colormap (bone) ; 


try 

%  OSC  =  visa ( ’ agilent ’ ,  * GPIBO : : 8 : : INSTR * ) ; 

o, 

o 

%  handles. OSC  =  OSC; 

%  set (OSC,  ' InputBuf ferSize ' ,  leS); 

%  OSC.Timeout=10; 

%  OSC . ByteOrder= ' litt leEndian ' ; 

o, 

o 

%  fopen (OSC) ; 
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%  fprintf (OSC, ' *RST;  : AUTOSCALE ') ; 

%  idn= ' Oscilloscope ' ; 

OSC  =  visa ( ' agilent ' ,  ' TCPIP : : 130 . 203 . 252 . 125 : : INSTR' ) ; 

handles. OSC  =  OSC; 

set (OSC,  ' InputBuf ferSize ' ,  le7); 

set  (OSC, ' timeout ' ,  20); 

OSC . Timeout=20; 

OSC . ByteOrder= ' littleEndian ' ; 


f open (OSC) ; 

fprintf (OSC, ' *RST;  : AUTOSCALE ' ) ; 
idn= ' Oscilloscope ' ; 

msg  =  msgbox ([' Connected  to  ',idn],'OSC  Connected'); 
osc=  handles. OSC; 

set (handles . vnaConnStatText , 'String', 'Connected', 'Foreground 
Color ' , ' g ' ) ; 

fprintf (OSC, ' *RST;  : AUTOSCALE ' ) ; 


fprintf (osc, 
from  Channel 

%  Set  up  VNA 
fprintf  (osc, 
reference  at 
fprintf (osc. 


' : WAVEFORM: SOURCE  CHANl ' ) ; 

1 

frequency  scale 

' : TIMEBASE : REFERENCE  CENTER ' ) ; 
center 

TIMEBASE: RANGE  le-5 ' ) ; 


100  ns  full  scale 
fprintf (osc, ' : TIMEBASE : DELAY 
fprintf (osc, ' : CHANNEL 1 : RANGE 
range  to  1.6  V  full  scale 


) ; 

5' 


CHANNELl : OFFSET  O'); 
WAVEFORM: TYPE?  RAW') 


fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
order  on  the 
pause (10) ; 
uiwait (msg)  ; 
catch  ME 

msg  =  msgbox (' Unable 
! ' , 'Connection  Error' ) ; 


ACQUIRE: POINTS  363000'); 
WAVEFORM :BYTEORDER  LSBFirst ' ) 
instrument  as  well 


to  connect  to  GPIBO 


end 


uiwait (msg) 


Specify  data 


%  Display 

Time  base  to 

%  Delay  to  zero 
%  Vertical 


Set  the  byte 


INSTR 


try 

s  =  serial (' COMIO ',' BaudRate ', 9600 ) ;  %  Connect  to 

Arduino  via  Serial 
handles . s  =  s; 

set  (s,  ' InputBuf fer Size ' , 512) ; 
fopen(s);  %  Open  serial  port 

set (handles . scanConnStatlext , 'String', 'Connected', 'Foregroun 
dColor ' , ' g ' ) ; 
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fgets(s);  %  Clear  junk  data 
set_xy_position (s, 2377, 2377, 0,300) ; 
msg  =  msgbox (' Connected  to  XY  positioner  !'); 
uiwait (msg) ; 

catch  ME 

msg  =  msgbox (' Error  connecting  to  XY  positioner 
! ' , 'Connection  Error' ) ; 
uiwait (msg) ; 

end 


%  Update  handles  structure 
guidata (hOb ject ,  handles); 


%  UIWAIT  makes  ScannerGUI  wait  for  user  response  (see 
U I RESUME) 


uiwait (handles . figure 1 ) 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 

%  Close  request  callback  for  exiting  with  Close  button. 


function  CloseGuiFcn (hOb ject ,  eventdata,  handles) 


fclose (instrfind) ; 
delete (instrfind) ; 


delete (hObject) ; 
clear  all; 

%  close  all; 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


ooooooooooooooo 


%  -  Outputs  from  this  function  are  returned  to  the  command 

line . 

function  varargout  =  ScannerGUI_OutputFcn (hOb ject ^ 
eventdata^  handles ) 

%  varargout  cell  array  for  returning  output  args  (see 
VARARGOUT) ; 

%  hObject  handle  to  figure 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
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MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Get  default  command  line  output  from  handles  structure 
varargout{l}  =  handles .output; 


function  menuBar_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  menuFile  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


function  menuFile_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  menuFile  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


function  menuEdit_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  menuEdit  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


function  editParams_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  editParams  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


function  edit Save_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  editSave  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 
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%  handles  structure  with  handles  and  user  data  (see 
GUI DATA) 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 

%  -  Executes  on  slider  movement. 

function  f Slider_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  fSlider  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 


MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


%  Hints:  get (hObject ^  *  Value  * )  returns  position  of  slider 
%  get (hOb ject ^  *  Min  * )  and  get (hOb ject ^  *  Max  * )  to 

determine  range  of  slider 


try 


NofFreqs  =  length (handles . freqs ) ; 
if (NofFreqs  >  1) 

set (handles . f Slider ^  *  Visible  *  ^  *  on  * ) ; 
f_ind  =  get (handles . fSlider^  *  Value  * ) ; 

set (handles . f Slide rText^  *  String  *  ^  numSstr (handles . freqs (f_ind 
) /10-9,2) ) ; 
else 

set (handles . f Slider, 'Visible ' , ' of f ' ) ; 

end 

axes (handles . imAxes) ; 
imagesc (handles . X,  handles. y, 
abs (handles . recon_data ( : , : , f_ind) ) ) ; 
colormap (bone) 
axes (handles . holAxes ) ; 
imagesc (handles . X,  handles. y, 
abs (handles .hol_data ( : , : , f_ind) ) / (2*pi) ) ; 
colormap (bone) 


catch  ME 
end 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
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%  -  Executes  on  slider  movement. 

function  zSlider_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  zSlider  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  Value  * )  returns  position  of  slider 
%  get (hOb ject ^  ’ Min  * )  and  get (hOb ject ^  *  Max  * )  to 

determine  range  of  slider 

z  =  get (handles . zSlider^  *  Value  * ) ; 

set (handles . z Slide rText^  *  String  * , numSstr (z^  3)  ) ; 

z  o  =  z  ; 

NofFreqs  =  length (handles . freqs ) ; 
if (NofFreqs  >  1) 

set (handles . f Slider ^  *  Visible  *  ^  *  on  * ) ; 
f_ind  =  get (handles . fSlider^  *  Value  * ) ; 

else 

set (handles .fSlider^  ’Visible*^  *off*); 

end 


X  =  handles.  X  handles,  dx; 
y  =  handles,  y  handles,  dy; 

if (NofFreqs  >  1) 

handles . recon_data  =  FourierSD (handles . hol_data^  []^  []^ 

handles. kz^  zo) ; 

axes (handles . ImAxes ) ; 
image sc (handles . x^  handles . y ^ 
abs (handles . recon_data {:  r  :  r  f_ind) ) ) ; 
axes (handles . holAxes ) ; 
image sc (handles . x^  handles . y ^ 
abs (handles . hol_data {:  r  :  r  f_ind) ) ) ; 

else 

handles . recon_data  =  Fourier2D (handles . hol_data^ 
handles. kz^  zo) ; 

axes (handles . ImAxes ) ; 

imagesc (x^  y^  abs (handles . recon_data) ) ; 

axes (handles . holAxes ) ; 

imagesc (x^  y ^  abs (handles . hol_data) ) ; 


end 


%  Update  handles  structure 
guidata (hOb ject^  handles) ; 
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%  -  Executes  on  button  press  in  loadDataFileBtn . 

function  loadDataFileBtn_Callback (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  loadDataFileBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

global  dims; 

file  =  get (handles . loadDataFileText ,  *  String  * ) ; 

if (exist  (file^  ’file*)) 

Root  =  load ( file^  * -mat *) ; 

fields  =  fieldnames (Root ) ;  %  get  parent  name 

Root  =  get  field  (Root  dields  { 1 })  ; 
fields  =  fieldnames (Root ) ; 

for  ii  =  1 : length ( fields )  %  cycle  through  each 

member 

handles  =  setfield(  handles^  fields{ii}^ 
get  field  (Root  dields  {  ii  } )  )  ; 

end 

handles . hol_data  = 

handles . hol_data/max (max (max (abs (handles . hol_data) ) ) ) ; 

NofFreqs  =  length (handles . freqs ) ; 

BW  =  handles . f_stop  -  handles . f_start ; 

handles  .R_max  =  (3e 8 "^NofFreqs)  /  (4"^BW)  ; 
handles . R_res  =  handles . R_max  /  NofFreqs; 

axes (handles . holAxes ) ; 
image sc (handles . x^  handles . y ^ 
abs (handles . hol_data ( : ^  : d ) ) ) ; 

axes (handles . imAxes ) ; 
imagesc (handles . x^  handles . y ^ 
abs (handles . hol_data ( : ^  : d ) ) ) ; 

try 

set (handles . zSlider^  *Value*^  0); 
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set (handles . zSlider,  'Min',  0); 
set (handles . zSlider,  'Max',  2); 

set (handles . zSlider,  ' SliderStep ' ,  [0.005  0.1]); 


if (NofFreqs  >  1) 

set (handles .fSlider, 'Visible', 'on') ; 

set (handles . f Slider, ' Value ' , 1 ) ; 

set (handles .fSlider, 'Min',1) ; 

set (handles .fSlider, ' Max ' , NofFreqs ) ; 

set (handles .fSlider, 'SliderStep',  [1/ (NofFreqs-1 ) 

0.1]); 

temp  =  get (handles . fSlider, 'Value')  *  (NofFreqs-1)  + 

1; 

else 

set (handles .fSlider, 'Visible', 'off') ; 

end 


catch  ME 

set (handles . zSlider 
set (handles . zSlider 
set (handles . fSlider 
set (handles . fSlider 
set (handles . fSlider 

end 


'Max' , 1) ; 

'SliderStep',  [0.01  0.1]); 
'Min' , 0) ; 

' Max ' , NofFreqs ) ; 
'SliderStep',  [0.1  0.1]); 


Nx  =  length (handles . x) ; 
Ny  =  length (handles . y) ; 


[handles. kx  handles. ky  handles. kz]  = 
calc_wavenums (handles . freqs,  Ny,  Nx) ; 
zo  =  Root . zo ( 1 ) ; 


%  hol_data_cart  = 

sphere2cart_gridinterp (handles . hol_data, handles . freqs ) ; 
if (NofFreqs  >  1) 

handles . recon_data  =  FourierSD (handles . hol_data, 
handles. kx,  handles. ky,  handles. kz,  zo) ; 
else 

handles . recon_data  =  Fourier2D (handles . hol_data, 
handles. kz,  zo) ; 
end 


%  handles . est_targ_dist  = 

est imate_target_distance (handles . recon_data) ; 
else 

msgbox (' Please  select  a  data  f ile .',' Unable  to  load 
file ' ) ; 
end 


%  Update  handles  structure 
guidata (hOb ject ,  handles); 
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%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  zSlider_CreateFcn (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  zSlider  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  slider  controls  usually  have  a  light  gray 
background . 

if  is equal (get (hOb ject  ^  *  BackgroundColor  * ) ^ 
get  (0^  ’ def aultUicontrolBackgroundColor * ) ) 

set (hOb ject ,  *  BackgroundColor  *  ^  . 9  .9  . 9 ] ) ; 

end 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  f Slider_CreateFcn (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  fSlider  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  slider  controls  usually  have  a  light  gray 
background . 

if  is equal (get (hObject  ^  ’ BackgroundColor  * ) ^ 
get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hObject ,  ’ BackgroundColor  *  ^  . 9  .9  . 9 ] ) ; 

end 


%  -  Executes  on  button  press  in  renderSdBtn. 

function  render3dBtn_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  renderSdBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 
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figure; 

isosurface (abs (handles . recon_data) ) ; 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 

%  -  Executes  on  button  press  in  startScanBtn . 

function  start ScanBtn_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  startScanBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 


MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


pause  (2 ) ; 


x_max  =  str 2 double (get (handles . scanWidthText  ^  *  String  * ) ) ; 
y_max  =  str 2 double (get (handles . scanWidthText ,  *  String  * ) ) ; 
angle  =0;  %  =  90;  %  H-pol  or  V-pol 

m  =  str 2 double (get (handles . Xs ample s Text ,  *  String  * ) ) ; 
n  =  str 2 double (get (handles . Ys ample s Text ,  *  String  * ) ) ; 


handles . x_max  =  x_max; 
handles . y_max  =  y_max; 
handles. m  =  m; 
handles. n  =  n; 

handles. X  =  linspace ( 0 ^ x_max^ m) ; 
handles. y  =  linspace ( 0 ^ y_max^ n) ; 

handles . angle  =  angle; 

axes (handles . holAxes ) ; 
image sc (zeros (m^  n) ) ; 

[file  p_data  p_ind]  =  create_position_data (m^  n^  x_max^ 
y_max^  angle) ; 

handles . p_ind  =  p_ind;  %  indices  for  scan  position 
handles . a_data  =  check_position_data (p_data) ; 

handles . hol_data  =  run_scan (hOb ject ^  handles); 

handles . lastSavePath  =  save_hol_data (handles ) ; 

axes (handles . ImAxes ) ; 


%  Update  handles  structure 
guidata (hOb ject^  handles) ; 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 
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handles ) 


%  -  Executes  on  button  press  in  browseBtn. 

function  pushbutton4_Callback (hOb ject^  eventdata^ 

%  hObject  handle  to  browseBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUiDATA) 


function  f ileTextbox_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  fileTextbox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUiDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
fileTextbox  as  text 

%  St r 2 double (get (hObject  ^  *  String  * ) )  returns  contents 

of  fileTextbox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  f ileTextbox_CreateFcn (hObject^  eventdata^  handles) 

%  hObject  handle  to  fileTextbox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  iSPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hObject ^  * BackgroundColor *) A 

get  (0^  ’ def aultUicontrolBackgroundColor * ) ) 

set  (hObject  ^  *  Backgrounded or  * ,  ' white  * ) ; 

end 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

9'9'9'9'9'9'9'9'9'9'9'9^9^9^9^ 

ooooooooooooooo 


-  Executes  on  button  press  in  browseBtn. 
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function  browseBtn_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  browseBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

[filename^  filepath]  =  uiget file  (*  "^ . mat  *  ^  *  Select  a 
file  *  ^  *  C : \Holograms\ * ) ; 
if (filename  ~=  0) 

set (handles . loadDataFileText ,  *  String  * ,  [filepath  filename ] ) ; 
end 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooooooooooooooo 


function  ScanWidthBox_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  ScanWidthBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
ScanWidthBox  as  text 

%  St r 2 double (get (hObject ,  *  String  * ) )  returns  contents 

of  ScanWidthBox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  ScanWidthBox_CreateFcn (hObject^  eventdata^  handles) 
%  hObject  handle  to  ScanWidthBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hObject ^  * BackgroundColor *) A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hObject ,  *  Backgrounded or  *  ^  *  white  * ) ; 

end 
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function  ScanHeightBox_Callback (hOb ject^  eventdata^  handles) 
%  hObject  handle  to  ScanHeightBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  ’ String  * )  returns  contents  of 
ScanHeightBox  as  text 

%  St r 2 double (get (hOb ject  ^  *  String  * ) )  returns  contents 

of  ScanHeightBox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  ScanHeightBox_CreateFcn (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  ScanHeightBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  8c8c  isequal  (get  (hObject  ^  *  BackgroundColor  *)  A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hObject  ^  *  BackgroundColor  *  ^  *  white  * ) ; 

end 


function  XsamplesBox_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  XsamplesBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
XsamplesBox  as  text 

%  St r 2 double (get (hObject  ^  *  String  * ) )  returns  contents 

of  XsamplesBox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  XsamplesBox_CreateFcn (hObject^  eventdata^  handles) 

%  hObject  handle  to  XsamplesBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
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CreateFcns  called 


%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hOb ject ^  ’ BackgroundColor * ) A 

get  (0^  ’ def aultUicontrolBackgroundColor * ) ) 

set (hOb ject  ^  *  BackgroundColor  *  ^  *  white ’ ) ; 

end 


function  YsamplesBox_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  YsamplesBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hOb ject ^  *  String  * )  returns  contents  of 
YsamplesBox  as  text 

%  St redouble (get (hOb ject  ^  *  String  * ) )  returns  contents 

of  YsamplesBox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  YsamplesBox_CreateFcn (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  YsamplesBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hObject ^  ’ BackgroundColor *) A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hObject  ^  *  BackgroundColor  * ,  ' white  * ) ; 

end 


%  -  Executes  on  button  press  in  SingleFreqRadio . 

function  SingleFreqRadio_Callback (hObject^  eventdata^ 
handles ) 

%  hObject  handle  to  SingleFreqRadio  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hint:  get (hObject ^  *  Value  * )  returns  toggle  state  of 
SingleFreqRadio 
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%  -  Executes  on  button  press  in  MultiFreqRadio . 

function  Mult iFreqRadio_Callback (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  MultiFreqRadio  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hint:  get (hObject ^  *  Value  * )  returns  toggle  state  of 
MultiFreqRadio 


function  stopFreqBox_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  stopFreqBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
StopFreqBox  as  text 

%  St r 2 double (get (hObject  ^  *  String  * ) )  returns  contents 

of  StopFreqBox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  stopFreqBox_CreateFcn (hObject^  eventdata^  handles) 

%  hObject  handle  to  stopFreqBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  8c8c  isequal  (get  (hObject  ^  *  BackgroundColor  *)  A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hObject  ^  *  Backgrounded or  * ,  ' white  * ) ; 

end 


function  startFreqBox_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  startFreqBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 
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%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
startFreqBox  as  text 

%  St r 2 double (get (hOb ject  ^  *  String  * ) )  returns  contents 

of  StartFreqBox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  startFreqBox_CreateFcn (hOb ject^  eventdata^  handles) 
%  hObject  handle  to  startFreqBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  8c8c  isequal  (get  (hOb  ject  ^  *  BackgroundColor  *  )  A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hOb ject  ^  *  Backgrounded or  *  ^  *  white  * ) ; 

end 


function  FsamplesBox_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  FsamplesBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
FsamplesBox  as  text 

%  St r 2 double (get (hObject  ^  *  String  * ) )  returns  contents 

of  FsamplesBox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  FsamplesBox_CreateFcn (hObject^  eventdata^  handles) 

%  hObject  handle  to  FsamplesBox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hObject ^  * BackgroundColor *) A 

get  (0^  ’ def aultUicontrolBackgroundColor * ) ) 

set (hObject  ^  *  Backgrounded or  * ,  ' white  * ) ; 
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end 


%  -  Executes  on  button  press  in  abortScanBtn . 

function  abort ScanBtn_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  abortScanBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


function  scanWidthText_Callback (hOb ject^  eventdata^  handles) 
%  hObject  handle  to  scanWidthText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
scanWidthText  as  text 

%  St r 2 double (get (hObject  ^  *  String  * ) )  returns  contents 

of  scanWidthText  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  scanWidthText_CreateFcn (hObject^  eventdata^ 
handles ) 

%  hObject  handle  to  scanWidthText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hObject ^  * BackgroundColor *) A 

get  (0^  ’ def aultUicontrolBackgroundColor * ) ) 

set  (hObject  ^  *  Backgrounded or  * ,  ' white  * ) ; 

end 


function  scanHeightText_Callback (hObject^  eventdata^ 
handles ) 

%  hObject  handle  to  scanHeightText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 
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%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
scanHeightText  as  text 

%  St r 2 double (get (hOb ject  r  *  String  * ) )  returns  contents 

of  scanHeightText  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  scanHeightText_CreateFcn (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  scanHeightText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hOb ject ^  * BackgroundColor * ) A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hOb ject  ^  *  Backgrounded or  *  ^  *  white  * ) ; 

end 


function  XsamplesText_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  XsamplesText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
XsamplesText  as  text 

%  St r 2 double (get (hObject  r  *  String  * ) )  returns  contents 

of  XsamplesText  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  XsamplesText_CreateFcn (hObject^  eventdata^  handles) 
%  hObject  handle  to  XsamplesText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hObject ^  * BackgroundColor *) A 

get  (0^  ’ def aultUicontrolBackgroundColor * ) ) 

set (hObject  ^  *  Backgrounded or  *  ^  *  white  * ) ; 
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end 


function  YsamplesText_Callback (hOb ject^  eventdata^  handles) 

%  hObject  handle  to  YsamplesText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  ’ String  * )  returns  contents  of 
YsamplesText  as  text 

%  St r 2 double (get (hOb ject  ^  *  String  * ) )  returns  contents 

of  YsamplesText  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  YsamplesText_CreateFcn (hOb ject^  eventdata^  handles) 
%  hObject  handle  to  YsamplesText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  8c8c  isequal  (get  (hObject  ^  *  BackgroundColor  *)  A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hObject  ^  *  Backgrounded or  * ,  ' white  * ) ; 

end 


function  edit 5_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  fileTextbox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
fileTextbox  as  text 

%  St r 2 double (get (hObject  ^  *  String  * ) )  returns  contents 

of  fileTextbox  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  edit 5_CreateFcn (hObject^  eventdata^  handles) 

%  hObject  handle  to  fileTextbox  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
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%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hOb ject ^  * BackgroundColor * ) A 

get  (0^  ’ def aultUicontrolBackgroundColor * ) ) 

set (hOb ject  ^  *  BackgroundColor  *  ^  *  white  * ) ; 

end 


%  -  Executes  on  button  press  in  singleFreqRadio . 

function  singleFreqRadio_Callback (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  singleFreqRadio  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hint:  get (hOb ject ^  *  Value  * )  returns  toggle  state  of 
singleFreqRadio 


%  -  Executes  on  button  press  in  multiFreqRadio . 

function  mult iFreqRadio_Callback (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  multiFreqRadio  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hint:  get (hObject ^  *  Value  * )  returns  toggle  state  of 
multiFreqRadio 


function  startFreqText_Callback (hObject^  eventdata^  handles) 
%  hObject  handle  to  startFreqText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
StartFreqText  as  text 

%  St redouble (get (hObject  ^  *  String  * ) )  returns  contents 

of  StartFreqText  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 
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properties . 

function  startFreqText_CreateFcn (hOb ject ,  eventdata, 
handles ) 

%  hObject  handle  to  startFreqText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hOb ject BackgroundColor ') , 
get  (0,  ' def aultUicontrolBackgroundColor ' ) ) 

set (hOb ject , ' BackgroundColor ' , ' white ' ) ; 

end 


function  stopFreqText_Callback (hOb ject ,  eventdata,  handles) 

%  hObject  handle  to  stopFreqText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUI DATA) 

%  Hints:  get (hObject String ' )  returns  contents  of 
StopFreqText  as  text 

%  str2double (get (hObject String ') )  returns  contents 

of  StopFreqText  as  a  double 


%  -  Executes  during  object  creation,  after  setting  all 

properties . 

function  stopFreqText_CreateFcn (hObject ,  eventdata,  handles) 
%  hObject  handle  to  stopFreqText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  &&  isequal (get (hObject ,' BackgroundColor ')  , 
get  (0,  ' def aultUicontrolBackgroundColor ' ) ) 

set (hObject, ' BackgroundColor ' , ' white ' ) ; 

end 


function  FsamplesText_Callback (hObject ,  eventdata,  handles) 

%  hObject  handle  to  FsamplesText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 

52 


DISTRIBUTION  A:  Distribution  approve(d  for  public  release. 


177 


MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  ’ String  * )  returns  contents  of 
FsamplesText  as  text 

%  St r 2 double (get (hOb ject  ^  *  String  * ) )  returns  contents 

of  FsamplesText  as  a  double 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  FsamplesText_CreateFcn (hOb ject^  eventdata^  handles) 
%  hObject  handle  to  FsamplesText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  8c8c  isequal  (get  (hOb  ject  ^  *  BackgroundColor  *  )  A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hOb ject  ^  *  BackgroundColor  * ,  ' white  * ) ; 

end 


%  -  Executes  during  object  creation^  after  setting  all 

properties . 

function  loadDataFileText_CreateFcn (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  loadDataFileText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  empty  -  handles  not  created  until  after  all 
CreateFcns  called 

%  Hint:  edit  controls  usually  have  a  white  background  on 
Windows . 

%  See  ISPC  and  COMPUTER. 

if  ispc  8c8c  isequal  (get  (hObject  ^  *  BackgroundColor  *)  A 

get  (0^  * def aultUicontrolBackgroundColor * ) ) 

set (hObject  ^  *  BackgroundColor  * ,  ' white  * ) ; 

end 


%  -  Executes  on  button  press  in  browseBtn. 

function  pushbutton6_Callback (hObject^  eventdata^  handles) 

%  hObject  handle  to  browseBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 
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function  loadDataFileText_Callback (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  loadDataFileText  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 

%  Hints:  get (hObject ^  *  String  * )  returns  contents  of 
loadDataFileText  as  text 

%  St r 2 double (get (hOb ject  ^  *  String  * ) )  returns  contents 

of  loadDataFileText  as  a  double 


% - If  Enable  ==  *  on  * ^  executes  on  mouse  press  in  5  pixel 

border . 

%  -  Otherwise^  executes  on  mouse  press  in  5  pixel  border 

or  over  renderSdBtn. 

function  render3dBtn_ButtonDownFcn (hOb ject^  eventdata^ 
handles ) 

%  hObject  handle  to  renderSdBtn  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of 
MAT LAB 

%  handles  structure  with  handles  and  user  data  (see 
GUIDATA) 


figure; 

isosurface (handles . recon_data) ; 
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function  set_xy_position (s,  x,  y,  z,  delay) 


stepsPerInch  =  520;  %Define  scanner  resolution 
stepsPerDegree  =  67; 

%f get s ( s ) ; 

%sit  and  wait  for  device  to  say  it  has  reset  to  (0,0,0) 
%while ( s . BytesAvailable  ==  0)  ;  end 

%fgets(s);  %  clear  serial  buffer 

x_value  =  num2str(x) ;  %Write  the  x-value  of  that  position  as 
a  string. 

y_value  =  num2str(y) ;  %Write  the  y-value  of  that  position  as 
a  string. 

z_value  =  num2str(z) ;  %Write  the  z-value  of  that  position  as 
a  string. 

delay_value  =  num2str (delay ) ;  %Write  the  y-value  of  that 
position  as  a  string. 

string2Arduino  =  strcat (x_value, ' ; ' , y_value, ' ; ' , z_value, ' ; ' , 
delay_value, ' ; ' ) ;  %Combine  x,  y  and  z  together,  separated  by 

;  '  s  . 

query ( s , St ring2Arduino ) ; 

%disp (string2Arduino) ; 

disp ([' Current  Point:  x  =  ',  num2str ( (x/stepsPerInch) ) ,  ',  y 

=  ',  num2str ( (y/stepsPerInch) ) ,  ',  z  =  ', 

num2str ( (z/stepsPerDegree) ) ] ) ; 
disp ( '  ' ) ; 

%This  sits  and  wiats  for  data  to  be  available. . . 

%  while ( s . BytesAvailable  ==  0)  %While  the  buffer  is  empty 
(no  data  has  been  sent  back) 

%No  operation  NOOP 

%  end  %Cont inuously  wait  (infinite  loop) . 

%  fgets(s);  %Gets  ok  signal  from  Arduino  to  proceed. 

end 
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function  [ saveFilePath]  =  save_hol_data (h) 

root_dir  =  h.root_dir; 
hol_data  =  h.hol_data; 

saveYesOrNo  =  questdlg ( *  Keep  data? ’Save 
data ’ , ’Yes’, ’ No ’ , ’ Yes ’ ) ; 

if (  strcmp ( saveYesOrNo, ’ Yes ’ )  ) 

%  Create  data  structure  for  saving 
%  Root . 

Root . m  =  h . m; 

Root . n  =  h . n; 

Root.x_max  =  h.x_max; 

Root.y_max  =  h.y_max; 

Root .X  =  h . x; 

Root.y  =  h.y; 

Root.dx  =  h . x_max  /  h.m; 

Root.dy  =  h . y_max  /  h.n; 

Root . zo  =  0 ; 

Root . img  =0;  %  needed  to  run  through  simulation  model 

Root. angle  =  h. angle; 

Root.p_ind  =  h.p_ind; 

Root . hol_data  =  h.hol_data; 


p  =  0; 

saveFilePath  =  [root_dir  ’hol_scan_’  numPstr (p, ’ %04d ’ ) 

’ . mat ’ ] ; 

while (exist (saveFilePath, ’ file ’ )  ) 

p  =  p  +  1; 

saveFilePath  =  [root_dir  'hol_scan_' 
num2str (p, ' %04d' )  '.mat']; 

end 

save (saveFilePath, ' Root ' , ' -mat ' ) ; 

end 

uiwait (msgbox ( [ ' Data  saved  as:  '  saveFilePath] ,' Data 
Saved ' , ' modal ' ) ) ; 

end 

function  [hol_data]  =  run_scan (hOb ject ,  handles) 
root_dir_turntable  =  ' C : \Holograms3\ ' ; 

%  osc=  handles . OSC; 

%  fprintf (OSC, ' *RST;  : AUTOSCALE ') ; 
osc=  handles. OSC; 
s  =  handles. s; 
a_data  =  handles . a_data; 


56 


DISTRIBUTION  A:  Distribution  approve(d  for  public  release. 


181 


m  =  handles. m; 
n  =  handles. n; 

s2  =  serial (' COM4 BaudRate 9600 ) ;  %  Connect  to 

Arduino  via  Serial 


set  ( s2 ,  ' InputBuf  f erSize ' , 512 ) ; 
fopen(s2);  %  Open  serial  port 
x=-400; 

delay=20; 

x_value  =  num2str(x) ;  %Write  the  x-value  of  that  position  as 
a  string. 

%y_value  =  num2str (y ) ;  %Write  the  y-value  of  that  position 
as  a  string. 

%z_value  =  num2str ( z ) ;  %Write  the  z-value  of  that  position 
as  a  string. 

delay_value  =  num2str (delay ) ;  %Write  the  y-value  of  that 
position  as  a  string. 

string2Arduino  =  strcat (x_value, ' ; ' ,  delay_value, ' ; ' ) ; 
%Combine  x,  y  and  z  together,  separated  by  ; ' s . 

%query ( s2 , St ring2Arduino ) ; 
x=400; 

x_value  =  num2str(x) ;  %Write  the  x-value  of  that  position  as 
a  string. 

%y_value  =  num2str (y ) ;  %Write  the  y-value  of  that  position 
as  a  string. 

%z_value  =  num2str ( z ) ;  %Write  the  z-value  of  that  position 
as  a  string. 

delay_value  =  num2str (delay ) ;  %Write  the  y-value  of  that 
position  as  a  string. 

string2Arduino  =  strcat (x_value, ' ; ' ,  delay_value, ' ; ' ) ; 
%Combine  x,  y  and  z  together,  separated  by  ; ' s . 

%query ( s2 , st ring2Arduino ) ; 

%  fprintf (osc, ' *RST;  : AUTOSCALE ') ; 

o, 

o 

%  fprintf (osc, WAVEFORM: SOURCE  CHANl ' ) ;  %  Specify 

data  from  Channel  1 

o, 

o 

%  %  Set  up  VNA  frequency  scale 

%  fprintf (osc, ' iTIMEBASE: REFERENCE  CENTER');  %  Display 
reference  at  center 

%  fprintf (osc, ' iTIMEBASE: RANGE  75E-9');  %  Time  base  to 

100  ns  full  scale 

%  fprintf (osc, ': TIMEBASE : DELAY  O');  %  Delay  to 

zero 

%  fprintf (osc, ' iCHANNELl : RANGE  1.6');  %  Vertical 

range  to  1.6  V  full  scale 
%  fprintf (osc, ' iCHANNELl : OFFSET  O'); 

%  fprintf (osc, ': WAVEFORM: TYPE?  RAW'); 

%  fprintf (osc, ': WAVEFORM iBYTEORDER  LSBFirst ' ) ;  %  Set  the 

byte  order  on  the  instrument  as  well 
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o, 

o 

o,  o, 
o  o 


%  if(length(f)  >  1) 


o,  o, 
o  o 

fprintf (vna, ['STAR  '  num2str (min ( f ) 

*10-9) ] 

)  ; 

o, 

o 

Set 

o,  o, 
o  o 

start  freq 

fprintf (vna, ['STOP  '  num2str (max ( f ) 

*10-9) ] 

) ; 

o, 

o 

Set 

o, 

o 

start  freq 

fprintf (osc,  ' LINFREQ;  ' )  ; 

%  Set  to 

linear-sweep  mode 

%  %  fprint f (vna^  [ * POIN  *  num2str (length (f) )  *;*]  ); 

%  %  fprintf (vna^  *MANTR1G;  * ) ;  %  Enable  manual 

sweep  trigger 
%  else 

%  %  fprintf (vna, ’ SENS : SWEEP : TYPE  CW’);  %  Put  analyzer 

in  CW  mode 

%  fprintf  (osc^  [*  CWFREQ  *  numPstr  (f  "^10^9)  ]  )  ; 

%  end 

o, 

o 

%  Get  the  nr  of  trace  points 
%  Get  the  reference  level 
%  fprintf (osc, ' : WAVEFORM : DATA? ' ) ; 

%  Wavef ormRawData  =  fscanf  (osc) ; 

%  temp  =  str2num (Wavef ormRawData) ; 

%  holLength=length ( temp) 
rotations=10; 
scan_iterations=l ; 

%  meas_data  =  zeros (  length (a_data) , length ( f)  ) ; 
fprintf (osc, ' : WAVEFORM : DATA? ' ) ; 

Wavef ormRawData  =  fscanf  (osc) ; 
temp  =  str2num (Wavef ormRawData) ; 

hol_data  =  zeros (n,  m, scan_iterations,  length (temp) ) ; 
axes (handles . holAxes ) ; 

x_ind  =  handles . p_ind (:, 2 ) ;  %  indices  switched  for 

imagesc ( ) . . . 

y_ind  =  handles . p_ind (:, 1 ) ; 

X  =  a_dat  a ( 1 , 1 ) ; 
y  =  a_data (1,2) ; 
angle  =  a_data ( 1 , 3 ) ; 
delay  =  a_data(l,4); 
set_xy_position (s, x, y, angle, delay) ; 
pause 

X  =  a_dat  a ( 1 , 1 ) ; 
y  =  a_data (1,2) ; 
angle  =  a_data(l,3); 
delay  =  a_data(l,4); 

set_xy_position (s, x, y, angle, delay) ; 
pause ( 10 ) 
for  c=l : rotations 

for  i  =  1 : length (a_dat a) 
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Wait  till 


X  =  a_data (i, 1) ; 
y  =  a_data (i, 2) ; 
angle  =  a_data(i,3); 
delay  =  a_data(i,4); 

set_xy_position (s, x, y, angle, delay) ; 
pause ( 1 ) ; 

%query (pxa, ' OPC? ; SING; ' ) 
op.  complete 

%  fprintf (osc, ' : WAVEFORM : DATA? ' ) ; 

%  Wavef ormRawData  =  fscanf  (osc) ; 

%  length ( str2num (Wavef ormRawData) ) 

%  %  hol_data (x_ind (i ) , y_ind (i ) ,  :)  = 

str2num (Wavef ormRawData) ; 

o,  o, 
o  o 

%  for  j=l : scan_iterations 
%  fprintf (osc, ' : WAVEFORM : DATA? ' ) ; 

%  Wavef ormRawData  =  fscanf  (osc) ; 

%  temp  =  str2num (Wavef ormRawData) ; 

%  hol_data (x_ind ( i ) , y_ind (i ) , c, j , : )  =  temp; 

%  pause (2 ) 

%  end 

for  j=l : scan_iterations 
fprintf (osc, ' : WAVEFORM : DATA? ' ) ; 

Wavef ormRawData  =  fscanf  (osc) ; 
temp  =  str2num (Wavef ormRawData) ; 
hol_data (x_ind ( i ) , y_ind ( i ) , j , : )  =  temp; 
pause (2 ) 
end 

%  figure ( 1 ) 

%  imagesc (abs (hol_data ( : , : , c, 1 ) ) ) ; 

%  axis  image; 

%  guidata (hOb ject ,  handles . holAxes ) ; 

pause  (0.01) ; 

data=hol_data (x_ind ( i ) , y_ind ( i ) , 1 , : ) ; 
data=squeeze (data) ; 
dt=10*10-"-9; 
f s=l / dt ; 

N=length (temp) ; 
df=l/ (N*dt) ; 
freqs=(0:N/2) *40000000; 

Fx=fft (temp) ; 

Fx2=Fx (1 : (N/2+1) ) ; 

%  figure (2 ) 

%  stem ( freqs ( 1 : 500 ), Fx2 ( 1 : 500 ) ) 

%  figure ( 3 ) 

%  stem (abs (Fx2 ( 1 : 31 98 ) ) ) 
mvariable=abs (Fx2 (1:3198) ) ; 
for  k=3 : (length (mvariable) -3) 

meanvariable (k) = (mvariable (k-2) Imvariable (k- 
1) Tmvariable (k) Imvariable (k+1) Imvariable (k+2) ) /5; 
end 
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meanvariable ( 1 ) =0 ; 
meanvariable ( 2 ) =0 ; 

meanvariable=20*logl0 (meanvariable) ; 

%f igure ( 4 ) 

%plot (meanvariable (1:500) ) 

%  activate  the  grid  lines 

%  Disconnect  an  clean  up 

end 

f close ( s ) ; 

s  =  serial (' COMIO BaudRate 9600 ) ;  %  Connect  to 

Arduino  via  Serial 

set  (s,  ' InputBuf ferSize ' , 512) ; 
fopen(s);  %  Open  serial  port 

fgets(s);  %  Clear  junk  data 

x=4  0; 

x_value  =  num2str(x) ;  %Write  the  x-value  of  that  position  as 
a  string. 

%y_value  =  num2str(y);  %Write  the  y-value  of  that  position 
as  a  string. 

%z_value  =  num2str(z);  %Write  the  z-value  of  that  position 
as  a  string. 

delay_value  =  num2str (delay ) ;  %Write  the  y-value  of  that 
position  as  a  string. 

string2Arduino  =  strcat (x_value, ' ; ' ,  delay_value, ' ; ' ) ; 
%Combine  x,  y  and  z  together,  separated  by  ; 's. 
query ( s2 , St ring2Arduino ) ; 
pause ( 120 ) 

p  =  0; 

saveFilePath  =  [root_dir_turntable  *hol_scan_* 
numPstr (p^  *  %04d’ )  *  .mat  * ] ; 

while (exist (saveFilePath^  *  file  * )  ) 

p  =  p  +  1; 

saveFilePath  =  [ root_dir_turntable  'hol_scan_' 
num2str (p, ' %04d' )  '.mat']; 

end 


save (saveFilePath, ' hol_data ' , ' -mat ' ) ; 
end 

f close  ( s2 ) ; 

guidata (hOb ject ,  handles); 
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function  [file,  pos_data,  pos_ind]  = 

create_pos_data (NofXpts,  NofYpts,  Xrange,  Yrange,  angle) 

delay  =  300;  %  delay  between  measurements,  milliseconds 

X  =  linspace ( 0 , Xrange, NofXpts ) ; 

Y  =  linspace ( 0 , Yrange, NofYpts ) ; 

file  =  'xy_pos_data.txt'; 

delete (file) ; 

fid  =  f open (file, ' w ' ) ; 

pos_data  =  zeros ( length (X) *length (Y) ,  4); 

pos_ind  =  zeros ( length (X) *length (Y) ,  3); 

n  =  1; 

for  i  =  1 : length ( angle ) 
for  j  =  1: length (Y) 
if(  mod(j,2)  ) 

for  k  =  1 : 1 : length (X) 
pos_ind(n,l)  =  k; 
pos_ind(n,2)  =  j; 
pos_ind(n,3)  =  i; 
pos_data (n, 1 )  =  X(k); 

pos_data (n, 2 )  =  Y(j); 

pos_data (n, 3 )  =  angle; 
pos_data (n, 4 )  =  delay; 

fprintf (fid, ' %2 . 3f  %2.3f  %2.3f  %4.0f 
\n',X(k),Y(j), angle, delay) ; 

n  =  n  +  1 ; 

end 

else 

for  k  =  length (X) :-l:l 
pos_ind(n,l)  =  k; 
pos_ind(n,2)  =  j; 
pos_ind(n,3)  =  i; 
pos_data (n, 1 )  =  X(k); 

pos_data (n, 2 )  =  Y(j); 

pos_data (n, 3 )  =  angle; 
pos_data (n, 4 )  =  delay; 
fprintf (fid, ' %2 . 3f  %2.3f  %2.3f  %4.0f 
\n',X(k),Y(j), angle, delay) ; 

n  =  n  +  1 ; 

end 

end 

end 

end 

fclose (fid) ; 
end 
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function  [Arduino_data]  =  check_position_data (position_data) 

%Define  the  limits  in  inches  of  the  scanner 
x_limit  =  35.6; 
y_limit  =  35.6; 
z_limit  =  360; 

%Find  the  indices  of  values  that  are  outside  of  the  scanner 
bounds 

ixl  =  find (posit ion_data (: ^ 1 )  >  x_limit) ; 

iyl  =  find (posit ion_data (: ^ 2 )  >  y_limit); 

izl  =  find (posit ion_data (: ^ 3 )  >  z_limit) ; 

%make  sure  the  value  are  positive. . . . 
ix2  =  find (posit ion_data (: ^ 1 )  <  0); 

iy2  =  f ind (position_data ( : ^ 2 )  <  0); 

iz2  =  f ind (position_data ( : ^ 3 )  <  0); 

%put  the  bad  data  indicies  together 
ix  =  [ ixl  ix2 ] ; 

iy  =  [iyl  iy2] ; 
iz  =  [izl  iz2]; 

%Find  the  x  and  y  values  that  are  outside  of  the  scanner 
bounds 

bad_x_vals  =  posit ion_data (ix^ 1) ; 
bad_y_vals  =  position_data (iy^ 2 ) ; 
bad_z_vals  =  position_data (iz^ 3) ; 

if  (  (length (ix)  ~=  0)  | |  (length (iy)  ~=0)  | |  (length (iz) 

~=0)  )  %If  there  is  at  least  one  value  outside  of  the  range 
for  c  =  1: length (ix)  %Display  x-values  and  indices  that 
exceed  range 

disp([*Your  x  value  of  *  num2str (bad_x_vals (c) )  * 

inches  in  point  number  *  ... 

num2str (ix (c) )  *  are  out  of  the  * 

num2str (x_limit )  *  inch  range.  Please  revise.*]) 

end 

for  c  =  1: length (iy)  %Display  y-values  and  indices  that 
exceed  range 

disp([*Your  y  value  of  *  num2str (bad_y_vals (c) )  * 

inches  in  point  number  *  ... 

num2str (iy (c) )  *  are  out  of  the  * 

num2str (y_limit )  *  inch  range.  Please  revise.*]) 

end 

for  c  =  1: length (iz)  %Display  z-values  and  indices  that 
exceed  range 

disp([*Your  z  value  of  *  num2str (bad_z_vals (c) )  * 

degrees  in  point  number  *  ... 

num2str (iz (c) )  *  are  out  of  the  * 

num2str ( z_limit )  *  degree  range.  Please  revise.*]) 

end 

error (* Exit ing  m-file*)  %Exit  the  m-file35 
else  %If  no  values  are  outside  the  range 
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disp('All  values  are  within  bounds  of  scanner!')  %Let 
the  user  know  the  scan  will  proceed 
end 

%Convert  position  data  to  step  information 
stepsPerInch  =  520;  %Define  scanner  resolution 
stepsPerDegree  =  67; 

step_data_x  =  round ( stepsPerInch*position_data (:, 1 )  ); 

%Convert  x  direction  to  steps 

step_data_y  =  round ( stepsPerInch*position_data ( : , 2 )  ); 

%Convert  to  steps 

step_data_z  =  round ( stepsPerDegree*position_data ( : , 3 ) ) ; 

%Convert  to  steps 

delay  =  position_data ( : , 4 ) ; 

%Add  in  delay  column 

Arduino_data  =  [step_data_x  step_data_y  step_data_z  delay] ; 
end 
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function  [kx  ky  kz]  =  calc_wavenums ( f reqs,  Ny,  Nx) 


freqs  =  freqs; 
c  =  3e8; 

k  =  2*pi*freqs/c; 

NofFreqs  =  length ( freqs ) ; 
kx  =  zeros (Ny, Nx, NofFreqs ) ; 
ky  =  kx; 
kz  =  kx; 

for  i  =  l:NofFreqs 
for  j  =  1 : Ny 

kx(j,:,i)  =  linspace(- 

end 

for  j  =  1 : Nx 

ky(:,j,i)  =  linspace(- 

end 

kz(:,:,i)  =  sqrt  (4*k  (i)  "^2 
ky ( : ,  : , i)  . ^2)  ; 
end 

end 
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#def ine 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 


xpinClk 
xpinDir  8 
ypinClk  5 
ypinDir  7 
zpinClk  4 
zpinDir  9 
pinStop  2 
pinsensor  10 


int  StepsX; 

int  StepsY; 

int  Steps Z; 

int  PauseBtwPoint s ; 

int  targetStepsX; 

int  targetStepsY; 

int  targetStepsZ; 

int  xPos  =  0; 

int  yPos  =  0; 

int  zPos  =  0; 

int  endStop; 

int  MoveDelay  =  500;  //This  controls  the  speed  of  the 
scanner . 
int  sensor; 


void  setup ( )  { 

pinMode (xpinClk, OUTPUT)  ; 
pinMode (xpinDir, OUTPUT) ; 
pinMode (ypinClk, OUTPUT) ; 
pinMode (ypinDir, OUTPUT)  ; 
pinMode (zpinClk, OUTPUT)  ; 
pinMode (zpinDir, OUTPUT) ; 
pinMode (pinStop, INPUT) ; 
pinMode (pinsensor, INPUT)  ; 
digitalWrite (pinStop,  HIGH)  ; 

Serial . begin ( 9  60  0 )  ; 

Serial . set Timeout ( 10 )  ; 

Serial . flush ( )  ; 

Serial . print In ( digit alRead ( sensor) )  ; 

SetToHome ( ) ; 

Serial_f lush ( )  ;  //  once  home  is  set,  set  the  serial  port 

to  let  matlab  know 

} 


void  loop ( )  { 

//This  is  the  main  code  that  will  run  for  each 
coordinate . 

while (Serial . available ( )  <  1) {  //Wait  for  a  data  point 

to  be  sent  over  serial 


} 

StepsX  =  Serial . parseint () ;  //X  value  in  steps 
StepsY  =  Serial . parseint () ;  //Y  value  in  steps 
StepsZ  =  Serial . parseint () ;  //Z  value  in  steps 
PauseBtwPoint s  =  Serial . parseint  () ;  //Delay  in 
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milliseconds 


if (StepsX  >  xPos) 

digitalWrite (xpinDir,  LOW)  ; 
else 

digitalWrite (xpinDir, HIGH) ; 
if (StepsY  <  yPos) 

digitalWrite (ypinDir, LOW) ; 
else 

digitalWrite (ypinDir, HIGH) ; 

targetStepsX  =  abs (StepsX  -  xPos) ;  //Calculate  target  X 
value  using  current  X  position. 

targetStepsY  =  abs (StepsY  -  yPos) ;  //Calculate  target  Y 
value  using  current  Y  position. 

targetStepsZ  =  abs(StepsZ  -  zPos) ;  //Calculate  target  Z 
value  using  current  Z  position. 

MoveXMot (targetStepsX) ; 

MoveYMot (targetStepsY)  ; 

MoveRotMot (targetStepsZ)  ; 
xPos  =  StepsX; 
yPos  =  StepsY; 
zPos  =  StepsZ; 

Serial_f lush ( ) ; 

} 


void  MoveRotMot ( int  stepsZ) 

{ 

int  k  =  1; 

if ( stepsZ  >  0 ) { 

digitalWrite ( zpinDir,  HIGH)  ; 
for (  k;  k  <= (stepsZ);  k++) 

{ 


digitalWrite ( zpinClk, HIGH)  ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (zpinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 

} 

digitalWrite ( zpinClk, HIGH) ; 

} 

else  { 

if  ( stepsZ  <  0 )  { 

digitalWrite (zpinDir,  LOW)  ; 
for (  k;  k  >=  (stepsZ);  k — ) 

{ 

digitalWrite (zpinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite ( zpinClk,  HIGH)  ; 
delayMicroseconds (MoveDelay)  ; 

} 
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digitalWrite (zpinClk, LOW) ; 


} 


} 


} 


void  Serial_f lush ( ) { 

while (Serial . available () >0 ) { 

Serial . read ( )  ; 

} 

delay (PauseBtwPoints )  ;  //Delay  to  give  positions 
settle . 

Serial . print In ( 1 ) ;  //Print  value  to  line  telling 
continue 
} 

void  SetToHome(){ 

//Set  X  Home 

endStop  =  digitalRead (pinStop) ; 
digitalWrite (xpinDir, HIGH)  ; 
while (endStop  !=  1){ 

endStop  =  digitalRead (pinStop) ; 
digitalWrite (xpinClk,  HIGH)  ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (xpinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 


} 

delay ( 500 ) ; 

digitalWrite (xpinDir, LOW)  ; 

for(int  j  =  l; j<=50;  j++)  { 

digitalWrite (xpinClk, HIGH) ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (xpinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 

} 

delay ( 1000 ) ; 

xPos  =  0; 

//Set  Y  home 

endStop  =  digitalRead (pinStop) ; 

digitalWrite (ypinDir, LOW)  ; 

while (endStop  !=  1){ 

endStop  =  digitalRead (pinStop) ; 
digitalWrite (ypinClk, HIGH)  ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (ypinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 

} 

digitalWrite (ypinDir,  HIGH)  ; 

delay ( 500 ) ; 

for(int  1=1; 1<=50;  1++)  { 

digitalWrite (ypinClk, HIGH) ; 
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delayMicroseconds (MoveDelay)  ; 
digitalWrite (ypinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 

} 

delay ( 100  0 )  ; 
y  P  o  s  =  0 ; 

//Set  Z  Home 
/* 

sensor  =  digitalRead (pinsensor )  ; 
while (sensor  !=  0) 

{ 

digitalWrite ( zpinClk, HIGH) ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (zpinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 
sensor  =  digitalRead (pinsensor )  ; 

} 

delay (2  00 )  ; 

digitalWrite ( zpinDir,  HIGH)  ; 
for(int  j  =  l; j<=450;  j++)  { 

digitalWrite ( zpinClk, HIGH)  ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (zpinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 

} 

zPos=0 ; 
delay ( 500 ) ; 

*/ 

Serial_f lush ( )  ; 

} 

void  e_STOP ( )  { 

digitalWrite (xpinClk, LOW) ; 
digitalWrite (ypinClk, LOW) ; 
while ( 1 ) { 

} 


void  MoveXMot(int  stepsX) 

{ 

endStop  =  digitalRead (pinStop) ; 
if (endStop  ==  1) 
e_STOP ( ) ; 

for (  int  k=0;  k  <=  stepsX;  k++) 

{ 

digitalWrite (xpinClk, HIGH)  ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (xpinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 

} 

} 
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void  MoveYMot(int  stepsY) 

{ 

endStop  =  digitalRead (pinStop) ; 
if (endStop  ==  1) 
e_STOP ( ) ; 

for (  int  k=0;  k  <=  stepsY;  k++) 

{ 

digitalWrite (ypinClk, HIGH)  ; 
delayMicroseconds (MoveDelay)  ; 
digitalWrite (ypinClk, LOW) ; 
delayMicroseconds (MoveDelay)  ; 

} 

} 
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